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Strength of reinforced brickwork is influenced by 
several variables such as size and properties of brick, 
properties of mortar and reinforcement, and proportions of 
the section. A large number of brick samples were studied 
for the statistical variations in tiiair dimensional and 
strength properties. Bimensions of bricks manufa,ctured from 
a single manufacturer ha,ve a very sma,ll variability with a 
tendency towards deteimiinistic value. The coefficient of 
variation in length, breadth and height of bricks was 
limited to 3*4 percent indicating a high consistency in the 
dimensional tolerance of the hand made bricks. Bensity of 
dry bricks of different individual manufacturers follows 
noimal distribution and the coefficient of., variation is 
found to be within 5 percent which again reflects a good 



finality control in brick making. Practically all properties 
such as dimensions, dry density, wet density, percentage of 
water absorption and compressive strength of bricks are 
found to follow normal distribution. Gement~sand motar 
cubes of three different mixes were cast in the laboratory 
and tested after 28 days to study the variability of mortar 
strength. Histograms and cumulative distributions of 
compressive strength of mortar for three different mixes are 
presented. Formal distribution appears to be a good model 
to represent the compressive strength of mortar, Ihe 
statistical variation of joint thickness of mortar in 
different existing buildings was also studied and found to 
follow normal distribution at 1 percent significance level. 

The strength of brick masonry depends on the quality of 
brick, mortar, joint thickness and joint layout. Due to non 
availability of field data, the strength of masonry was 
generated by Monte Carlo simulation on a digital computer 
from the randomly generated dsfa of brick and mortar 
strengths using their statistical properties. Histogram and 
cumulative distribution of simulated masonry strength are 
presented. The central region of the simulated data matched- 
with normal distribution but mismatching was observed at 
either tails. Strength of masonry may still be assumed 
as normally distributed random variable. 
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Chi-square test and Kolmogorov-Smirnov test 
generally used to examine whether the data of a 

particular random variable follovfs an assumed probability 
distribution. Usual Chi-square test with equal class 

interval for testing the suitability of a hypothesised 
probability distribution has limitations. It depends on 
three arbitrary decisions such as the choice of the starting 
point, number of classes and width of classes. Thus it may 
lead to doubtful acceptance or rejection of the null 
hypothesis. The three arbitrary choices can be made to a 
single variable if the test is conducted with classes of 

equal probability. The level of significance at which the 

null hypothesis can be a,ccepted is found to vary raMomly 
T^ith the choice of number of classes. Therefore, a series of 
tests should be conducted with different number of classes 
for a dependable decision. Kolmogorov-Smirnov test ( K- S 
test ) does not involve any arbitrary choice of paramecers 
like that of Ghi-square test. Standard tables available for 
use with K-S test are only valid when the paraneters of the 
null distribution is fully specified. Those tables cannot 
be used when the parameters of the null distribution is not 
specified but have to be estimated from the data. The 
critical values of K-S statistic for different level of 
significance are generated by Monte Carlo simulation and are 
presented for use vnth K-S test. 



Reliability analysis of reinforced brick beam(RBB) 
under flexure is presented^ Due to random variations in 
strengths of masonry and steelj* there exists a probability 
that a RBB section will fail as an over-reinforced sven 
though the section is designed as an under-reinforced based 
on deterministic analysis, A general method of analysis of 
probability of failure of RBB section is presented treating 
strengths of masonry and steel as random variables. 
Probability of failure (p^) of a RBB section is the sim of 
the probabilities of the section failing as an 
under-reinforced (pf^.^ failing as an over - reinforced 
(p^o^* Ihe computation of the above probabilities involves 
evaluation of multiple integrals. As the area of steel In a 
RBB section is increased, the probability of failing as an 
under-reinforced (Pf^) decreases and the probability of 
failing as an over-reinforced (pfo) increases but the 
probability of failure (p^) decreases for the same external 
moment, Monte Carlo simulation is also used to study the 
reliability of the section, 

A method of reliability/- based design of RBB section is 
presented. The method involves solution of an integral 
equation to find the area of steel or the dimensions-, of the- 
section for a preassigned reliability. Examples are given to 
illustrate the method, limit state design method based on 
semi -probabilistic approach is illustrated considering the 
coefficient of variation of resistance and load. Partial 



safety factors applied to strength of materials and load for 
a target reliability depend on the coefficient of variation 
of resistance and load» A set of partial safety factors is 
suggested for limit design of RBB sections. 

Aspects of durability of reinforced brickworlc are 
reviewed and discussed. Corrosion of reinforcement in 
reinforced brickwork appears to be a serious problem. 
Estimation of time of cracking due to corrosion of 
reinforcement is formulated, A method of reliability 
analysis, taking into account the decrease of moment capacity 
with time, is presented. Ihe load factors or partiaj. safety 
factors needed for the design of a HBB section for a given 
design life depends on the rate of deterioration and type of 
exposure. This is illustrated through an example. Conclusions 
and recommendations based on the discussions of the results 
are presented at the end. 
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nSKRODUCTIOlT 

1,1 Introduction 

Brick masonrj -construction is one of the oldest and it 
continues to dominate in residential building coiietruction. 
Brick units bonded together by mortar in a predetermined 
orientation used for resisting external loads is called 
brick masonry. Burnt clay building bricks are produced in 
India for building construction. The concept of reinforced 
masoncy structural member to resist tensile forces is not 
new to the Civil Engineering profession. Mild steel flats 
and rounds have been used for many years in some form or 
other. The use of high yield strength deformed (HYSD) bara 
in reinforced brick masonry has become more popular for its 
improved strength and bond capacity, 

Structuiral. safety is provided by means of factor of 
safety in the conventional working stress design philosophy. 

The safety of a structural member against external loads is 
ensured by limiting the actual stresses to a set of permissible 
stresses whereas the safety is ensured through proper choice 
of load factors in ultimate strength design. In both these 
design philosophies load and strength of materials are treated 
as deterministic. However, the real loads and strength of 
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niaterials are probabilistic rather than deterministic. 

To provide a meaningful safety to a structure, a semi- 
probabilistic design philosophy using different limit 
state has come up taking into consideration the variability 
in material strength and load. 

1 .2 literature Review 
1 .2*1 Brick masonry 

The general development of reinforced brick masoncy 
was reviewed by Srivastava (l )* and Dayaratnam et al. (2). 
Poster ( 3 ) has given an extensive review of the development 
of reinforced and prestressed brick work up to 1981, For 
the sake of completeness, the main developments sire reviewed 

I 

here. Some of the investigations, even though do not deal 
with RBM, are reviewed to establish certain relevant strength 
and stress relations. 

It was reported by Beamish ( 4 ) and Plummer et al.(5) 
that Sir M.I, Brunei is the first one to use the concept of 
reinforced brick masonry in building a chimney in 1813# 
Results of tests on RBM structural moabers was reported by 
Brebner (6), Brebner's work is believed to be the first 
organised research work carried out in India around 1920* 
Several tests on Reinforced Masonry beams were carried out 
by Parsons et al. (7) and Withey (8) during 1932-33* In 
1932 , the first draft code of the recommended practice and 

♦ aihe numerals in parenthesis indicate the reference numbers 
given at the end* 
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specification for RBM construction was developed (5)* Since 
then, several investigations were conducted on plain and 
reinforced brick masonry at different places. 

Different codes of practice (9»10,1l) outlined the 
methods of designing plain and reinforced masonry by permissible 
stress approach, Basic compressive stress of brickwork for 
different brick strength and different mortar type is given 
in tabular form. Strength of brickwork in compression can 
be calculated through basic compressive stress by allowing 
different reduction factors for slenderness ratio, shape of 
individual units and eccentricity, 

’Building Code Requirements for Reinforced Brick 
Masonry’ (12) published in i960, and ’Recommended Practice 
for Engineered Brick Masonry’ (I 3 ) published in 1969 are 
the outgrowth of different research works carried out in 
U,S,ji, on plain and reinforced brick masonry. Cutler 
et al, ( 14 ) explained how research informations reported by 
different countries and individuals have been taken into 
consideration for developing the Canadian Building Code for 
masonry. The structural design requirements for load bearing 
nonreinforced brick masonry walls contained in various 
standards and codes were compared by Cross et al, (15). 
Comparison of different standards and codes are difficult 
because of difference in test methods and requirements of 
brick and mortar, size and shape of the individual units 
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as well as the different approaches used to consider the 
various design factors. 

Hallar (16) reported test results of shear test on 
brick masonry. Several researches were carried out to study 
the ultimate shear strength of reinforced brickwork beams. 
Suter and Hendry (1?) concluded from experimental results, 

*the ultimate shear stress of reinforced brickwork beams 
increases only slightly with increasing amount of tensile 
reinforcements but increases significantly with decreasing 
shear-span to effective depth ratio ' , It is recommended 
that the significant increase in ultimate shear stress with 
decreasing shear span to effective depth ratio should be 
recognised in limit state shear provisions. Ihe effect of 
ratio of shear span to effective depth, the amount of 
tensile reinforcements and the compressive strength of 
brickwork on the ultimate shear resistance of RB beams was 
examined and limit state shear values were proposed (18), 

It has been observed that the failure of masonry, 
loaded in compression is initiated by vertical splitting of 
bricks, Veirious attempts have been made to relate the 
compressive strength of masonry to the tensile strength of 
bricks (19, 20). Stafford Smith (21) discussed existing method 
of testing brick masonry discs to establish the diagonal 
tensile strength of brickwork. It was shown through 
experimentaO. results and finite element analysis that the 
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diagonal tensile strength of bricicwork is approximately 
equal to either the tensile strength of mortar or the 
tensile strength of brick, -whichever is smaller (21). 

Based on the data reported by Talbot and Abrams (22), 
Parsons ( 23 ) developed an empirical formula for estimating 
the probable compressive strength of concentrically loaded 
masonry walls or columns constructed with hollow units, Oihe 
foimula relating strength of masonry with strength of brick 
and mortar can be stated as(23)t 


where f^ 

^b 

^m 

a33d 3c^ 


ivircfbf^) ( 1 . 1 ) 

strength of masoniy, 
strength of brick, 
strength of mortar, 

a factor which depends on the orientation of 
the individual units. 


Dayaratnam et al.(2) carried out regression analysis to 
fit Eq., 1,1 to the experimental results and recommended 
different range of values for k for calculating the strength 
of masonry, A formula for predicting strength of masonry 
from strength of brick and tjrpe of mortar is recommended by 
SC5PI(13). Hoath et al. ( 24 ) have shown that brickwork cube 
made with the traditional cement -lime-sand mortar have 
higher ratio of brickwoiii cube strength to mortar strength 
as compared to that of cement-sand mortar, Prancis 
et al, ( 25 )have studied the effect of joint thickness on the 
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compressive strength of brickwork. Experimental results 
showed that pi^ism compressive strength decreases as the 
average joint thickness increases. 

‘Ihe problem of masonry infilled frame was dealt by 
Garter and Stafford Smith(26). !Ehe mode of failure was 
found to be governed by the geometry of the structure and 
also by the relative values of the bond shear strength, 
internal friction and diagonal tensile strength of masonry. 
Stafford Smith et al,(27) has studied the distribution of stresses 
in brickwork walls subjected to vertical loads through 
finite element approach. It was found that the peak tensile 
stress values in a wall increases with height to length 
ratio of the wall and with the brick to moitar elastic 
modular ratio. When a structxiral frame of building is 
infilled with masonry, the diaphragm action of the wall in 
the frame produces substantial increase in the racking 
stiffness of the structure(26,27). Riddington et al.(28) 
have shown by finite elCTient approach that the failure loads , 

Iflie magnitude of diagonal tension and shear stresses are 
governed almost by the length to height ratio of the infill. 

The composite action of masonry walls with beams and frames 
has been studied by Stafford Smith(29). Design methods to 
take care of the composite action of wall with beams and 
frames were also discussed in the paper. 
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The compressive strength of a wall is related to 
compressive strength of brick and mortar, but it also 
depends on the height of the bricks used(30). Howard 
et al, (31 ) have pointed out the effect of manufacturing 
variables on the strength and durability of bricks. The 
effect of construction variables on the strength of brick 
masonry was also discussed. Grimm and Halsell(32) 
illustrated the difference between good and bad workmanship, 
and discussed quality control in brick masonry construction 
as. practiced in U.S.A. Thomas(53) discussed about the 
quality control and suggested regular testing of bricks and 
brickwork test cubes. ITariables affecting strength of brick 
masonry was reviewed and summarized by Grimm{34)« It was 
concluded that compressive strength of brick mason 2 :y prism 
is a function of compressive strengths of brick and mortar, 
workmanship quality, magnitude and direction of load, 
eccentricity and prism geometry, age and curing. Empirical 
equations relating strength of brick masonry with different 
variables are derived based on different research data. 
Based on the data obtained by different researchers, 
Plummer (35) commented, 'The wide variations between 
individual specimens of the same group frequently 
encountered in such data, are due to large number of 
variables, some of which are very difficult to control, to 
which all tests, which might be e3cpected to produce results 

comparable to those obtained in masonry construction, are 
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sensible. These Tariables include air content and flow of 
mortar, elapsed time between spreading mortar and placing 
brick in contact with it, suction of brick, pressure or 
tapping applied to' joint during forming, texture of brick 
surface and probably, to a minor degree at least, other 
factors which have not been identified , Hendry ( 36 ) 
critically discussed the strength and properties of load 
bearing briclavork materials with particular reference to 
workmanship. 

Based on statistical properties of masonry, Maochi(37) 
discussed about the choice of suitable factors of safety 
which also take into account the structural behaviour, 
Hendiy(38) pointed out, 'Research work which results in the 
reduction of the coefficient of variation by even a few 
percent is likely to be of great practical significance in 
terms of reducing the safety factor and thus achieving 
efficient use of materials', G-raphs of central safety 
factor needed for different variability of material strength 
when designed for a particular probability of failure were 
given (37,38). 

With a view to incorporate separate safety provisions 
to loads and strength of materials as a more rational approach 
to the design problem, the limit state format has nov/ been • 
adopted in different national codes (39,40,41) for design of 
reinforced concrete, Czechoslovakia has already developed 
their brickwork code in limit state format ( 42 ). 

After the adoption of limit state format of OP 110: 

1972 ( 39 ), the brick masonry code has been revised and 
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developed in the same line of reinforced concrete, and as 
a result BS 5628; Part 1(43) for plain brick masoniy is 
published. At the same time, a committee of engineers and 
architects produced SP 91 (44) to serve as a guide to the 
design of reinforced and prostressed masonry. 

Stress-strain curves of brick masonry for different 
cement-sand mortar are presented by Jain et al. (45). Based 
on experimental investigation, stress - strain relation of 
brick masonry prisms is derived and presented by Dayaratnam 
et al, (2). Ultimate load theory similar to reinforced 
concrete is found to be satisfactory in determining moment 
capacity of RB beams and slabs (1,2,45). Beard (46) 
suggested a formula for predicting ultimate strength of 
reinforced brick beam in flexure using a parabolic stress 
strain relationship with a falling branch and rejected the 
enhancement factor included in SP 91 (44). Sinha (47) 
adopted another approach of predicting the ultimate strength 
using a curvilinear and cubic parabolic stress block. 

Recently, Britij^ Standard Institution (331) has published a 
draft code BS 5628 ; Part 2(48) for reinforced and prestressed 
brick masonry in a limit state format, A simplified rectangular 
stress block is adopted in the draft code. Partial safety 
feu;tors 2,5 and 2.8 for brickyrork given in the draft code 
was criticized by Beard (49). Various clanses in the draft 
code (48 ) relating to the shear strength of reinforced brickwork 
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are discussed by Hendiy(50). 

1*2,2 Probabilistic concept of safety 

The basic concept of structiiral safety analysts was 
first discussed by Preudenthal(51 ) and later by Asplund.(52) 
and Pugsley(53)* Preudenthal(51 ) said, 'With, increasing 
perfection of desi^ methods, the element of 'ignorance* 
can be largely eliminated; but the element of 'urwertainity' 
is caused by circumstances that cain be changed, to a certain, 
extent, but can never be removed* Hence, the safety factor 
is a measure of imcertainity rather than ignorance' * The 
statistical frequency distribution of load and resistance 
must be considered in determining the reliability, 
5'reudenthal(54-) commented, ’ 'Because the design of a 
structure embodies uncertain predictions of the structural 
materials as well as the expected load patterns and 
intensities, the concept of probability must form an 
integral part of any rational design; any conceivable 
condition is necessarily associated with a numerical measure 
of the probability of its occurrence* It is by this measure 
alone that the structural significance of 'a specified 
condition can be evaluated ' ' * 

The concept of unserviceability and failure was 
defined anfl procedur®for computation of probability of 
failure were outlined by i'reudenthal(54)* Statistical data 
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on intermediate grade reinforcing steel and ultimate strengtli 
of concrete v/as presented by Julian(55). G-raplis of minimum 
required factor of safety versus probability of failure for 
normal and lognormal distributiors f or load and resistance 
were sbown(55). Methods of providing engineering safety and 
additional safety necessary for social purposes were 
formulated by Brown(56), Attempts for providir^ safety can 
be made by specifying (i) a minimum ratio of resistance to 
maximum design load,(ii) a maximum value of probability of 
failure and (iii) a minimum factor of safety coupled with a 
maximum probability of failure (56). Ihe concept of structural 
life along with functional and collapse failures was aiso 
discussed. 

Methods of safety analysis based on lognormal and 
extremal distributions of resistance and of load were 
developed and the effects of redundancy on the risk of failure 
and reliability function were examined(57). Ihe work of 
the Task Committee on factors of safety, ASCE, were summarized 
by Ereudenthal, G-arelts and Shinozuka(58). Probability of 
failure (p^) has been recognised as the probability of the 
occurence of the event (R-S) < 0 or (R/s) < 1 and can be 
computed as (for R and S independent) 

p^ = P(R < S) = // %(r) fs(s) dr ds (1,2) 

in which I) is the region of integration defined by the event 
(R < S), probability density functions 
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of resistance R and load S respectively. Reliability 
functions were derived for the cases when there is N 
repeated load application at equaiL inteirvals or at 
prescribed instance and the number of occurence of load is 
governed by the Poisson’s law, Numerical examples of trussed 
tower and aircraft agad-nst wind load were given to illustrate 
the method. Derivation of reliability functions was also 
extended for structures with multiple members. Bounds on 
reliability of structural system were presented by Cornell (59), 
Reliability theory of deteriorating structures was formulated 
by Kameda and Koike (60) considering resistance deterioration 
under repeated ramdom loads. 

An extensive review of the available literature till 
December, 1971 on structural saifety v/as presented by the 
Task Committee on Structural Safety (61). The choice of 
safety level in structural design is one of the most difficult 
job. The choice of safety level based on minimum expected 
loss criterion v/as formulated by Turkstra(62), Ang and 
Amin(63) formulated a method for the determination of 
reliability of a statically determinate system amd established 
a monotonically decreasing property for the hazard function. 

The method waa also extended for indeterminate system a n d 
numerical examples were given, Ang and Amin(64) proposed 
a new concept of structural safety where the lack of 
knowledge aufi information was distinguished from statistical 
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dispersions, The lack of knowledge and information was 
handled through a judgement factor whereas the statistical 
variables were treated probabilistically. Reasons and 
advantages of probabilistic studies of safety were discussed 
and the basic probabilistic procedures in the analysis of 
structural safety were introduced by Shah(65) using only 
mean values and standard deviations of resistance and load 
for an assumed probability distribation of (R-S) or (E/s), 

Sexanith and Nelson (66) discussed several limitations 
in the application of probabilistic concepts. The major 
difficulties were recognised as (i) the choice of the 
appropriate probability model, (ii) introduction of 
subjective elements into the probabilistic structure and 
(iii) interpretation of probabilistic information in a form 
which leads to rational decisions. Different techniques to 
avoid these difficulties were presented by other^67,68). 

The approximate probability measures associated with the 
existing code provisions was investigated by Benjamin and 
Dind(67) and reliability based procedures were illustrated 
using a set of assigned probabilities. The concept given in 
the paper is to relate probability measures to existing 
deterministic code provisions thereby allowing improvement 
of design procedures within the existing deterministic Qode 
format. It was shown that the design is relatively insensitive 
to these probabilities so that a code making body has 
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considerable freedom in assigning such measures of 
probabilities, A consistent first oMer code format based 
on first and second moments of all stochastic variables was 
presented by Cornell (68). Biis approach was extended to 
code formats of arbitrary higher order than first and second 
moments of stochastic variables by Lind(69) and the 
calibration of a partial safety factor format to Cornell 's(68) 
format as well as to Ang— imin(64) format was demonstrated* 

Systematic analysis of uncertainity through first 
order statistical analysis and the explicit use of probability 
of failure or survival as a measure of risk and safety were 
suggested by At^CTO), Practical methods of risk assessment 
in terms of failure probability and for the development of 
reliability based design criteria were described and 
developed, Hasof er and Iiind(71 ) presented an exact and 
invariant code format and described a natural measure of the 
second-moment reliability index of a design as the distance 
from the mean of basic variable vector to the failuxe region 
boundary when all the variables are measured in standard 
deviation units. It was shown that this measure does not 
depend in any way on the precise analyiiical form of the 
failure criterion. 

Systematic evaluation of risk and analysis of 
uncertainity were developed by Ang and Cornell (72). Under 
conditions of uncertainity, safety and serviceabili:fcy of 
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structures can be assured only in terms of the probability 
of survival (or conversely of failure). Uncertainity 
(including uncertsiinity associated with errors in estimation 
and imperfection in mathematical models) was expressed in 
terms of the coefficient of variation, formulation for 
multiple load cases was also presented. It was emphasized 
that the distribution of safety margin (R-S) or ln(R/s) is 
important in the calculation of probability of failure and 
not the distribation of the individual variates (7 0,72), At 
high risk levels, e.g. , p^ 2 10”^, the calculated failure 
probability p^ is not very much different regardless of the 
t 3 rpe of the distribution,however it depends significantly on 
the choice of the distribution of (R-S) or ln(R/s) for a low 
risk level (pf < 10”^). It was argued (70,72) that the 
relative measure of risk is all that is necessary and for this 
purpose a prescribed distribution of (R-s) or ln(R/s) should 
suffice. A distribution function? composed of a central 
gaussian portion and an exponential tail with a single point 
of discontinuity, was proposed by Iiind(73) as a convenient 
representation of statistical demand and capacity data to 
evaluate the reliability, regardless of any conjecture about 
the parent distribution. The convolution integral of the 
probability of failure was evaluated in closed form using 
simple statistics of the data, Sllingwood and Aug(74-) 
formulated a risk based design criteria and showed how the 
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level of risk is affected by the uneertainities of the 
design variables, fhe analysis of design -uneertainities 
arising from the statistical nat-ure of some variables and 
incomplete information regarding others was examined. 

Ravindra et al,(75) illustrated the design method proposed 
by Ang and Cornell (72) with practical examples and showed 
that the code parameters such as safe-ty indices can be 
selected to match the safety level of current designs. Only 
mean values and standard de-vLations were used to represent 
the design random variables and the final design criteria 
were developed on a probability-distribution free basis. 

Models of structtiral reliability of systems including 
determinate structures (weakest link system), indeterminate 
structures (fail safe system) and brittle members were 
described by Moses(76), TJsing the design format suggested 
by Ang and 0omell(72), Moses ( 76 ) introduced partial Scifety 
factor defined as the quantity by which the element safety 
factor should be changed to account for its presence in a 
structural assonblage so that the overall failure 
probability remains roughly equal to what is desired for 
the element, 

Costello and Chu(77) presented a formulation for 
calculation of failure probabilities of reinforced concrete 
rectangular beams v/ith and without compression reinforcoments 
assuming beta distribution for strengths of concrete and 
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steel, Warner and Kabaila(78) described a Monte Carlo 
simulation technique to determine the eunrulati-ve distribution 
functions of stochastic variables which are functions of 
several stochastic variables, The cumulative distribution 
function of ultimate strength of an axially loaded short 
reinforced concrete colunn was investigated using this 
technique and the results were compared with closed form 
solution. Procedures for determining the probability 
distributions of the safety margin and related quantities 
of several types of reinforced concrete members were 
presented by Sexsmith (79), Using Monte Carlo technique, 
Allen (80) made a detailed study on estimated probability 
distributions for the ultimate bending strength and 
ductility ratio (ratio of curvature at ultimate to curvature 
at yield) of rectangular concrete sections reinforced in 
tension only. It was concluded that there is a significant 
probability that an under-reinforced section designed 
according to ACI 318-65(81) wotild undergo a brittle 
compressive fail\zre. The variability of ductility ratio was 
found, much higher than that of ultimate bending strength. 
Reliability analysis of prestressed concrete simply 
supported beams designed on the basis of Indian Standards, 
American Concrete Institute and Recommended Code of Practice 
by GEB-PIP(82) was presented by Chandrasekar and 
Dayaratnam (83,. 84), Ranganathan and Dayaratnam (85,86 )used 
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Monte Carlo Sinnilation technique to evaluate the probabi3J.ty 
of failure of a prestressed concrete flanged section. The 
failure of the section at limit state of strength was divided 
into two cases ;(i) under-reinforced and (ii) over-reinforced 
and the probability of failure was evaluated based on the 
above two conditional failxire events. 

Moses and Kinser (87) discussed on optimum sizing of 
elCTients for multielement and multiload conditions with 
safety defined in terns of allowable probability of failure. 
Some examples (trusses and frames) of the optimal 
reliability based design were presented by Moses and 
Stevenson (88), Rao (89) demonstrated the feasibility of 
optimum cost design of under-reinforced concrete beams with 
a reliability based constraint. The design was considered 
to be safe vhen the expected tiltimate moment capacity 
exceeds the design moment by a certain number of standard 
deviations. Optimization of prestressed concrete beams 
subject to reliability constraint in addition to the normal 
constraints was presented by Ghandrasekar and Dayaratnam(90), 

1 ,3 Statement of the Problem 

The design of reinforced brick masonry in India is 
being carried out by working stress method. The safety of 
the structure is ensured by limiting the actual stresses 
less than or equal to permissible stresses. A semi- 
probabilistic or probabilistic design method which 
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incorporates the uiKjertainity in strengths of materials and 
load is more i^tional, k semi-probabilistic or probabilistic 
design criteria similar to that of reinforced concrete is to 
be developed to provide a rational safety margin in the 
design of reinforced brick masonry* The object of the 
present investigation is as follows: 

(a) Statistical analysis of physical and some 
impoirtant mechanical piaperties of brick and 
mortar* 

(b) Prediction of empirical probability laws of strengths 
of brick, mortar and masonry* 

(c) Presentation of a method of reliability analysis 
of reinforced brick beam in flexure* 

(d) Design of reinforced brick beams for a given 
probability of failiire, 

(e) Study of the variation of load factors with 
coefficient of variation- of resistance of the 
section and probability of failure* 

(f ) Durability study and analysis of probability of 
life expectancy of reinforced brick beams. 

The organization of the thesis is described in the 
following paragraphs* The general introduction of reinforced 
brickwork along with the relevant literatures on the general 
development of reinforced brickwork and probabilistic concept 
of saf ety are presented in CJhapter 1 • 
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Ciiapter 2 deals with, the statistical analysis of field di 
on brick and mortar* Dimensional and strength properties 
of bricks were measured in the laboratory and statistical 
analysis was carried out. Cement-sand mortar cubes of three 
different mixes were cast in the laboratory and tested 
after 28 days of curing. Statistical analysis of strength 
of mortar was also carried out to study the variability of 
mortar strength. Due to nonavailability of field data, 
strength of masonry was generated by Monte Carlo simulation 
using the deterministic formula relating the strengths of 
brick, mortar and masonry. Samples of masonry strength were 
generated from the specified distribution of strengths of 
brick and mortar and the statistical analysis was carried 
out. Chi-square test was conducted to test the suitability 
of normal and lognormal distributions to the data of 
different dimensional and strength properties. Finally, 
relation between characteristic and mean strength of masonry 
was established assuming normal and lognormal model for 
masonry strength. 

Reliability of chi-square test is discussed in 
Chapter 3. Usual chi-square test depends on three arbitrary 
decisions such as the choice of the starting point, number 
of classes and width of classes selected for the test and 
thus may lead to doubtful results. Some of these problems 
are illustrated through examples and a practical procedure 
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is suggested for a dependable decision* KolmogoroT-SniirnoT 
test (K-S test) does not involve any arbitrary cboice of 
parameters liJce that of chi-square test. Standard tables 
available for use with the Kolmogorov-Smirnov test are only 
valid ^en the null distribution is completely specified, 

The critical values of K-S statistic for different level of 
significance are generated by Monte Carlo Sirrulation for use 
with Kolmogorov-^mimov test when the parameters of the null 
distribution is not specified but have to be estimated from 
the data. 

Chapter 4 of the thesis deals with the reliability 
analysis of reinforced brick beam (EBB) at flexural 
strength. !Phe failure of a reiiaforced brick beam section is 
divided into two cases (i) under-reinforced and (ii)over- 
reinforced. A general method of analysis of probability 
of failure of EBB section is presented based on the above 
two failure events. Monte Carlo simulation is also used 
to study the reliability of the section and a comparative 
study is presented. 

A method of reliability based design of EBB section is 
presented in Chapter 5. The method involves solution of 
an integral equation to find the area of steel or dimensional 
properties of the section. Examples are given to illustrate 
the method. Design method based on semi-probabilistic approach 
is illustrated considering the coefficients of variation of 
resistance and load. 
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Aspects of durability of reinforced brickwork are 
reviewed and discussed in Gbapter 6 with particular 
reference to corrosion of reinforcement, Sstimation of 
time to cracking due to coirrosion of reinforcement is 
formulated. A method of reliability analysis and design 
considering strength deterioration of RB beam, is also 
presented. Conclusions and recommendations dra'vm based 
on the discussions of the results are presented at the end. 



GHAPIER 2 


SIATISTIGAL ANlIiYSIS OP PROPERTIES 
OP BRICK, MORTiR AKD MASOmi 

2,1 Introduction 

Brick and brick construction has been used in 
building activities for centuries and it continues to 
dominate among tiie commonly used building materials. 
Reinforced brick masonry coi^truction has got popularity 
in ttiis century because of its low cost compared to 
reinforced concrete. The design of reinforced brick 
masonry requires a knowledge of the properties of brick, 
mortar and masonry strength. The safety against load is 
ensured and reflected by permissible stresses in working 
stress design and by use of load factors and partial safety 
factors in ultimate strength design or limit state design 
philosophy. In the deterministic methods of design, the 
design variables are treated as deterministic even though 
there is a certain degree of variability in them. On the 
other hand, the semiprobabilistic methods of designs 
incorporate the statistical variations in some of the 
design variables. The statistical variations- in brickwork 
depend very much on the constructional practice and the 
degree of quality control. 
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The production of brick and brick masonry 
construction in India is highly labour intensive. The 
production of bricks in India is rural oriented by and 
large, and it appears to use hand moulding methods vdiich 
have been in practice for centuries. Statistical analysis 
of different parameters such as dimensional and strength 
properties of brick and brickwork is presented in this 
chapter. 

2.2 Statistical Analysis of Bricks 
2,2.1 Sample collection and notations 

For statistical analysis of brick properties, 
about 100 samples of bricks were collected from each of 
eight different manufacturers. Three mamifacturers were 
selected from Kanpur zone and the rest five were from 
Faroke zone in Kerala State, The set of bricks from a 
particular manufacturer collected at a time will be called 
as a fixed lot in subsequent discussions. About 100 
samples were selected at random from each manufacturer as 
the representative of different fixed lots corresponding 
notations used in the subsequent sections are aa follow^s: 
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2one Notation* Remarks 


Kanpur 


BK1 

Kanpur 


BK2 

Kanpur 


BK3 

Faroke 

(Kerala) 

BK4 

Faroke 

(Kerala) 

BK5 

Faroke 

(Kerala) 

BK6 

Faroke 

(Kerala) 

BK7 

Faroke 

(Kerala) 

BK8 


Fixed lot 


* Name of the manufacturer is not indicated here, 

A mixed lot is a combination of different fixed 
lots. Data of fixed lots of bricks of Kanpur zone was 
combined and called mixed lot 1 . Similarly data of 
Faroke zone was also combined as mixed lot 2 and the 
statistical analysis is carried out, 

The ‘Structural Engineering laboratory of I.I.T. 
Kanpur has been associated with the testing of bricks 
supplied by different manufacturers. Brick samples of 
about 6 bricks each called as set were supplied by 
different manufacturers at different time to the laboratory 
for testing. The data of those test results were taken 
for statistical analysis and will be called a random lot. 
The random lot consists of 19 sets (total of 115 samples) 
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and can be considered as the representative sample of 
random collection of bricks from different manufacturers. 
Ihe summary of mixed lot and random lot is given below: 


lot type 

Zone 

Rotation 

Remarks 

Combination of 

BK1 , BK2 and BK3 

Kanpiir 

BK9 

Mixed lot 1 

Combination of Bi;4, 

BE:5,BK6,BK7 and BK8 

Paroke 

BK10 

Mixed lot 2 

Random collection 

Uttar 

Pradesh 

BK11 

Random lot 


2,2.2 Parameters of interest 

Ihe parameters which were measiired directly firom the 
brick specimens are called measured parameters, fhe following 
parameters were chosen to study their statistical variations : 


i) 

length of brick 

= 1 (in mm) 

ii) 

Breadth of brick 

= B (in mm) 

iii) 

Height of brick 

= H (in mm) 

iv) 

Dry weight of brick 

= W^(in kg. ) 

t) 

Saturated weight of 

brick = Wg(in kg.) 

Vi) 

Compressive load 

= P (in kg. ) 
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From the above measured parameters, the following 
parameters were computed: 


i) 

Area (flatwise) 

= A (in cm^) 

ii) 

Volume 

= Y (in cm^) 

iii) 

Dry density 

II 

iv) 

Wet or saturated density 

= r^iln kg/m^) 

v) 

Percent of water 
absorption 

= p (percent) 

Ti) 

Compressive strength of 
brick (flatwise) 

= f^(in ET/mm^) 


The measuring and computing proeedxires of different 
parameters are discussed briefly in the following section, 

2,2,3 Measuring procedures 

Iiength, breadth and height of a brick specimen were 
measured up to 1 mm accuracy by slide calipers. Each 
dimension was measured at four different sections and the 
mean of the four was taken as the value of the parameter 
for that specimen. Each brick specimen v/as dried in a 
ventilated oven at a temperature of 110°0 for 24 hours and 
then cooled to room temperature. Dry weight (W^) v/as t.aken. 
The specimen was then immersed completely in clean water at 
room temperature (27 + 2*^0) for 24 hours to saturate the 
specimen. After 24 hours, the specimen v/as removed and all 
excess surface water was wiped out by a damp cloth. The 
wet wei^t (Wg) was taken and the percent of water absorption 
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was calculated by the following formtila (91 ) 

- % 

p (percent) = x 100 (2.1) 

Compressive strength test v/as done according to the 
specification given in Indian Standard Code (92), Progs 
and all voids were filled with 1 ;1 cement sand mortar and 
the sample was cured under damp jute bags for 24 hours 
and then tinder water for the following 3 days. After 3 days 
immersion in water, the specimen was removed and traces of 
moisture was wiped out, -The specimen was then placed on a 
compression testing machine with flat faces horizontal and 
mortar filled face upwards. Maximum compressive load (P) 
was recorded at which the specimen failed. Flatwise 
compressive strength is then calculated as 


P 



where A = flatwise cross-sectional area of the specimen 
- IB. 

2,2,4 Statistical analysis and histograms 


Mean and standard deviation are the two most important 
statistics of a random variable. If ^ 

n different observations of a random variable X, then 


1 ^ 

Z 


n 

3a (x..-X_)^ 1/2 

(S.D.) sx =[.S — — ] 

1=1 n-1 


( 2 . 2 ) 


(2.3) 
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!I!he above formula gives an unbiased estimate of standard 
deviation and the co-efficient of»variation (dj^) is 
given by 

(C.0.Y.)6^ = — (2.4) 

The variability of a random variable is reflected through 
the coefficient of variation. Skewness and kurtosis, are 
calculated as 

81 =^ (2.5) 

^ = l ( 2 . 6 ) 


Coefficient of skewness and coefficient of kurtosis are 
computed as follows 


Cs 


gf / s| 


(2.7) 



§2 / 


(2.8) 

where g^ 


skewness of 

2 


S2 


kurtosis of 

Z 


^s 

== 

coefficient 

of skewness 


Cjj. 

=: 

coefficient 

of kurtosis 



The approximate frequency distribution of a random variable 
is represented by histogram. The approximate number of 
intervals for drawing the Histograms is selected as (93) 

k = 1 + 3.3 log^Q(n) (2.9) 



30 


where k — the number of class intervals 
and n = the total number of samples* 

The statistical analysis of length (It), breadth(B), 
height (H), area (a), volume (v), dry density (r^), wet 
density percent of water absorption (p %) and flatv/ise 

compressive strength (f^) of the brick data for eight 
different fixed lots are given in Table 2,1 throng Table 2.8. 
The statistical analysis of mixed lots BK9 and BK10 are given 
in Table 2.9 and Table 2,10 respectively. The statistical 
analysis of random lot BK1 1 is given in Table 2,11. Units 
of the parameters are not indicated in the tables and should 
be read as given in section 2,2.2. 

Histograms and the cumulative probability distribution 

of I, 3, H, A, V, p ( %) and f^, of fixed lots BK1 

and BK8 are shown in Big. 2,1 through Big. 2.18, The variation 

of mean value and coefficient of variation of compressive 

strength of brick of different sets of random lot BK1 1 is 
shown in Big, 2,19. Big. 2.20 shows the histogram of the random 

lot BK11, A summary of statistical analysis of t?/o fixed lots » 
two mixed lots and random lot for parameters I, K,p {% ) and 
f^ is also given in Table 2,12 for ready reference, 

2,2,5 Chi-square test 

Histogram of the data represents approximate frequency 
distribution. Several nonparametric statistical tests are 
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lable 2.1 ; Statistical Analysis of Fixed lot BKI 
n = 1CX) 


Parameter 

Me ail 

o(;^ ) 



Range 

Bower Upper 

L 

230.7 

1.2 

0.061 

2.403 

224,7 

236.5 

B 

110.8 

2.0 

-0.127 

3.028 

104.0 

116.0 

H 

64.5 

2.4 

-0.283 

2.647 

60.2 

68.0 

A 

255.6 

2.8 

0.056 

3.216 

236.1 

273.2 

T 

1649.8 

4.3 

-0.201 

2.230 

1483.8 

1784.5 


1579.2 

3.1 

0.731 

3.546 

1485.1 

1747.8 


1790.3 

2.4 

0.684 

4.980 

1 691 . 8 

1968.7 

P 

13.41 

14.7 

-0.146 

3.040 

8.2 

18,43 

fb 

23.69 

23.2 

0.379 

2.767 

13.92 

39.35 


Table 2.2 

: St at i st i c al Analyst s 
n = 95 

of Fixed 

lot BK2 


Parameter 

Mean 

&{% ) 

Cs 


Range 

Lower Upper 

L 

227.5 

1.5 

-0,596 

3.586 

215.5 

234.8 

B 

108.1 

2.1 

-0.298 

2.591 

102.3 

112.5 

H 

62.1 

2.8 

0.440 

2.641 

58.3 

66.5 

A 

245.8 

2.9 

-0.611 

4.067 

220.4 

260.1 

Y 

1526.7 

3.9 

0.077 

3.080 

1357.6 

1681.9 


1717.1 

4.9 

0.782 

4.002 

1558.9 

2023.4 

7 

w 

1947.4 

5.9 

0.359 

3.469 

1777.7 

2202.4 

p 

13.47 

14.2 

-0.715 

4.277 

7.53 

18.19 

fb 

21,49 

24* 6 

0.809 

4.141 

8.26 

39.69 


Table 2,3 
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: Statistical Analysis of Fixed lot BK3 

n = 100 


Parameter 

ifean 

6(fO 

Cs 


lower 

Range 

Upper 

I 

227.3 

1.5 

-0.309 

2,226 

219.8 

234.3 

B 

109.4 

1.9 

-0.749 

4.121 

101 .8 

114.0 

H 

63.3 

3.4 

0.047 

3.041 

57 .3 

68.8 

A 

248.7 

3.0 

-0.546 

3.053 

226.1 

263.3 

V 

1573.6 

5.3 

-0.139 

2.723 

1373.0 

1776.4 

7 

d 

1557.9 

5.0 

0.163 

2.535 

1398.3 

1768.7 

r 

w 

1841 .1 

3.2 

0.596 

3.351 

1719.8 

2039.6 

p 

18.31 

17.7 

0.268 

1.975 

13.58 

26.76 


H.57 

22.5 

0,795 

3.629 

8.05 

25.58 

Table 2.4 

. : Statistical 

Analysis 

of Fixed 

lot BK4 



n = 

100 





Parameter 

Ifean 

8 (? ) 


Ok 


Range 






Lower 

Upper 

L 

235.8 

0.3 

0.424 

2.637 

234.0 

238.0 

■ B 

116.3 

0. 6 

-0.030 

3.556 

114.2 

117.9 

H 

73.6 

1.6 

0.097 

2.150 

71.4 

75.9 

A 

274.1 

0.7 

0.006 

2.941 

268.4 

278.5 

¥ 

2016.6 

1 .8 

0.245 

2.622 

1945.7 

2101.9 

7 

1611.1 

1.7 

-O.6I9 

3.8^ 

1509.8 

1668.2 

r 

1899.6 

1.6 

-0.522 

4.671 

1779.2 

1976.6 

w 

p 

17.92 

6.1 

.0.179 

6.534 

14.13 

22.22 

^b 

7.60 

16.9 

-0.437 

3.301 

3.84 

10.99 
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Table 2.5 ; Statistical Analysis of fixed lot B£5 

n = 97 


Parameter 

Mean 

6 (% ) 


Gk 

Range 

Lower Upper 

1 

231.4 

0.7 

0.888 

2.728 

230.0 

235.0 

B 

119.2 

1.4 

- 0.998 

7.915 

110.7 

123,5 

H 

79.7 

1.8 

-2.513 

14.926 

71.2 

82.4 

A 

275.8 

1.9 

0.270 

3.238 

260.1 

290.2 

Y 

2197.3 

2.5 

0.352 

4.482 

2012.3 

2368.2 

y 

d 

1174.0 

2.2 

1.976 

10.815 

1700.1 

1982.8 

7 

T*T 

2036.0 

2.1 

1.913 

11.279 

1927.4 

2261.1 

W 

p 

14.78 

5.1 

- 0.399 

4.546 

12.05 

17.05 


14.16 

16.3 

0.187 

2.097 

9.52 

19.53 


Table 2 , 

6 : Statistical Analysis of fixed 
n = 100 

lot BK6 


Parameter 

Mean 

6 (.%) 

Os 

Ok 

Range 


1 




Lower 

Upper 

1 

. 231.0 

0.5 

0.645 

3.023 

229.0 

234.0 

B 

117.8 

0.7 

0.576 

3.331 

116.0 

120.5 

H 

78.8 

1.2 

0.307 

2.665 

77.0 

81 . 3 

A 

272.2 

1.1 

0.651 

3.462 

266.3 

231.0 

¥ 

2146.0 

1.8 

0 

CD 

0 

2.996 

2064.0 

2251.0 

■^d 

1795.1 

1.3 

- 0.253 

2,859 

1724.4 

I 846.6 

r 

w 

2038.6 

1.1 

-0.126 

2.419 

1976,6 

2082.1 

p 

13.57 

4.7 

-1.188 

6.123 

10.71 

14.63 

^b 

10.85 

15.1 

0.202 

2.421 

7.50 

14.39 
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Sable 2.7 : Statistical Analysis of Fixed lot BK7 

n = 100 


Parameter 

Mean 

6 (% ) 

Gs 


Range 






lo; 7 er 

Upper 

L 

252.7 

0.9 

- 0.081 

3.070 

227.0 

238.0 

B 

113.2 

1.4 

-0.372 

2.946 

108.9 

117.0 

H 

70.5 

2.1 

0.277 

3.224 

67.5 

74.8 

A 

263.4 

2.0 

-0.224 

3.101 

249.7 

275.0 

¥ 

1857.9 

2.9 

0.060 

3.022 

1722.9 

2003.5 

y 

d 

1606.0 

2.0 

0.334 

5.520 

I5I8.4 

1693.3 

W 

1 885 . 1 

1.8 

-0.115 

3.077 

1788.3 

1973.4 

p 

17.39 

5.1 

-2.411 

11.928 

12.68 

18.94 


8.13 

18.6 

0.754 

3.130 

5.54 

12.89 

Table 2.8 

; Statistical Analysis of 

Fixed 

lot BK8 ■ 



n = 

100 





Parameter 

ifean 

6 (% ) 


Gk 

Range 






lower 

Upper 

1 

232.8 

0.7 

0.170 

2.031 

230.0 

236.0 

B 

113.8 

1.3 

0.681 

3.141 

111.0 

119.0 

H 

76.1 

1.6 

0.168 

2.193 

73.5 

78.5 

A 

265.0 

1 .8 

0.602 

2.482 

257.6 

279.7 

¥ 

2015.8 

2*4 

0.566 

3.075 

1931.4 

2153.5 

'"d 

1607.6 

1 .6 

-0.174 

2.987 

1532.5 

1659.9 


• 

0 

CT^ 

1 .5 

-0.243 

4.196 

1797.2 

1981.3 

P 

18.47 

. 6.3 

5.420 

21.262 

16.12 

25.16 

^b 

8.58 

17.6 

0.324 

2.582 

5.54 

12.40 
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Table 2,9 • Statistical Analysis of !fi.xed lot 1(BK9) 

n = 295 


Parameter 

Ifean 

6(f. ) 



Range 

lov/er Upper 

1. 

228.5 

1.6 

-0.423 

3.121 

215.5 

236.5 

B 

109.4 

2.3 

-0.265 

3.217 

101.8 

116.0 

H 

63.3 

3.3 

-0.045 

2.483 

57.3 

68.8 

A 

250,1 

3.3 

-0. 204 

3.504 

220.4 

273.2 

Y 

1584.3 

5.5 

0.062 

2.446 

1357.6 

1784.5 

7 

d 

1 61 6, 4 

6.2 

Om 666 

3.705 

1398.3 

2023.4 

r 

w 

P 

1858.1 

4.8 

0.768 

3.255 

1691 .8 

2202.4 

15.09 

22.3 

0.820 

3.583 

7,53 

26.76 

^b 

19,89 

31.0 

0.589 

3.005 

8.05 

39-69 

Table 2, 

.10 : Statistical Analysis of Mixed lot 2(BK10) 

n = 497 

Parameter 

■ Sfean 

6(fc ) 

Gs 


lower 

Range 

Upper 

L 

232.7 

1 .0 

0.254 

2.001 

227.0 

238,0 

B 

116.0 

2.5 

-0.092 

2.527 

108.9 

123.5 

H 

75.7 

4,8 

-0.315 

1.991 

67.5 

82.4 

A 

270.1 

2.4 

-0.307 

3.010 

249.7 

290.2 

Y 

2045.8 

6.2 

-0,227 

2.389 

1722.9 

2368.2 

^d 

1678.2 

5.5 

0.430 

1.765 

1509.8 

1982.8 

r 

1952.3 

3.9 

0.395 

2.328 1779.2 

2261.1 

w 

P 

16.43 

12.9 

-0.049 

2,776 

10.71 

25.16 

^b 

9.84 

29.8 

0.885 

3.323 

3.84 

19.53 
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lable 2,11 : Statistical Analysis of Random lot(BEll) 

n = 1 1 3 


Parameter 

Ifean 

6(%) 



Range 

Lower Upper 

L 

2??. 2 

2.1 

-1 .012 

9.723 

201 .0 

241 .0 

B 

109.2 

2.8 

-0.052 

2,940 

102.0 

116.0 

W 

2815.5 

4.8 

0.630 

2,925 

2542.0 

3220,0 

P 

12.05 

22.8 

0.291 

3.299 

6.00 

19.40 


22.18 

29.8 

0.451 

^.488 

00 

• 

38.85 


I 


Frequency 



L (mm) 

(a) Histogram 



LCmm) 

(b) Cumulative distr ibutipn 


Fig.2^1 Histogram of length of brick of 
fixed lot 8K1 
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Relative frequency 




2%0 


I 


n = 100 -0-25 

Lm=232*8mfn 

6 =0*7 % 

- 0*20 



231-7 233-4 J 

L(mm) 

(a) Histogram 


235-1 236-8 


-normal 
log normal 


JQ 

5 0-4 

0 

01 



230-0 231-7 233-4 235-1 236*8 

LCmm) 

(b) Cumulative distribution 

Fig.2-2 Histogram of length of brick 

of fivoH lot 


Relative frequency 




Probability Frequency 



B (mm) 

(q) Histogram 
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Fig.2'3 Histogram of breadth of brick of 
fixed lot 3K1 


Relative frequency 




Probability Frequency 



Relative frequency 





Probability Fr«qu«ncy 
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H(mm) 

(a) Histogram 



Fig. 2-5 Histogram of height of brick of 
fixed lot BK 1 


Relative frequency 




Probability Fr«qu«ncy 



n=100 

HrT5=^76»1mm 



1 

|6 =1-6% 




lOh 



» 

0*8 h 


7S*U 

H (mm) 

Cc) Histogram 


0*05 


normal 
log normal 


0*0l L_ 

74*0 75-14 76 

H (mm) 

(b) Cumtilativo distribution 

Fig 2-6 Histogram of height of brick 
of fi Ksd lot BK 8 


Rolativo frequency 





Probability Frequency 



A (cm^) 


(a) Histogram 



Ftg.2*7 Histogram of area of brick of 
fixed lot BK \ 


Relative frequency 




Probability Frequency 
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n=100 

A = 265-0 cn? 
8 = 1 - 8 % 


0-4 

0-3 


20 


10 


01 1 1 I 1 i i ■ -V— I 

258-0 264-0 270-0 276-0 2820 

ACcrr?) 



0-1 


0-0 


(a) Histogram 



A(cm^) 

(b) Cumulative distribution 

Fig. 2-8 Histogram of area of brick of 
fixed lot BK 8 


Retotive frequency 










Vm=20t5-8cfT? 
6 =2*4% 



1931-0 1994-4 2057-8 2021-2 2084-6 

V(cm^) 

(a) Histogram 



•normal 

lognormal 


n 

jo 0-4 
o 

t. 

a. 



19310 1994-4 2057-8 2021-2 2084-6 

V(cm^) 

(b) Cumulative distribution 

Fig.2-10 Histogram of volume of brick 

, 1 ^# #1 ii # 01^0 


Relative frequency 





ProbabiUty Frequency 
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>d(kg/n^) 


(a) Histogram 



14850 1559*86 1634*72 1709-58 1784*44 
7ij<kg/f!?) 

(b) Cumulative distribution 

Fig.2-11 Histogram of dry density of brick 
of fixed tot BK1 


Relotive frequency 



ProbabUHy Frequency 




(b) GumulatWe distribution 


Fig.242 Histogram of dry density of brick 
of fixed lot BK8 


4S 


Relative frequency 





Probability Frequency 



1692-0 1770-86 1849-72 1928-58 2007-44 

r^Ckg/m^) 

(b) Cumulative distribution 

Fig. 2T3 Histogram of wet density of brick 


Relative frequency 




Probability Frequency 


40 


•4 


50 



(kg/m^) 


(a) Histogram 



(kg/rr?) 


(b) Cumulative distribution 

Fig 2*14 Histogram Of wet density of brick 
of fixed lot BK8 


Relative frequency 





Probability Frequency 



80 10-85 13'7 16*55 19A 

P (%) 

(a) Histogram 



Fig. 2-15 Histogram of water absorption of 
brick of fixed lot BK1 



PfobabllHy Fr«qu«ncy 


n =100 
D-= 18-47% 


(a) Histogram 


normal 
log normal 


/ 




^ 0 

1 

1 ! 

1 1 


(b) Cumulative distribution 

Fig 2 t6 Histogram of water absorption 
of brick of fixed lot BK 8 


Relative frequency 







Probability Froquency 


fj^^=23'69 N/m 
S-23*2% 


26’28 35*A2 


(a) Histogram 


normal 
log normal 


(N/mm*") 

(b) Cumulative distribution 

Fig. 2-17 Histogram of compressive strength 

fixed lot BK1 


Relative frequency 









Probability Frequency 


(a) Histogram 


normal 
log normal 


( H/mm > 

(b) Gumulati ve distribution 

Fig.218 Histogram of compressivo 
^ of fixed lot BK8 


Relative frequency 









(N/m‘rn?) 





Probability Frequency 


n = m 

f|^^=22-18N/mm* 
S =29*8% 


14*13 41*84 


(a) Histogram 


normal 
log normal 



34*13 41*84 

) 


(b) Cumulotive distribution 

Fig 2*20 Histogram Of compressive strength 
of brick of random lot BKIt 


Retotive frequency 
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/ 


2,12 j Summary of Briclc Properties 









Parameter 

Mean 

6(^ ) 

Range 





Lower 

Upper 




1 , ,1 II , i 



typical 

Is 

230.7 

1.2 

224.7 

236.5 

Pixed i-ot 

BK1 : 

H 

64.5 

2.4 

60.2 

68-0 


P 

d3.41 

14.7 

8.2- 

18.43 


fb 

23.69 

■ 23.2 

13.92 

39.35- 

gc . 






Typical 

Is 

232.8 

0,7 

230.0 

236.0 

Pixed lot 

BK3 

H 

76.1 

1.6 

73.5 

78.5 


P 

18.47 

6.5 

16.12 

25.16 



8.58 

17.6 

5-54. 

12.40 

Mixed lot 1 

L 

228.5 

1.6 

215.5 

236.5 

(m9) 

H 

63.3 

3.3 

57.3 

68,§ 


P 

15.09 

22.3 

li53 

■26,76 


f-u 

19*89. 

31 . 0 . . . 

,8.0§ 

■ 39,69 


-^b • 



,,,U . ... . ^ XT 7 ^ 


rm^- i: 


. -a* ; -t. 

Mixed lot 2 

1 

232.7 

1.0 

227.0 

238.0 

(WO) 

H 

75.7 

4.8 

67.5 

82.4 

- 

P 

16.43 

12.9 

10.71 

25*16 


• ^b ■ 

, .9.84 

. . 29.8 . . 

. . 3.84 

19.53 

aandoffl lot 

Ii 

227.2 

2*1 

2 « 1 .| 

241.0 

. (BK11) 

H , 


-i. 




P 

12.05 

22 ^ , 

6.00 

I9.4i 


fb 

22.18 

29.8 

11.38 

38.85 
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available in the literature to fit a mathematical model to 
the observed data. Chi-square goodness of fit test is 
applied to find the suitability of an assumed mathematical 
model in the present investigation. The value of the 


ohi-square is computed from the formula (94) 

2 _ ; (0,-E,)2 


( 2 . 10 ) 


where 

y2 

^ — the value of ohi— square 

Oj_ = the observed frequency in the ith interval 
®i “ expected frequency corresponding to the 

assumed mathematical model in the ith interval 
k = the number of categories or class intervals 
considered in the test. 

The number of degrees of freedom is given by 


V = k - r - 1 (2.1l) 

where 

^ = number of degrees of freedom 
and r = number of parameters estimated from the data. 
OJhe chi-square value computed from Eq. 2.10 for an assumed 
distribution is checked with the chi-square values corres- 
ponding to different significance levels at which the assumed 
distribution will be accepted. The maximum significance 
level at which the test is satisfied j called p- level (93) 
is also computed. Hormal and lognormal distributions are 
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considered in tlie present investigation to fit the observed 
data. 


formal Distribution 

The probability density function of a Formal 
Distribution is given by (93) 

1 4 X— Xm 2 

fy(x) — esp r ~ i - °° < X < “> (2,12) 

sfaiT 2 s M “ “ 

where 

= the mean value of random variable X 
s = standard deviation. 


The probability distribution function is given by 


Mx) = / 


1 


8 ^; 
x-Sm 


0 ( 


s 


exp 

) 




] dx 


(2.15) 

(2.H) 


where 


0 (,) = the cuaiulative distribution function of a 

standardized normal random variable. 

The parameters of the normal distribution are estimated from 
the data throu^ Eq. 2,2 and Eq, 2.5 respectively. 


Iiogporm-al Distribution 

The probability density function of a ."lognormal 
distribution is given by (93) 
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fx(x) 


1 1 In X - In m .2 

— exp ) ] 

x V an tin ^ 

0 < 2 < “ (2.15) 


where m and are the two parameters of the lognormal 
distribution, The probability distribution function is 
given by 




/ 7= 

o xY^ 

ln(x/m) 


0 ( 


^In 


exp [ - 


) 


1 In X - In m 

sir- 


2 

) ] dx 


( 2 . 16 ) 

(2.17) 


The parameters of the lognomal distribution aie expressed as 

aj = In (6^ + 1) (2.18) 

In 

m = (*■ ) (2.19) 


2.2,6 Statistical results 

Chi-square test is carried out to fit normal distribution 
and lognormal distribution to the data of different parameters 
of the samples. The parameters of normal distribution are 
estimated from the data by Eas. 2.2 and 2.3. Ihe parameters 
of lognormal distribution m and estimated from 

the data by Eqs. 2.18 and 2.19. The maximum significance 
level ( p-level ) at i^feich the test is satisfied is computed 
for the two distributions and is given in Table 2.13. 

Computed p— level less than 0.5 percent is not listed m the 
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table. Usually a null hypothesis that the data of a 
particular random variable follow an assumed distribution 
is accepted at a level of significance 1 %, 5 %and 10 %, 

The p— level is computed to see whether a parameter like 
length, breadth etc. follows normal or lognormal distribution 
at a particular significance level in all the cases. Each 
parameter is discussed below seperately. 

Length 

Mean length of bricks of different fixed lot varied 
in the range of 227 to 236 mm and the coefficient of variation 
varied from 0.3 to 2 % . The lower and upper range is 
observed to be within + 5 mm about the mean value of a 
particular lot. Samples of BK4, BK5, BKS, BK7 and BK8 could 
not be fitted to normal or lognormal distribution. The 
data is clubbed around the mean values which signifies that 
there is a deterministic trend rather than probabilistic, 
Lo^aormal seems to be a better fit over normal distribution 
at 1 % significance level for the data of random lot while 
normal distribution appears better for fixed lots BK1 , BK2, 

BK3 and mixed lot BK9. Length of brick can be taken as 
normally distributed for practical purposes. The local 
standard of Kanpur zone is 230 mm. The probability of 
getting length of brick less than 230 mm can be obtained 
as (from mixed lot 1), 
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, , 250-228*5 

P (I. < 230) = 0 ( ^ ) = 0 (0.41) 

3*o56 

= 0.658 

wh-icli is same as 658 in 1000 bricks (asstunlng noisaai 

distnibution) . 'The probability of getting length of brick 

less than 230 mm (from mixed lot 2), assuming normal, is 

, ^ ^ ^ 230-232.7 

P (I> < 230) = 0 ( ) = 0(-l.l6) = 0.123 

2.327 

which is same as 123 in. 1000 bricks. 

Breadth 

Mean breadth of bricks of different fixed lots varied 

from 108 mm to 119 mm and the coefficient of variation varied 

from 0.5 to 2.8 percent. Normal distribution fitted well at 

5 percent level of significance except in two cases. Normal 

also fitted the random lot at 5 percent level, llthough 

lognormal appears better in case of random lot, normal 

distribution seems to be a better fit in most of the cases. 

Thus breadth of brick can be taken as normally distributed 

random variable. Assuming the IocelI standard of breadth of 

bricks 1 lOMU (for mixed lot 1 ) 

110-109.4, 

P (B < 110 ) = 0 ( -) = jZJ(0.24) = 0.595 

2.5162 

which is same as 595 in 1000 bricks. 

Por mixed lot 2, the probability of getting breadth of bricks 
less than 110 mm is 
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P (B < 110) = 0 ( ) = 0 (-2.25) = 0.012 

whicli is same as 12 in 1000 bricks. 

Height 


Mean height of bricks of different manufacturers of 
Kanpur zone (BK1, BK2 and BK3 lot) varied from 62 mm to 
64»5 mm and the coefficient of variation vairied from 2*4 to 
3.4 percent. The coefficient of variation of height of 
bricks appears to be higher as compared to length and breadth 
of bricks. Mean height of bricks of different manufacturers 
of Paroke zone varied from 70,5 mm to 79,7 mm and the 
coefficient of variation varied from 1.2 to 2.1 percent. 
Normal fits the data at 1 percent level except in few cases. 


The local standard of height of bricks is 63 mm in 
Kanpur zone, and the probability of getting less than 65 mm 
(from mixed lot I) is 


P(H < 63) = 0 


63-63,3 

63. 3x0.033 


= 0 (-0.144) 


= 0.444 


which is same as 444 in 1000 bricks. 


The local standard of height of bricks in Faroke zone is 
75 mm, the probability of getting less than 75 mm (from 
mixed lot 2) is 


P(H < 75 ) = 0 ( 


75-75.7 

75.7x0.048 


) = 0 (-0.193) 

= 0.424 


which is seone as 424 iii 1000 bricks. 
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Area 

The cross-sectional area of bricks is observed to 
have higher variability as compared to the individual 
variability of length and breadth. The coefficients of 
variation of area of bricks are found to be 3*3 percent 
and 2,4 percent for mixed lot 1 and 2 respectively. In 
both the cases, a marginal negative skewness is observed, 
Uoimal and . lognormal distribution fits the data of all 
the fixed lots (except BK6 and BKS) at 1 percent significance 
level. The local standard of bricks of Kanpur zone, is 

1 = 230 mm, B = 110 mm and H = 63 mm 
A = 230x110 = 25300 mm^ = 253 cm^ 

2 

Probability of getting area of bricks less than 253 cm is 

253-250.1 

P(A < 253 ) = 0 ( ) = 0 (-0,35) 

250.1 X 0.003 

= 0,363 

which is same as 363 in 1000 bricks (assuming nozmal). 
Similarly for Faroke (Kerala) zone, the local standard is 

1 = 230 mm , B = 110 mm and H = 75 mn 

A = 25300 mm^ = 253 cm^ 

253-270.1 

Thus, p(A < 253) = 0 ( ) - 0 (-2.64) 

270.1 X 0.024 

= 0,004 

which is same as 4 in 1000 bricks (assuming nomal). 
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Yolmae 

The coefficient of variation of volume of bricks 
are found to be 5.5 percent and 6,2 percent for mixed lot 1 
and 2 respectively, Formal distribution fits the data of 
all the fixed lot at 1 percent significance level whdreas 
lognormal fits at 5 percent. lognormal distribution 
is found to be a better fit at 1 percent significance level 
in some of the cases. 


The volume of a brick as per local standard of 
Kanpur zone is 

V = 250 X 110 X 63 = 1593900 mm^ = 1593.9 cm^ 

Probability of getting volume of bricks less than 1593.9 cm , 

assuming normal distribution, is given by (for mixed lot 1) 

^ 1593.9 - 1584.5, 

P (Y < 1593.9) = ^ (:: ) = ^ (0.11) 

— '\r 


1584.3 X 0.055 


= 0.543 


which is same as 543 in 1000 bricks. 


Similarly, volume of a brick as per local standard! of 
Paroke (Kerala) zone is 

¥ = 230 X 110 X 75 = 1897500mi?^= 1897.5 cm^ 

1897.5 - 2045.8 


and 


P (¥ < 1897.5) = 0 ( 


2045.8 X 0.062 
which is same as 121 in 1000 bricks. 


) = 0 (-1.17) 

= 0.121 
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Dry Density 

Ilie CO efficients of variation of dry density of 
bricks are found to be 6,2 percent and 5.5 percent for 
mixed lot 1 and 2 respectively, SToimal distribution fits 
the data of all the fixed lots (except BK5) at 10 percent 
significance level. Lognormal distribution also fits at 
10 percent significance level except BK5 lot* In some of 
the lots lognomal appears a better fit over nomal 
distribution and vice versa, 

?/et Density 

The coeff icienis of variation of wet density of 
bricks are found to be 4.8 percent and 5.9 percent for 
mixed lot 1 and 2 respectively which is slightly lower 
as that of dry density. Mixed lots did not fit to any 
of the distribution whereas fixed lots, except BK5,fit 
normal and lognormal distributions at 1 percent significance 
level. 

Percent of Water Absorption 

The mean value of different fixed lots varied 
between 13.41 to 18.47 with coefficient of variation varying 
from 4.7 percent to 17.7 percent. Percent of water absorption 
of mixed lot 1 has a mean value 15.09 with a coefficient of 
variation 22.3 percent whereas mixed lot 2 has mean 16,43 
and coefficient of variation 12.9 percent. Mean of percent 
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water absorption is found to be 12,05 (ranging from 6,0 to 19.4) 
■VTitb. a C,0,Y, = 22,8 percent for the random lot, Beth 
normal and lognormal did not fit even at 0.5 percent 
significance level except in three cases. 

Brick Strength 

Mean compressive strength of brick of fixed lots BK1 , 
BK2 and BK5 is in the range of 14.57 M/mm^ to 23,69 
with C.O.T, varying from 22,5 percent to 24.6 percent. Mean 
brick strength of fixed lots BK4, BK;5, BE? and BK8 is in 
the range of 7.6 B/mm^ to 14.16 E/mm^ v/ith C.O,?, varying 
from 15.1 percent to 18,6 percent. Mean value of mixed 
lot 1 and mixed lot 2 are 19.89 E/mm and 9.84 E/mm which 
clearly represents that the quality of bricks of Kanpur zone 
are better than that of Baroke (Kerala) zone. The 
coefficient of variation of mixed lot 1 , mixed lot 2 and 
random lot are found to be around 30 percent which signifies 
variability in brick manufacturing from manufacturer to 
manufacturer, Eoimal distribution fits the data of all 
the fixed lots at 2.5 percent significance level while log- 
noimal could not be fitted to the data of fixed lot BK4 
even at 0.5 percent significance level, Eoxmal did not fit 
mixed lots and random lot at 1 percent signiiicance level. 

Thus, normal distribution can be accepted at 2.5 percent 
significance level to represent the probability distribution 
of compressive strength of brick. 



2«5 Statistical Inalysis of IThickness of Mortar Joint 

Ihe thickness of mortar joint has two aspects in 
masonry construction, (i) aesthetic aspect and (ii) strength 
aspect. Too much variability of joint thickness in a wall 
does not give a good appearence to the observer or ovwier 
of the building. The thickness of mortar joints play a 
great role on the strength of brickwork. Tests carried out 
by Francis et al. (25) showed that prism compressive strength 
is a function of average joint thickness and reduces as the 
joint thickness increases. The strength of prism decreases 
at a faster rate with increase in joint thickness for 
perforated bricks as compared with the solid bricks. So, it 
is expected that the variability of joint thickness will 
affect the strength of brickwork although no such theoretical 
relationship is established. To study the variability of 
joint thickness, data is collected from different existing 
structures and statistical analysis is carried out, and is 
presented in the following sections. 

2.3.1 Collection of field data and notations 

Three different buildings of exposed brick joint 
wall were selected for study. Horizontal joint (called bed 
joint) thickness and vertical joint thickness were measured 
at different locations of the wall selected at random, 100 
measurements were carried out for each bed joint and vertical 
joint thickness of each building. Thickness was measured by a 



slide calipers from exposed surface of tlie wall upto an 
accuracy of 0,02 mm, Tlie three sets of data for bed joint 
thickness and vertical joint thickness were combined for 
the statistical analysis of the combined sets, fhe following 
notations are used in the subsequent discussions. 


Building Ho. 

notation 

Remarks 

Building 1 

MTH1 

Bed joint 


MTV1 

Vertical joint 

Building 2 

MTH2 

Bed joint 


MTV2 

Vertical joint 


MTH3 

Bed joint 

Building 3 



■MTV3 

Vertical joint 

Bed joint data of three 
buildings combined 

MTHC 

Bed joint 

Vertical joint data of 
three buildings combined 

MTVG 

Vertical joint 


One hundred measurements of bed joint thickness of mortar 
joint were taken after breaking some brick columns. These 
are the actual joint thicknesses inside the columns and 
are refered in the text as MTH4. 

2,3,2 inalysis of data 

The data of bed joint thickness and vertical joint 
thickness are analysed and the result is given m Table 2, 


Table 2,14: Statistical Analysis of Thickness of Mortar Joint 


designation n 

Mean 

) 

^s 


: 

lower Upper , 

MTH1 

100 

13.48 

12.4 

0.552 

2.879 

10.50 

17.60 

MTH2 

100 

13.02 

12.8 

0.261 

2.586 

9.82 

17.00 

MTH3 

100 

12.28 

14.0 

0.124 

2,986 

8.80 

16.50 

MTVl • 

100 

12.21 

14.2 

0.770 

3.255 

9.50 

17.00 

MTV2 . 

100 

11.78 

19.5 

0.421 

2.499 

7.50 

17.50 1 

MT73 • 

100 

12.10 

15.1 

0.075 

2.434 

8.50 

16.00 

MTHG . 

300 

12.93 

13.5 

0.245 

3.005 

8.80 

17.60 

MTIG 

300 

12.03 

16.4 

0.540 

2.767 

7.50 

17.50 

MTH4 

100 

14.90 

11.5 

-0.086 

2.736 

10.60 

19.00 
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The coefficient of variation of bed ^oint thicloiess varied 
between 12 to I 4 percent in all the buildings whereas that 
of vertical joint thickness varied from 13 to 20 percent. 

The lower and upper range has a deviation of + 5 eeq from 
the mean in all the cases. Chi— square test for normal and 
lognormal distribution is conducted and the result of 
chi-square test is given in Table 2.15. ill the data 
satisfies log-normal distribution at a level of significance 
of 2,5 percent whereas all of them satisfies normal distri- 
bution at a level of significaoje of 1 percent. Histogram 
and the fitted distribution of two typical bed joint thickness 
(MTHI) and (MTH 4 ) are given in Pig. 2.21(a) and Pig. 2.21(b) 
respectively, 

2,4 Statistical Analysis of Strength of Mortar 

Mortar is used as the binding material between brick 
units. The mix proportion of cement sand mortar in load 
bearing brick work and reinforced brickv/ork varies from 1 :3 
(cement :sand) to 1:6, The water cement ratio varies in the 
range of 0.7 to 1 .0 and it is usually governed by the require- 
ment of mortar. Although there is considerable emphasis on 
the quality maintenance of bricks, there is hardly any 
qualitative test carried out in the field on the strength 
of mortar. The cement sand mortar needed for two to four 
hours of working period is prepared at one time and the 
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Table 2,15: Chi-square Test Result of Thickness of Sfortar Joint 


5>esig- ic-mnal Ifonaal 

nation 




s 

X2 

if 

P -level 
(% ) 

m 

^n 

X2 

i 

p-leve 
( fo ) 

MTH1 

13.48 

1.67 

16.80 

6 

1.0 

00 

m 

0.123 

10.55 

6 

10.0 

MTH2:: 

13.02 

1.67 

4.35 

6 

62.0 

12,91 

0.128 

4.18 

6 

65.0 

MTH3~ 

12.28 

1.71 

10.84 

6 

7.5 

12.16 

0.159 

14.01 

6 

2.5 

MTV1 

12.21 

1.73 

16.54 

6 

1.0 

12.09 

0.141 

9.72 

6 

13.0 

MTV2_. 

11.78 

2.29 

12.83 

6 

2.5 . 

11.56 

0.193 

8.53 

6 

20.0; 

MTV3‘ 

12.10 

1.83 

7.87 

6 

24.0 

11.96 0.1‘50 

13.40 

6 

2.5 

MTHO ■ 

12.93 

1.75 

15.86 

8 

2.5 

12,81 

0.135 

15.52 

8 

2 . 5 ; 

MTYG 

12.03 

1.97 

8.62 

8 

37.0 

11.87 

0.163 

5.88 

8 

66,0; 

MTH4 

1 4.90 

1.71 

2.21 

6 

89.0 

If. 80 

0.115 

4.90 

6 

55.0 


Frequency . Frequency 
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Fig 2'21 Histogram of thickness of mortar joint 


Relative frequency Relative frequency 
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consistency of mortar is maintained by constant addition 
of water and remixing the material. This is a common practice 
observed in India even in large brickwork construction. The 
strength of masonry depends on the quality of bricks and 
mortar used in the construction, A study of the variability 
of mortar strength is needed in order to understand the 
behaviour of brickwork strength. lue to the absence of 
data from field, mortar cubes were cast, cured and tested in 
the laboratory after 28 days. The statistical analysis of 
mortar strength data is presented in subsequent sections, 

2,4-.1 Sample collection 

Mortar cubes of 70,7 mm size were cast for three 
different mixes i.e,, 1:3 , 1 54 and 1:5 (cement :sand by 
weight) and tested after 28 days of water curing in the 
laboratory. Six cubes were cast at a time for each mix 

is called a set. A fresh mix has been done for casting 
every six cubes. The amount of water is kept constant at 
each set of a particular mix. 

2,4.2 Histogram and statistical analysis 

The variation of mean value and coefficient of 
variation of compressive strength of mortar fj^CH/mm ) of 
different sets of each mix is shown in Hig, 2,22 through 
Fig, 2.24, 411 the sets of a particular mix is combined and 

statistical analysis is done. The statistical analysis of 



(N/mm?) 8 (•/«) 
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Set number 

(b) Variation ot mean value 


Fig. 2-22 Set variation of strength of mortar 
{1:3 mix) 



(%) s 


10 


0 



Set number 

<a) Variation of coefficient of variation 



1 5 9 13 17 21 25 

Set number 

(a) Variation of mean value 


Fig.2-23 Set variation of strength of mortar 
(1:4 mix) 




( N / m m 


2 


SaX number 

(a) Variation of coefficient of variation 


xr-o-a>^ P'-o-o-^ 

0 - « o^O'o-o-^ 

0 - 0 ^ ^ 0-0 


y-^ 9 13 17 21 

Set number 


(b) Variation of mean value 


Rg.2'24 Set variation of strength of mortar 
(1:5 mix) 
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the three mixes are given in Table 2.16. Chi-sq.uare tests 
for normal and lognoimal distribution is carried out for 
each mix. The results of chi-sq.uare test aud the mgylTmim 
significance level (p-level) at which the test is satisfied 
is given in Table 2,17. The histograms aud the fitted 
distributions of three mixes are shown in Pig. 2.25 throu^ 
Pig. 2.27. 

The coefficient of variation of the strength of 
mortar cubes varied from 10 percent to 18 percent. Data 
of all the three mixes are found to be slightly positively 
skewed. The variation of the mortar strength in the field 
cubes is likely to be more by twenty percent and appears to 
be consistent with that of the concrete cubes. Normal 
distribution fits the data at 10 percent level of significance 
in all the cases, whereas lognormal distribution fits at 
5 percent level of significance. Normal distribution appears 
to be a better fit over lognormal distribution. 


2,4.5 -Analysis of probability of failure 


The characteristic strength of mortar (fjjjjj-) of a 
particular mix for an accepted probability of failure (p^) 
can be computed as 


Vf ■” ^ — ^mk^ 


( 2 . 20 ) 


where 


and 


^mk 

Pf = 


characteristic strength of mortar 
accepted risk or probability of failure. 
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Table 2.16: Statistical Analysis of Strength of Mortar 


Mix 

Designation 

n 

Mean 

:6(? ) 



Range 

lower Upper 

1:3 

168 

21.54 

13.3 

0.290 

3.207 

15.20 

31.64 

1 :4 

168 

14.29 

10,0 

0.073 

2.969 

10.05 

17.89 

1 :5 

156 

10.42 

18.1 

0.128 

2.252 

6.35 

15.00 


Table 2.17 J Ghi -Square Test Result of Strength of Mortar 


Mix Sbrmal log Hbrm^ 

Design— . ■ ' ■. ■ ■ ■ ■” 

nation s tC ^ p-level m ^ 

( %. ) ^ ^ ) 

1:3 21.54 2.86 4u78 6 57.0 21.35 0.132 4.14 7 76.0 

1:4 14.291. 443.32 7 85.0 14.22 0. 1005.86 7 55.0 

1:5 10,42 1.89 9.90 7 19.0 10.26 0.179 12.44 7 7.5 
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Fig.2*25 Histogram of strength of mortar {|:3mix ) 


Frequency 
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Fig.2*26 Histogram of strength of mortar (1:4 mix) 



Frequency 
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Fig 2*27 Histogram of strength of mortar (V«5mix) 



Assuming strength of mortar follows normal distribution, 
Eq., 2*20 can be written as 


84 - 


Pf = 


0 ( 


^mk *■ 


■mm 


s 


( 2 . 21 ) 


or 


^mk ^mm 


( 2 . 22 ) 


where fjjjj^ = mean strength of mortar 
s = standard deviation 

k = 0 (p^) 

and 0 ( . )= inverse cumulative distribution function of 

a standardized normal random variable. 

If strength of mortar follows lognormal distribution, 

Eq.* 2.20 can be written as 


p = ) 

‘^In 


‘•mm 


or 


‘■mk 


Y(6 +1) 


exp [ fc ir(ln(6^+l))] 


(2.23) 

(2.24) 


where ^ 

^ “ f 

imn 

= ln(6^+1 ) 

2 

and m = fmm exp ( - ) 

For an accepted probability of failure of 5 percent 
(i.e, , Pf ~ 0.05) 

k = 0 "^(Pf) = - 1-645 
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and Eq.s» 2,22 and 2,24 becomes 

^mk ^mm ~ ‘'•^4-5 s (2,25) 

^ " fT ^+V) ''*6^5f(ln(6^+l))] (2,26) 

Eqs, 2,25 and 2,26 represents the relationship between 
characteristic strength and mean strength for noimal and 
lognormal distribution respectiTely, 

2,5 Statistical Analysis of Masonry Strength 

2.5.1 General 

Brick masonry is defined as an assembly of brick 
units bonded together by mortar in a predetermined 
orientation. Hie strength of brick masonry depends on 
several parameters. Some of the parameters are given below: 

i) the strength of individual brick xuiits 

ii) the strength of mortar 

iii) the thickness of mortar joint 

iv) layout and orientation of bricks. 

Since strength of brick and strength of mortar are 
themselves random variables , the strength of masonry also 
becomes a random variable. Ihe variability of masopry 
strength depends upon several factors j some of which are 
listed below: 

i) variability of strength of bricks 
ii) variability of strength of mortar 
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vaxiability of joint thickness 
iv) variability of workmanship. 

The effect of workma i ^hip and other variables on 
strength of masonry has been discussed by Hendry (56) in a 
great detail. To study the statistical behaviour of masonry 
strength, sufficient data is needed which presently is not 
available in India. One of the main reason is that 
regular testing of prisms is not carried out to ensure the 
quality of brickwork and workmanship even in large 
construction projects. Deterministic formula derived from 
experimental investigations (1 , 2) and SGPI recommendation 
(13,95) are used to simulate masonry strength to study the 
statistical behaviour. 

2.5.2 Masonry strength 

Prom the experimental investigations undertaken at 
I.I.T. Kanpur, the strength of brick masonry can be obtained 


as (1 , 

2). 




'w = 


(2.27) 

where 

f = 

J-W 

strength of masonry 



^b = 

strength of brick 



^m = 

strength of mortar 


and 


a coefficient which depends 

on the layout of 


brick and joints. 



Value of is found to be 0.275 and 0.42 for loading 
perpendicular and parallel to bed joint respectively. The 
dean, standard deviation and coefficient of variation can 
be estimated (96) from Eq^. 2.27 as (see Apendix A) 


^wm 

= k f (f^ f ) 
w ^ ^ Dm 

(2.28) 

®w 

^ r- , f ^bm 2 ^mm 2 \l/2 

= 0-5 v 7- ■ < + -^b ) 

^mm ibm 

(2.29) 

and 6 

w 

= 0.5 T (6f + 6^ ) 

( 2 . 30 ) 

where 




^TOm^^bm^^TtiTTi mean value of f^,f|j and fjjj^ 
respectively 

s^,s^^,Sjjj = standard deviation of and f^ 

respectively 

and “ coefficient of variation of 

and fjjj respectively. 

In absence of test results, the strength of masonry can be 
obtained by the follo?nng formula (changed to SI \mits) as 
per the recommendations given by SCPI (I3) 

f = A (2.758 + B f^) (2.51) 

w 

f = assumed compressive strength of masonry at 
w 

28 days in l/mm 

f^ = average compressive strength of brick in 
N/mm^, but not to exceed 96.53 N/mm^ 


where 
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A = a coefficient eq.iial to 2/3 v/ithout inspection 
and 1,0 with inspection 

and B = a coefficient ec[ual to 0,2 for type N mortar j 

0,25 for type S mortar and 0,3 for type M mortar. 

Type N, S and M mortar should be according to the specifications 
given in ASTM 0270-68(95) and the compressive strength of 
mortar corresponding to different types are given below: 


Mortar Type 

Average Comp, strength 
fm 

B 

IT 

750 psi (5.17 Vffim^) 

0.20 

S 

1800 psi (12,41 F/mm^) 

0.25 

M 

2500 psi (17.24 N/mm^) 

0.30 


where 

= average laboratory compressive strength of 2 inch 
cube at 28 days. 

An equation of the foim B = a + pfjj^ is fitted to find the 
relationship between the coefficient B and compressive 
strength of mortar f^* The value of oc and P thus obtained is 

a = 0.155 and P = 0.0082 
The equation relating B and fjjj is then given by 

B = 0.155 + 0.0082 fjQ (2.32) 

Then Bq, 2,31 can be expressed as 

f = A (2.758 + 0.155 fb + 0.0082 f^) 

^ (2.33) 
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The mean, standard deviation and coefficient of variation 
of oan be estimated from Sq.« 2,33 by (see ippendix A ) 


fm “ A(2.758 + 0,155 f^jjjj + 0.0082 f^ 

\ (0.155+0.0082 + (0.0032 


(2.34) 

(2.35) 

( 2 . 36 ) 


2.5.3 Monte farlo simulation 


Monte Carlo simulation tecbniiue is applied to 

generate samples of masonry strength by the deterministic 

foraiula given in Sq. 2.27 and Eq. 2.33. One of the usual 

objectives in using the Monte Carlo simulation technique 

is to estimate certain parameters like mean, standard 

deviation etc. and the probability distributions of random 

variables whose values depends on several other random 

variables of specified distributions, let I be a random 

variable which is a function of several independent random 

variables with probability density fimction fx . • I is 
J 3 

given by 

Y = g (X^, X 2 ,X^, X^) ( 2 . 37 ) 

Monte carlo simulation is carried out in the following way (78) 

i ) A value X^j^ is generated for each random variable X^ 

from their specified distribution, Ex.(x^) or 
probability density fj.(x^) independently. 
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The kvh function value is calculated from the 
3q. 2.37 as 

= g 22Jj, 

and stored in a array. 

iii) Step (i) and (ii) are repeated for k = 1 to SF 
(usually the required number of experiments is 
decided before hand) times. The array of Y vfill 
then contain if simulated samples of Y, 

iv) The generated samples of Y is then used to compute 
the statistical moments and chi-square test is 
conducted to fit a particular distribution to the 
data. 

Strength of brick f^^ and strength of mortar fj^ are assumed 

independent and normally distributed random variables in 

m 

the Monte Carlo simulation. The mean and standard deviation 

of each variables are as follows 

^bm “ 23.69 » % “ 5.439 H/mm^, 6^ = 0,232 

^mm N/mm^, = 2.362 N/mm^, = 0.133 

30,000 samples of strength of masonry f^ are generated by 


simulation from the randomly generated data of and 
using their individual statistical properties as given above. 
Deterministic formulae given by Eqs. 2.27 and 2.33 are 
used, Eig. 2.28 illustrates the histogram and cumulative 
distribution of the simulated masonry strength using Sq. 2.27 


with k^ = 0.4* 


ahe simulated data did not follow normal 



Probability Relative frequency 




Fig 2-28 Histogram ot simulated strength 
of masonry using Eq,2*27 
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distribution even at 0,1 percent significance level, (Etie 
central region of tbe simulated data matcbed well with 
noimal distribution, but mismatching was observed at either 
tails. 

Ihe histogram and cumulative distribution of the 
simulated masonry strength using Eq. 2.55 with A = 1,0 is 
shown in Eig. 2.29. Simulated data based on Eq, 2.35 followed 
normal distribution at 0.5 percent significance level. Mean, 
standard deviation and coefficient of variation computed from 
Eqs. 2.28 to 2,30 and Eq». 2.34 to 2.56 are given in Table 2.18, 
Corresponding values estimated from the simulated data are 
also presented in Table 2. IS, 

2,5.4 Characteristic strength of masonry 


Characteristic strength of masonry can be defined 
as the strength under which the test results falljvfith an 
accepted probability p^ . The relationship between charac- 
teristic strength and mean strength depends on the choice 
of assumed distribution. If masonry strength follows 
normal, then the relationship is given by 


^wh ^wm ^ 


( 2 . 38 ) 


and for f following lognormal distribution is given by 
w 




^ f(62+l) 


^ exp [ k fin (6^+1 ) ] 


(2.39) 


- ^ (Pf) 


where k 



Probability Relative frequency 



3*0 6-4 9*8 13*2 16*6 20*0 

fyy (N/mrr?) 


(a) Histogram 



(b) Cumulative distribution 


Fig.2-29 Histogram of simulated strength 
of masonry using Eq. 2*33 
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lable 2.18 : Mean, Standard Deviation and Coefficient of 
Yariation of Strengtii of Masonry 


Estimated 

Yalues 

f =k V (fuf ) 

f^=A(2. 758+0. 

155fm+O.OOS2 

By 

simulation 

From Eqs. 
2.28 to 
2.30 

By 

simulation 

From Sq.s, 
2.34 to 

2.36 


8.958 

9.036 

10.615 

10.614 

sjN/am^) 

1.259 

1.207 

1.923 

1.903 

6w 

0.1406 

0.1335 

0.1812 

0.1793 


ISiote : 


f ^ - N (23.69, 5.489) 

N/im^ 

, 6j, = 0.232 

-* I (21.54, 2.862) 

N/mm^ 

, = 0.133 


= 0,4 and A = UO 
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The above two equations are rearranged as 
^wk 

= (1 + k 6^) (2.40) 

f 1 

= , 7-0 exp[ k fln(62+i)] (2.41) 

^wm T(6^+1) 

S'ig. 2*30 shows the graph of the above equations for an 
accepted probability of 5 percent, 2.5 percent and 1 percent 
(Pf ~ 0*05 , 0,025 and 0,01) and for coefficient of variation 
^^Si^g from 0 to 0,30, As tho coefficient of variation 
increases, the ratio between characteristic to mean strength 
decreases. Por a specified characteristic strength the 
required mean value predicted by nomal distribution is 
always higher than that predicted by lognormal distribution. 
Thus prediction through Eq, 2,38 of normal distribution 
will be always on safer side as compared to lognormal 
distribution, BS 5628: Part 1 (43) defines characteristic 
strength of masonry at 5 percont risk. IS Code (4I ) and 
GBB-PIP committee (82) recommends the prediction of 
characteristic strength of concrete at 5 percont acceptable 


risk 


(2.42) 


where fjj. = characteristic strength of concrete, 
= mean strength of concrete 
s = standard deviation. 



UJM 



Fig.2-30 Relation between characteristic and mean 
strength of masonry 
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2.6 Discussions and Conclusions 

2.6.1 General 

Brick samples v^ere collected from different brick 
manufacturers of Kanpur zone and Karoke (Kerala) zone to 
study the q.uality of bricks manufactured in the coimtry* 

(The thickness of mortar joint from different existing 
buildings v/as measured and statistical analysis has been 
carried out. Cubes of three different mixes were cast and 
tested to study the variability of strength of mort.ar. 
Masonry strength was simulated in a digital computer 
through different deterministic formula and the probability 
behaviour of masonry strength is studied. Bricks and 
ma,sonry is discussed separately in the following sections. 

2.6.2 Bricks 

The coefficient of variation of brick dimensions of 
different brick manufacturers has been found to be vrithin 
3 percent. Statistical variations of different dimensions 
of precast elements of prestressed concrete T beams under 
controlled conditions were studied by Ranganathan (85). It 
was reported that the coefficient of variation in width 
and ddpth of laboratory made beam sections was within 

1.6 percent whereas that corresponding to the thickness 
and cover was in the range of 5 to 8.6 percent. Hence, it 
can be concluded that the consistency maintained in the 
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dimensions of the hand made bricks is satisfactory and can 
be compared well with that maintained in concrete construction. 

Bricks of Kanpur zone have been found much stronger 
than those of Paroke, 'Ihe coefficient of variation of 
compressive strength of bricks of different manufacturers 
has been found to vary from 15 to 24.6 percent, A similar 
statistical analysis oaanried out on several sots of 
concrete cube samples supplied to the laboratory for quality 
analysis, gave coefficient of variation of 23.9 and 18.8 
percent for all tho random sets of Ml 5 and M20 concrete 
respectively (97). The corresponding values of M25 and 
M30 concrete were found to be 12 and 1 3 percent respectively. 
The coefficient of variation of compressive strergth of 
brick can be compared with that of M15 concrete. Tho cube 
strength results of M15 concrete which wore analysed are of 
the nominal mix propoirtioned concrete. Bricks from different 
manufacturers put together resulted in coefficient of 
variation of strength as high as 29 to 31 percent. The 
same degree of variability is also seen in random lot. 

2.6.3 Masoniy 

The strength of masonry depends on several factors 
such as strength of individual brick unit, strength of 
mortar and thickness of bed joint, the type of supervision 
given in the cor^truction, etc. Effect of uncertainity or 
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variability of different parameters are responsible for the 
variability of masonry strength* Assuming that IK) single 
parameter affects the variability of strength of masonry too 
much, the strength of masonry can be assumed as noiraally 
distributed by the central limit theorem. However, actual 
data of masonry strength from the field is to be collected 
to estimate the variability of masonry strength and its 
pro bability behaviour, 

2,6,4 Conclusions and recommendations 

The following conclusions and recommendations are 
arrived at in the present investigation; 

1 • Dimensions of bricks manufactured from a particular 

manufacturer have a very small coefficient of 
variation ranging from 1 to 5 percent and have tendency 
towards deterministic values. The dimensions of 
brick can be assumed to be normally distributed for 
practical purposes, 

2, Density of dry bricks of a manufacturer follows 
normal distribution and the coefficient of variation 
is found to be 5 percent which reflects a good, 
quality control over the mixing Jind compaction of 
soil in brick making, 

3, Percent of water absorption did not fit normal or 
lognormal distribution and the coefficient of varie.tion 
is found to be anywhere between 13 to 22 percent. 



Bricks of Kanpur zone are mucli stronger than those 
made in some places of Baroke (Kerala). 

Strength of brick follows nomal distribution in 
most cases specially for a particular manufacturer. 

The thickness of mortar ^oint follows nomal 
distribution at 1 percent significance leTsl and 
lognormal at 2.5 percent significance level. 
Compressive strength of mortar (laborato 3 ::y specimens) 
follows normal distribution at 10 percent significance 
level and has a coefficient of variation 
of 10 to 18 percent. The coefficient of variation of 
field cubes is expected to be about 20 pcixent more 
than that observed in the laboratory an d could be 
anywhere between 12 to 22 percent. This is comparable 
with the quality of M15 concrete. Kield de^ta need 
be collected. 

Actual field data on masonry strength need to be 
collected for a better xmderstaiding of the 
statistical behaviour of masonry strength. Simulated 
data based on deterministic formula 

f = k t(fufm) 

'*“'57 'V7 * ^ Dm 

derived from limited experimental investigation did 
not follow normal distribution even at 0.1 percent 
significance level. Simulated data of masonry 
strength based on the deterministic formula 
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^vr “ ■^^2.758+0.155f^+0,0082f|ji^) 
followed normal distribution at 0.5 percent 
significance level. Strength of masonry may be 
assumed as normally distributed random variable. 
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REIiIiBIH!!?? OP GHI-S®IARE 
5,1 Introduction 

load and strength of laaterials like concrete, masonry, 
reinforcement etc. have been recognised as random variables. 
The probability of failure of a structural member under the 
action of external load can be calculated as 

Pf = P ( R < S) (3.1) 

where H and S represent resistance load respectively. 
Probability of failure p^, is sensitive on the choice of 
the probability distribution of R and S. To compute p^ ^the 
data of load S may be collected from a load survey on 
existing buildings and then an empirical probability 
distribution is to be fitted to the data. Similarly, de-ta 
of resistance B, is to be collected and then an en^irical 
probability distribution model is to be fitted. After 
arriving at the probability distribution of the individual 
variable, tne can easily compute the probability p^, by 
Eq. 3.1 assuming R and S are independent raniom variables. 

In most cases, only the data is available from field, 
but the probability distribution, along with the parameters 
of the distribution is unknown to the engineer or the 
statistician. Suppose, data of a rrjidom vexiablo is collected 
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f 2?om the field and one has to fit a probability distribution 
model to the data, Por obvious reasons, one will face the 
following questions: 

(i) Since there are many standard probability 
distributions available, which distribution should 
be tried to fit the data? 

(ii) How the parameters of the distribution could .be 
estimated from the data and what sort of statistical 

tests are to be done to fit a particular probability 
distribution? 

One has to exercise his engineering or scientific 
judgement to select a particular probability distribution 
considering the actual physical phenomena* Two or more 
probability distributions can also be selected as equally 
competent by considering the actual phenomena* For example, 
the strength of masonry is a random variable. !Ilhe variability 
of masonry strength depends on several factors, which are 
themselves random variables. So, if the uncertainity of 
the individual factors are assumed to give an additive 
effect on the randomness of strength of msisonry one c^ 
select normal distribution as a model to fit the masonry 
strength data through central limit theorem. On the otherhand, 
if one assumes that the randomness or variability of masonry 
strength comes from the multiplicative effect of several 
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individual factors, lognormal distribution can also be 
selected as the probability model since log-noiaal 
distribution comes from the multiplicative effect of 
individual random variables. Ihus, normal and lognormal 
distributions are the two competent distributions to fit 
the data as far as the engineering judgement is concerned. 

After arriving at the decision regarding the 
probability distribution one has to test whether the 
distribution can be accepted as a mathematical model to 
represent the particular random variable. If there are more 
than one competent distributions, decision has to be taken 
as to which one of them represents the data in a better way, 
This leads to non-parametric hypothesis testing problems. 
Several non-parametric methods are available in the 
literature for testing hypothesis. Chi-square test and 
Kolmogorov-Smirnov test (K-S test) are the two popular 
tests used in hypothesis testing, although both of th^i 
suffer from some disai vantages. 

The bhi -square statistic computed in the chi-square 
test depends on the way how the groups are formed and the 
niBuber of class intervals chosen etc. Thus, from the same 
data, one can accept a particular null hypothesis at a 
particular significance level whereas another parson by 
grouping the data into different classes m^jr reject the null 



105 


li3T)othesis at that particular significance level. Some of 
these problems of chi-stjuare test have been highli^ted in 
this chapter and a practical recommendation is given for 
conducting this test. Kolmogorov-Smirnov test does not 
group the data like «hi— square test and axbitary decisions 
like number of classes, width of class etc. are not involved 
in this teat. Hence, from a particular data and a particular 
null hypothesis, the decision of acceptance o.r rejection of 
null hypothesis at a particular significance level is 
unique and does not depend on the person concerned* The 
standard tables used for K - S test are valid when the 
null distribution is completely specified but if one or 
more parameters have to be estimated from the data then 
the tables are no longer valid. Tables are given for the 
Kolmogorov-Smirnov statistic for different null distributions 
(such as uniform, exponential, normal) v/hen the parameters 
are not specified but have to be estimated from the data. 
These tables are generated by Monte Carlo simulation and 
•are presented in this chapter. 

3,2 Hypothesis Testing and Significance level 

let “ independent observations of a 

random variable X with unknown probability distribution 
function To test whether the observations come from a 

population following probability distribution function Jq, 



a null liypothesis that the random variable follows a 
probability distribution is set up against an alternative 
hypothesis that it does not follow Fq. Mathematically it 


can be expressed as (93,98) 
Ho ’ = ^0 


: F^ = not F^ 

A statistic D is then formed under the null 
hypothesis !I?he null hypothesis is accepted at a 

patrticular level of significance a if the observed value 
of the statistic D is less than the critical value If 

the distribution function of D is known then the critical 
value Dq; can be found as 


P (B > ) = a 

or 

where 5*5^ (•) = inverse distribution function of D. 

Two types of errors are involved in hypothesis 
testing and they are called type I and type II errors. 
They are as follows J 

Type I : One may wrongly reject when it is true. 

Type II i One may wrongly accept \’\hen it is false* 

The probability of committing a type I errer is 
called the level of significance u (also called size of 
the test), and the probability of coumiitting a type II 
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error is a fanction of ttie alteimative hypothesis and is 
usually denoted by p. (l-p) is called the power of the 
test of the hypothesis against the alternative hypothesis 
• The complete specification of is essential, since 
power is a function of 

Usually, a fixed value of a is chosen ( a = 0.01, 0.05, 
or 0.10) before hand and the test is conducted. On the 
otherhand, the maximum significance level at which 

the null hypothesis can be accepted, may be computed without 
prior choice of a. Then, the null hypothesis at any 
significance level less than can be accepted. The 

maximum significance level at which the null hypothesis 

can be accepted is called p-level and can be computed as 

P = w = ^ ® ^ 

where is the observed value of the statistic D. 

3.3 Reliability of Chi-square Test 
3.3.1 Chi-s quare test 

In the usual Pearson chi-square test, the observations 
are grouped into k classes and the chi-square statistic is 
computed as 



i=1 % 


where 

— Chi— square statistic 
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k number of classes 
0^ = observed frequency of the ith class 
= expected frequency of the ith class 

corresponding to null hypothesis HQ(i.e^l‘Q), 

The statistic has a distribution with k -1 
degrees of freedom as n becomes large (96), When the 
parameters of the null hypothesis Pq are estimated from the 
data by the maximum likelihood method, the above statistic 
is distributed as x with v degrees of freedom, "vdiere v 
is defined as 


V = k - r - 1 (3.5) 

where r is "^he number of parameters estimated from the data, 
♦ 

The chi-square method gives a numerical measure of 

the difference between the observed and the theoretical 

distribution, let denotes the observed value of X^ • 

0 

For a given significance level a , the null hypothesis is 
accepted if 

whore X'^ is defined by the following equation 


,2 

‘a,v 


P (X^ > ) = a 


(3.6) 


The maximum significance level at i^ich the null hypothesis 
can be accepted is computed as 

‘ P = 




(5.7) 
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The computed value of dhi-square statistic depends on the 
number and width of classes. Care has to be taken that 
^ 3 for all classes. 

3»3»2 Choice of class interval 

The way the data is grouped to form the classes 
plays an important role in the 6hi— square test. IPwo schoaes 
are generally adopted in practice for the choice of class 
Intemrals: 

(i) equal length of interval 
(ii) equi-probable cells 

Stiruges ( 99 ) suggested a formula for the choice of claiss 
interval which can be computed as 

C = — (3.8) 

1 + 3.222 log^Q n 

where 

C = class interval 
R = the range of the obseirvations 
and n = number of observations. 

Benjamin and Cornell (93) also reconmended the same 
foimiula given by Struges (99) ns a practical guide to choose 
the number of classes (k).for ohi-square test and is given by 

= 1+3.2 log |0 ^ 


k 


(3.9) 
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Mann and Wald (100) proposed a formula for the 
optimal choice of the number of classes for chi-square test 
with equi -probable cells where the population parameters 
are not estimated from the data but are based on staMard 
theory or past experience* The number of classes can be 
calculated as 



( 5 . 10 ) 


where c is determined such that 

1 “ - ^/Z 

— / e dx = a 
YTn c 

where a is the level of significance. 


(3.11) 


The usual chi-square test is unreliable (101 ) for 
a continuous variate when equal width classes are formed, 
since it involves three arbitary decisions such as starting 
point, width and number of classes. Gumbell (101 ) has 
shown through examples that from the same observations, 
different statisticians, equally well trained and equally 
careful may obtain different probabilities of the critical 
region (i.e. different p— level) and may proclaim anyone of 
these as final. Thus, chi-sq\iare test does not lead to a 
decision whether a h 3 rpothesis has to be rejected or not. 
When equi— probable cells are used in the chi-square test, 
the decision about the starting point and width of classes 
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does not arise, Williams (102) suggested that the number of 
classes given by Marni and Wald (100) may be halved for 
practical purposes without greatly affecting the power of 
the test, 

Watson (103) recommended that at least 10 classes 
of not very unequal probability content should be used* 

Hamdan (104) showed that a number of classes between 10 to 
20 is adequate to ensure a reasonably sensitive test, and 
the optimum choice of class boundaries corresponds to eqtial 
class width of about 0,4 standard deviation being sli^tly 
more pov/erful than the equal class probability partition. 

The paper gives a direct approach to the problem of boundary 
determination in the case of the ohi -square test for location 
of the normal distribution, Dahiya et al,(l05) discussed 
about the number of classes to be used in testing for 
normality against certain families of alternatives. They 
recommended a range of choice of k for several different 
alternatives in testing for normality, 

3.3.3 Variation of significance level 

Chi-square test when used with equal width classes 
has some limitations as it depends upon three arbitrary 
decisions discussed earlier. If a certain choice of the 
intervals gives a good fit, it cannot be concluded that a 
broader or narrower classification gives the same or better 
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fit. It may so happen that a different classification may 
lead to rejection of the null h3rpothesis« Observed value of 
Chi-square may change considerably through different sets 
of starting point, although number of degrees of freedom 
remains constant, Ihe arbitrary choices of the three 
variables can be made to a single variable, if cells of 
equal probability are used instead of cells of equal width. 
Best choice of the number of classes fc, given by Mann a-nr) 
Wald (lOO) is valid for small levels of significance and for 
large number of observations, but the question of the best 
choice of k for small number of observations and largo levels 
of significance is not yet solved, Uumber of classes 
needed for dhi-square test as computed from 5,9 and 
3.10 for different sample sizes are shown in fable 3.1, 

Data of strength of mortar of three different mixes 
are taken up for numerical experimentation to study the 
effect of number of classes on the chi-square test result. 

The chi-square test is carried out using cells of equal 
probability for different number of classes. Ifumber of 
classes considered here are 5 to 40. i. t 3 ?pical variation 
of chi— square statistic with number of classes for the data 
of 1 :5 mix is shown in Fig. 3.1(a), when the null hypothesis 
is normal distribution, and the corresponding p— li^vcls at 
which the normal distribution can be accepted is shown 
in Fig. 3.1(b). It can be seen from Fig. 3.1(a) that there 
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Table 3.1 

J Sf-umber of ■ 

Classes 

for the Chi-square 

Test 

n 

Prom 

Eq.3.9 


From Eq. 3.10 



level of significarice 



0.01 

0.05 

0.10 

50 

7 

16 

18 

20 

100 

8 

21 

24 

26 

150* 

8 

24 

28 

31 

200 

9 

27 

31 

35 

500 

10 

40 

45 

50 

1000 

11 

52 

60 

66 

2000 

12 

69 

79 

87 

5000 

13 

99 

114 

126 

10000 

14 

131 

150 

166 


-level (7.) 


30i 


114 



No. of classes 

(a) Observed values of statistic 



(b) Maximum significance level 
Fig. 3-1 Variation of chi-square and p-level w.th 
number of classes for normal distribution 
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is no steady increase or decrease of chi-square value with 
the increase in number of classes* A random fluctuation 
in p— level is observed* IThe p— level varied from 26 to 99 
percent as can be seen from Fig, 3.1 (b). Similar fluctuations 
in p— levels are found in other two sets of data (1 j 4 auri 1 :5 mix). 
Mean, standard deviation, coefficient of variation and range 
of p-level for three sets of data are given in Table 3 . 2 , 

The maximum significance level at which null hypothesis 
can be accepted is not only a function of the data but also 
a function of the number of classes chosen. Sometimes a 
drastic change in p-level between any two consecutive number 
of classes is also observed. For example, if lc =16 is used, 
the p-level becomes 93.5 % which signifies a rery good fit 
of normal distribution to 1 53 mix data • On the other hand, 
if k =17 is used, the p-level becomes 64 . 5 %, Similar 
variations are observed in all the cases which clearly 
show that the test should not be done by just taking a 
particular value of k and come to a final decision about the 
acceptance or rejection of the null hsrpothesis. 

If the null hypothesis is assumed to be lognormal, 
the value of also changes with the choice of number of 
classes. The variation of and p-level with the number 
of classes for 1 ;3 mix data is shown in Fig, 5 * 2 * The p-level 
varied from 46*7 %to 99*7% with a mean value of 82*8 5^ 

Mean, standard deviation, coefficient of variation and range 
of p-level for three different mixes are given in Table 3 . 3 * 



Table 3.2 ; 


P— level ( %) for No anal Distribution 


Sample 

designation 

Mean 

S.D. 

C.O.Y. 

Rar^e 




lower 

Upper 

1:3 

mix 

78.316 

16.746 

0.2138 

26.15 

98.96 

1 :4 

21 . 608 

17.987 

0.8324 

0.92 

67.76 

mix 






1:5 

mix 

11,302 

12.467 

1.1031 

0*60 

42.89 


Table 3.3 J p-levels (%) for lognormal Distribution 


Sample 

designation 

Mean 

S •!)• 

a.p.T. 

Range 




lower 

Upper 

1:3 

mix 

82,839 

14.584 

0.1760 

46.73 

99.72 

1 :4 
mix 

30.702 

22.583 

0.7356 

0.68 

75.54 

1 :5 
mix 

8.898 

10.621 

1.1936 

0.25 

56.32 




Variation of chi-sc|uar€ and p-lsv€l with 
number of classes for lognormal distribution 
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3.3.4 Characteristic significance level 

As discussed earlier, it is found that the maximum 
significance level varies randomly with number of classes 
used in chi-square test. Choice of different value of k can 
lead to random results. In order to overcome this difficulty, 
a series of chi-square tests with different values of k may 
be conducted to decide the most reliable p-level. 

The p-levels observed for different values of k for 
different mixes did not follow any particular distribution, 
althou^ a certain grouping tendency has been observed. 

Let pj^ be defined as characteristic p-level such that the 
probability of p-level falling below pjj. is p* , Assuming p 
follows normal distribution, the characteristic p-level 
Pj^, can be computed as 

Pk = Pm ^ h ^ 


v/here Pjji = mean of the p-levels 

s = standard deviation of the p-levels 

k^ = 0"'’ (pp 

and 0“^(*) is the inverse of the cumulative standardized 
normal distribution function. 

If the probability of falling below p^ is accepted to 
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be 5 percent, the characteristic p-level can be computed as 

Pk “ •" s = pjj^ (1-1.645 5) (3.13) 

where 6 = coefficient of variation 

k^ = j2^‘‘’'(0.05) = -1.645 

Coefficient of variation of p— level is found to be 
as high as 120 % in some cases as shown in Tables 3* 2 and 3.3. 
Characteristic p-level computed from Sq. 3.13 will become 
negative for any value of 6 greater than 60.79 percent. 

Theoretically, p-level cannot be negative and at 
the most can vary within 0 to 100 %. The theoretical or 
conceptual difficulty of pjj. being negative can be avoided 
it lognormal distribution is assumed. It should be noted 
that the limiting value of p-level is 100 % whereas lognormal 
distribution permits value of p-level greater than 100 %, 
However, lognormal distribution can still be used to compute 
characteristic p-level for practical, purposes. Proscribing 
lognormal distribution for p-level, the characteristic 
p-level can be computed as i 

Pk = ^1 ] (3.14) 

where 6 = coefficient of variation 

k^ = 0”'' (pp 

For an accepted probability of 5 %, 3.14 becomes 

T3, = e:^ [ -1,645 TClnd+d^))] (3.15) 

tci+ 62) r 



Hhe characteristic p-level defined by Eg., 3*15 
reflects the random variation of p-level through mean and 
coefficient of variation. The characteristic p-levels 
computed by Sq. 3.15 for the three sets of data are given 
in Table 3.4 for nomal and lognomal distributions, 

3 . 3 « 5 . Acceptance of the better fitted distribution 

'When there are more than one competent distribution, 
the problom often encountered is how to judge which 
distribution is a better fit to the data. For example, frcmi 
ei^ineering judgement both normal and lognormal distributions 
can be taken as the representative model of strength of 
mortar. It has to be decided which model is better. 

For the same number of classes, values can be 
computed for normal and lognormal distributions and the 
distribution which gives a less value of chi-square can be 
taken as a better fit compared to the other competent 
distribution. The computed value of ^^2 is found to be 

22.05 and 15,8 for normal and lognormal distributions from 

Figs, 5,1 3,2, when the number of classes k = 3I is used. 

This leads to a conclusion that lognormal distribution is 

O 

a better fit. On the other hand, values of X are 16,4 and 

22.5 for normal and lognormal distributions respectively 
for k = 32, wdiich leads to the opposite conclusion that 
normal distribution is a better fit. Thus, direct comparison 
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Table 3.4 : Characteristic p~level (%) 



mix 

1:4 21,608 17.98? 5.05 30.702 22.585 8,38 

mix 


1:5 11.302 12.467 1*75 8.398 10,621 1.22 

mix 
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of chi sq^uaz*© at a particular chosen nuEiber of classes may 
not lead to a proper conclusion, She p-leirels for nomal 
and lognormal distributions for different values of k in 
the case of data of 1:3 mix are shown in Pig. 3.3(a). 

The firm line shows the variation of p— levels for normal 
^distribution whereas the dotted line corresponds to 
lognormal distribution, !Ehe p-level of lognormal distribution 
is higher than that of normal distribution f fov k = 5 to 31 • 
Thus, if the numoer of classes is anyvdiere between 5 and 31 , 
lognormal distribution can be claimed to be a better fit. 

On the other hand, if the test is caarcied out with number 
of classes anywhere between 32 to 40, normal distribution 
appears to be better fit over lognormal distribution. Here , 
two distinct zones are observed, Por k varying from 5 to 31* 
p-level s for normal distribution are alv/ays higher than 
those of lognormal and for k varying from 32 to 40, p-levols 
for normal distribution are always lower than those of 
lognormal distribution as seen in Pig, 3.3(a). 

In case of data of 1 :4 mix, no such distinct zones 
are observed. The variation of p-lovcl for normal and 
lognormal distributions is shovm in Pig, 3.3(b), The 
p-levels of normal and lognormal fit are 7,5%“ and 26%for 
k = 10, which leads to the acceptance of lognormal 
distribution as a better fit. But, for k — 11, p— levels 
are 26 % and 16 %for normal and lognormal distributions 



(•/.) p- level (•/.) ' p>(cve( (*/•) 


log norma I 


normal 


* 1 i t 1 1 1. 

5 20 

,1 

25 30 

XJ 1 1,,,l 1 i 

35 

No« of 

dosses 


(q) 1:3 

mix 



o lognormal 

normal 


•m *XH*7*1 


kcli 


No. of classes 

(b) 1*.4 mix 


normal 


lognormal 


No. of dosses 

(c) V.5 mix 

Fig.3‘3 Variation of p4€/el with number of 
classes for normal ana lognormal 
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i^especliively, wliich. leads "fco tiie acceptance of as 

better fit* A change in the choice of k frcas 10 to 11 has 
changed the decision about the acceptance of the model* 
Similarly, k = 12 will again change the decision iidiich 
has been foimed by taking k = 11. fhe decision about the 
better fit of a particular distribution keeps on changing 
randomly with the choice of number of classes k as can be 
seen from Fig. 3.3(b), 

Data of 1 :5 mix has shown a similar type of behaTiour. 
The variation of p-level with number of classes is shown 
in Fig, 3, 3(c), In some regions, normal distribution 
appears a better fit as compeired to the lognormal distribution. 
The variation of p-level is snail in the range of k varying 
from 7 to 23, but p-level is found to vary drastical l y for 
k varying from 23 to 40. 

In the three typical c^es discussed above, it has 
been observed that the p-level for a particular distribution 
depends on the particular choice of the value of fc and the 
decision that a particular distribution is better than the 
other distribution also depends on the choice of k. Thus, 
acceptance of a particular distribution by chi-square test 
Is highly dependent on the choice of k. 

A qualitative analysis of the reliability of chi- 
squarc test has been carried out in the present investigation. 



125 


It is clear that one should not take a particular value of 
number of classes and carryout the ohi-square test to d-ecide 
which one of the two competing distributions is a better 
fit, fo judge which distribution is a better fit, a 
practical method is described as follows, First, the 
characteristic p-levels for the two competing distributions 
are to be found out after performing a series of cSii-square 
tests with different values of k and then tho distribution 
which gives a hi^er characteristic p— level should be 
accepted as a better representative model. Comparing the 
characteristic p-levels for normal and lognormal 
distributions given in Table 5*4 » it can be concluded that 
lognormal distribution is better fit for the data of 1 :5 
and 1 :4 mixes, whereas normal distribution is better for 
the data of 1 :5 mix. 

3 .4 Kohnogoro v-Saii 2 ?nov Tost 

Kolmogorov-Smimov test (K-S test) is an alternative 
to chi-sq.iiare test. In this test, a statistic D, defined 
by the following formula »is computed to test whether the 
sample belongs to a completely specified distribution; 

D = maxi 5 * 0 ( 2 ) -- Sj^(x)| (3.1^) 

X 

where »F«(x) is the ccmpletely specified continuous 
o 

cumulative distribution function and Sjj^(x) is the sample 
cumulative distribution function. She statistic D is a 
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distribution free statistic and depends only on tbe sample 
size n (106,107)« llthough. both chi-s«iuare and K— S tests 
suffer from some disadvantages, K-S test iias the following 
major advantages over the usual chi-square test(93,96,98,107) : 

(a) for small sample sizes, the validity of chi-square 
test becomes questionable vdiereas K-S test can be 
used with small sample size, 

(b) often it appeeirs to be a more powerful test than the 
chi-square test for any sample size, 

(c) data is not grouped into classes and hence choice of 
class intervals etc., does not complicate the 
procedure, 

(d) this test considers the effect of all individual 
data points in the sample where as the individual 
effect of every data point is lost by grouping 
the data into classes in chi-square test, 

Smirnov (106) has given a table of critical values 
of statistic 1 for use with K-S test, Massey (107) has 
given a table of the statistic h for different sample 
sizes n = 1 to 35 for levels of significance 1,5,10,15 and 
20 percent. A more extensive table was given by Miller(108). 

ill these tables are valid only vfeen the null 
distribution is completely specified. When the parameters 
have to be estimated from the sample, then the statistic I) 



IS 210 longer a distribution free statistic and the tables 
are no longer valid. If the test is used with the paraoieters 
esbimated from the sample > the critical values of 3D are to 
be decreased than those given in the standard tables (106, 107, 
108) as suggested by Massey (107). Idlliefors (109) has given 
a table to be used with the Kolmogorov-Smirnov statistic for 
testing whether a set of obseirvatiotis is fi*oni a nonaal 
population when the mean and the standard deviation (obtained 
by n-1 as the denominator) ace estimated fixim the data. Ihe 
table was obtained by Monte Carlo simulation and is given for 
1» 5, 10, 15 and 20 percent significance levels. In the 
present work, similar tables are generated by Monte Carlo 
simulation for the follov/ing distributions for all significance 
levels using similar procedure as given by lillief ors(109) 
for normal distribution: 

(i) uniform distribution 

(ii) normal distribution 

(iii) exponential distribution. 

The procedure is illustrated in the following steps: 

(i) Generate n number of values frcmi a particular 
distribution F(x), with some values assigned to 
the parameters.- 

(ii) Estimate the parameters from the sample of size n. 

(iii) Compute the statistic 

D = max (F*(-x) - S^(:k)\ 

where, F*(x) is the cumulative distribution function 
i*(x) with parameters estimated frxmi the data and 
(x ) is the sample distribution fumjtion. 
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Repeat steps (i) to (iii) if times (1000 or more) 
ajid store the values of B in an array. 

3?he distribution of the statistic B is then estimated 
and the critical values Bcc different significance 
level are picked up such that 
P (B y Bq^ ) = tt , 

U niform distribution: 

Ihe parameters of uniform distribution for eve 2 :y 
sample are taken as the leftmost and ri^tmost values 
respectively, !Ehe critical values of B are given in 
Table 3.5 for different sample size. 

Uormal distribution; 

The parameters of the normal distribution are 
estimated as = 3^ and o = s, standard deviation with 
denominator n— 1, The critical values of 3) for different 
sample size are given in Table 3.6. These values mateii 
well with the values given by Lilliefor3(109), although 
there are differences observed in the third decimal place, 

Exponential distribution; 

The parameter of exponential distribution is 
estimated from the mean of the sanple* The critical values 
of B for exponential distribution are given in Table 3.7. 


(iv) 

(v) 



Table 3.5 


I Critical Talues of D for Uniform 
Distribution 


level of 
signifi- 



•Saaple size 

n 



cance a 

4 

6 

10 

20 

25 

30 

over 30 

0.01 

0.669 

0.572 

0.454 

0.339 

0.301 

0.278 

1 . 52 ^fn 

0.05 

0.584 

0.487 

0.403 

0.285 

0.260 

0 

• 

0 

1,281/fn 

0.10 

0.532 

0.449 

0.368 

0.259 

0.235 

0,216 

1,178/fn 

0.15 

0.496 

0.427 

0.359 

0.244 

0.221 

0.202 

1 . 104 /fn 

0.20 

0.483 

0.407 

0.324 

0.232 

0.209 

0.192 

1 . 049 /fn 

0,25 

0.471 

0.387 

0.309 

0.220 

0.200 

0.184 

1 . 007 /Tn 

0,30 

0.459 

0.368 

0.294 

0.210 

0.194 

0.175 

0.955/fn 

0.35 

0.442 

0.348 

0.280 

0.200 

0.186 

0.169 

0.916/fn 

0.40 

0.425 

0.332 

0,270 

0.192 

0.177 

0.162 

0.876/fn 

0.45 

0.414 

0.324 

0.262 

0.185 

0.169 

0.156 

0.842/fn 

0.50 

0.395 

0.317 

0.255 

0.178 

0.163 

0.152 

0,807/Vn 

0.55 

0.376 

0.309 

0,244 

0i171 

0.157 

0.147 

0.779/fn 

0.60 

0.363 

0.301 

0.235 

0.164 

0.152 

0.141 

0.751/fn 

0.65 

0.347 

0.292 

0.227 

0.160 

0.146 

0.135 

0 . 724 /Yn 

0.70 

0.331 

0.279 

0.218 

0.153 

0.141 

0.130 

0.696/fn 

0.75 

0.315 

0.268 

0.209 

0.146 

0.135 

0.125 

0.671/fn 

0.80 

0.290 

0. 256 

0.197 

0.141 

0,128 

0.118 

0.635/fn 

0.85 

0.265 

0.240 

0.185 

0.132 

0.118 

0.110 

0,594/fn 

0.90 

0.250 

0.219 

0.175 

0.122 

0.110 

0.100 

0,547/Yn 

0.95 

0.250 

0.196 

0.158 

0.111 

0.102 

0.091 

0.489/Yn 

0.99 

0.250 

0.167 

0.126 

0.095 

0.087 

0.080 

0.386/fn 


lable 3*6 . Critical falues of D for Homial Distribution 


IieTel of 
signifi- 
cance a 


Sample si 3 

e n 

— 1 



4 

6 

10 

20 

25 

30 

over 30 

0.01 

0.415 

0.374 

0.299 

0.226 

0.201 

0.137 

1.017/fn 

0.05 

0.376 

0.326 

0.261 

0.194 

0.172 

0.157 

0,890/fn 

0.10 

0.351 

0.299 

0.238 

0.176 

0.160 

0.147 

0.807/Y'ii 

0.15 

0.326 

0.282 

0.226 

0.168 

0.151 

0,138 

0.748/fn 

0.20 

0.307 

0.268 

0,216 

0.161 

0.143 

0.131 

0.721/Yn 

0.25 

0.294 

0.259 

0.206 

0.153 

0.137 

0.126 

0,687/Yn 

0.30 

0.285 

0.249 

0.200 

0.146 

0.132 

0.120 

0.656/fn 

0.35 

0.277 

0.240 

0.195 

0.142 

0.128 

0.117 

0.635/fn 

0.40 

0.271 

0.233 

0.188 

0.138 

0.123 

0.113 

0.618/Yn 

0.45 

0.263 

0.225 

0.181 

0.134 

0.120 

0.110 

0.599/fii 

0.50 

0.258 

0.217 

0.176 

0.130 

0.116 

0.107 

0.580/fn 

0.55 

0.252 

0.210 

0.171 

0.126 

0.113 

0.103 

0.565/fn 

0.60 

0.245 

0.204 

0.164 

0.122 

0.109 

0.100 

0.551/Yii 

0,65 

0.239 

0.197 

0.158 

0.118 

0.106 

0.097 

0.536/Yn 

0.70 

0.230 

0.191 

0.154 

0.114 

0 

• 

0 

0.094 

0.514/fn 

0.75 

0.222 

0.186 

0.149 

0.110 

0.099 

0.091 

0.501/fn 

0,80 

0.212 

0.178 

0.143 

0.107 

0.094 

0.088 

0.482/fn 

0.85 

0.200 

0.169 

0.133 

0.102 

0.090 

0.084 

0.458/Yn 

0.90 

0.187 

0,163 

0.127 

0.096 

0.086 

0.080 

0.434/Yn 

0.95 

0.168 

0.150 

0.119 

0.089 

0.079 

0.073 

0.395/fn 

0.99 

0.144 

0.126 

0.103 

0.076 

0.068 

0.062 

0.348/irn 



Table 3 *7 


Critical Values of D for Sxpcnential 
Distribution 


Level of Sample size n 

Sianifi- 

cance a 4 6 


0.01 

0.548 

0.471 

0.589 

0.05 

0.482 

O.4O6 

0.322 

0.10 

0.441 

0.377 

0.297 

0,15 

0.41 8 

0.353 

0.279 

0.20 

0.402 

0.334 

0.263 

0.25 

0.386 

0.318 

0.252 

0.30 

0.371 

0.303 

0.241 

0.35 

0.355 

0.294 

0.254 

0,40 

0.340 

0.283 

0.225 

0.45 

0.327 

0.276 

0.216 

0.50 

0.513 

0.265 

0.209 

0.55 

0.303 

0.256 

0.201 

0.60 

0,295 

0.249 

0.195 

0.65 

0,286 

0.241 

0.189 

0.70 

0.279 

0.231 

0.181 

0.75 

0.267 

0.222 

0.175 

0.80 

0.258 

0.213 

0,166 

0.85 

0.244 

0.206 

0.156 

0.90 

0.230 

0.194 

0.146 

0.95 

0.207 

0.174 

0.134 

0.99 

0.172 

0.149 

0.113 


0,27 3 0.25 2 0.229 1.240yYn 

0.231 0.212 0.195 1.07l/fn 
0.213 0.191 C.178 0.979A^»- 

0.200 0.1 80 0,167 0.91 3/fn 

0.137 0.170 0.157 0.863/fn 

0.177 0.163 0.150 0.832/fn 

0.171 0.155 0.144 O.SOl/Vn 

0.166 0.150 0.140 0.769/fn 
0.159 0.145 0.135 0.743/Yii 

0.154 0.140 0.130 0.717/Yn 
0,150 0,156 0.126 0,689/ir3a 
0.145 0.133 0.122 0.665/fn 
0,141 0.128 0,118 0.645/Yn 
0.136 0.125 0.114 0.619/fn 
0.132 0.119 0.109 0.596/Yn 

0.128 0.115 C.IO5 0,574/Y^ 
0.123 0,110 0,100 0,546A^^ 
0.116 0.105 0.095 0.516/Yn 

0.109 0.099 0.089 0.490A’n 

0.099 0.090 0.081 0.444A'’3a 
0.085 0.071 0.071 0,380/Yn 
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It is observed that the distribution of D is 
dependent on the choice of the null distribution when the 
parameters are estimated from the sample itself, but independent 
of the numerical values of the parameters when the parameters 
are location or scale parameters of the distribution, Ihe 
distribution of D with different levels of significance for 
the three null distributions is shown in Fig. 3,4* These 
curves are drawn for n = 40, since for n > 30, fn D is foiind 
to be almost stable. Bie critical values obtained from 
Smirnov’s Table (106) when the null distribution is completely 
specified, is also shown in Fig, 3»4* It can be noticed 
that the nature of the curve changes with the type of null 
distribution when the parameters are estimated from the 
data and the D values thus obtained are always smaller than 
those given by Sminiov(106), Whon the parameters estimated 
from the data are other than parameters of location or 
scale, the distribution of D is found to be a function of 
the numerical values of the parameters as well as the t 3 ?pe 
of the distribution. This pheiujmenon is observed for some 
of the distributions such as Wcibull, beta, gamma etc. 

The critical values given for normal distribution oan 
be used for lognoimal distribution, After taking natural 
logarithm of the sample data points, the K-S test can be 

p erf oimed for noimal distribution. 



80 95 99 


1 5 to 20 40 60 

Signlflcanc* Uvtl oc {*/•) 

F}g.3*4 Critical values for use with K-5 test 
when the parameters are estimated 

from the sample 
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3*5 G-oodness of ?it 

To test whether a given sample of n observations 

follows an assumed distribution, the E-S statistic D is 

computed and then compared with the critical values given 

in the tables. For example, to test v\hether the data of 

1:3 mix follows normal distribution, parameters of the normal 

distribution are estimated from the data, E-S test, yields the 

value of statistic D = 0.0366. From Fig, 3*4 it is found that 

with = 0*4744 , normal distribution can be accepted at 

a maximum significance level of 82 percent. It is to be 

noted that the p-level found by this test is unique for 

a particular data and for a particular null distribution as 

2 

opposed to the X goodness of fit test vihich changes with 
the number of classes. One of the major advantages of 
E-S test is that there is no ambiguity about the p-level 
of the test, IBo test yhether the same data follows 
lognormal distribution, mean and standard deviation are 
estimated after taking natural logarithm of the data points. 

E-S test for normal distribution is then perfor-’ied giving 
the value of statistic D = O.O 37 I. Hius lognormal distribution 
can be accepted at a maximum significance level of 
80,5 percent. The- E-S test results of the three 
different mix data are given in Table 3«8. Nov/, comparing 


Table 5*8 * K-S Test Result of Three Mix Data of Mortar Strength 
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tile p~levels, the distribution vdiich gives a higher p-level 
can be accepted as a better fit. It can be noticed from 
Sable 3,8 that nonnal distribution appears to be a better 
fit to the 1 :3 mix data whereas lognonaal distribution 
appears to be better in case of 1:4 and 1 j5 thttc data. Had 

the parameters been completely specified or known frcHn 

pas't experience or from common sense, the direct ojmparison 
of D values for the two- distributions could bo done and 
the distribution which gives a lower value of D would be 
the better fit over the other distribution. It can be 
observed from Tables 3*4 and 3.8 that lognormal distribution 
appears to be better fit according to chi-square test 
whereas normal distribution appears to be better fit 
accoMing to K-S test on 1 :3 mix data. The decision 
about the acceptance of the better fitted distribution 
has changed in this case. Similar result is observed 
for 1 :5 data. However, both chi-square and K-S test 
have shown that lognormal distribution is a better model 
to represent 1 !4 mix data than normal distribution. 



CHiPISR 4 


reliability analysis 0? RSI15E0RGSD BRICK BEAMS 
4.1 Introduction 

Strength of masonry (f^) is a r.andom variable 
since it depends on other variables such as brick and mortar 
strength, thickness of mortar joint, joint layout etc, as 
discussed in Chapter 2, Macchi (37) adopted noimal 
distribution for strength of masonry. It was reported by 
Foster (3) that truncated normal and lognormal distributions 
were argued by Beech as better models to represent masonry 
strength. From the statistical analysis of simulated 
samples presented in Chapter 2, strength of masonry is taken 
as normally distributed random variable -in computation of 

probability of failure. Statistical analysis of strength 
of high yield strength deformed bars are presented in 
the following section. 

Reliability analysis of reinforced brick beam 
section at ultimate flexural strength is presented in this 
chapter. Moment capacity of a RBB section is a function of 
strengthsof masonry and steel, area of tensile steel and 
geometric properties of the section. Strength of masonry (f^) 
and strength of steel (fy) are considered as random variables 
whereas rest of the parameters are tireated deterministically. 
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A general formulation for computation of pi^obability of 
failure of EBB section for deterministic and probabilistic 
load is presented considering under— reinforced and over- 
reinforced failures, A Monte Carlo simulation approach is 
also presented to study the probability distribution of 
moment capacity of EBB section for probabilistic variations 
of strengths of materials. Probability density function 
f^Cx) and probability distribution function are 

denoted as fx^^) respectively to avoid mixing up 

of different symbols through out the thesis, i.e., 
differentiation between random variable X and its state 
variable x is not done. 

4.2 Statistical Analysis of Strength of Steel 

High yield strength deformed (HYSP) bars of 8 mm 
diameter were selected to study the variability of ultimate 
strength. Samples were cut from two different lots of 8 mm 
dia HYSD bars and tested for ultimate strength in Universal 
Jesting Machine. Samples of two different lots denoted as 
S1 and S2 can be taken as the representative of two fixed 
lots. Statistical analysis of data of tiltimate strei^h f-y^ 
is given in Jable 4.1. After combining the strength data 
of St and S2 (denoted as S3), statistical analysis is carried 
out and given in Jable 4.1 •’ Sanple S5 can be taken as 
the representative of mixed lot supplied by different 
contractors at a project site. Many builders and engineering 



Table 4.1 s Statistical Analysis of Strength of HYSD Bars 
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supervisors had sent sets of about 3 to 6 HYSD bars of 
diameter varying from 8 mm to 28 mm for qualit;r analysis 
at I.I.T. Kanpur. The sample containing these sets can 
be considered completely random and is denoted as random 
lot. Results of statistical aialysis of ultimate strength 
f.^, and proof strength fy of random lot are given in 
fable 4.1. 

fhe coefficient of variation of ultimate strength is 
found to be v/ithin 1,68 percent for two fixed lots which 
signifies a very good quality cont 2 X)l in a particular 
production batch, fhe mean values of the two lots are found 
to be 558,71 N/mm^ and 606.28 H/mm^. A higher coefficient 
of variation of 4.38 percent is observed in the mixed lot. 
This again reflects a consistent quality control between 
different batches. Coefficient of variation observed in 
the random lot is 12.28 percent for the ultimate strength 
data and it is •"'qual to 11,42 percent for the proof strength* 
Histograms of ultimate strength and proof strength of HIBD 
bars of random lot are shown in Rig. 4.1 and Rig. 4*2 
respectively. Chi-square test is Cvarried out to the 
strength data. Ultimate strength of steel followed nomal 
distribution at 0,5 percent significance level 'sdiile normal 
distribution can be accepted for proof strength of steel 
at 1 percent sigiiificance level. In both cases normal 
distribution appears as better representative model than 
lognormal distribution. Ihe coefficient of variation of 


Fr*qu«ncy 


U.f 



F(g.4*t Histogram of ultimate strength 
of steel ' 
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proof strength of HYSD bars could vary between 5 to 12 
percent in actual construction. 

4*3 Equations for Determination of Ultimate Strength of 
RB3 Section 

The equations for calculating moment capacity are 
derived from the following assumptions: 

(i) Plane sections remain plane after bending, 

(ii) Tensile strength of masonry is ignored and tensile 
force is resisted by reinforcing steel alone, 

(iii) A rectangular stress block is assumed. The width and 

depth of stress block are k^f^^ ^2^ where 

f^Y ~ strength of masonry 
a = depth of neutral axis 
k^jkg = constants. 

(iv) Maximum strain in brickwork is limited to 

0.003 and maximum strain in steel at failua^ is 
given by 

f^ 

= 0.002 + ^ (4.1) 

1.15 Eg 

where fy = yield or proof strength of steel 
Eg = Young's modulus of steel. 

Details of stress block are shown in Eig, 4.3» When 
the failure is governed by the yielding of tensile reinforce- 
ment, the section is said to be under— reinforced. Moment 
capacity for an under-reinforced section is given by 


(a) Section 


(b) Strain 


Cc) Stress block 


Fig.4'3 Details of stress block 


350 



175 


200 


3nos.8 ^ 


65 30, 65 ,30, 65 ,30, 65 , 2. 

4* ' — H ' **H •f*4* (15O-0miT> ) 


All dimensions arc ir^ mm 
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= Ag^ fy od (4.2) 

wiiere = area of tensile steel 

fy = yield or proof strength of steel 

a 

3 d = lever arm = (d - “ ) 

The depth of stress block is calculated by equating total 
compressive force 0 and total tensile force I, 


C = k^f^ ab = 1 = Ag^fy (4.3) 

which results in 

a = ~ ^ — iL- (4.4) 

and 3 = (1 > 0.5 — • — ) (^-5) 

!Ihe section is said to be over-reinforced when the failure 
is governed by crushing of masonry, !2he strain in masonry 
reaches its limiting value whereas strain in steel is 
less than its limiting value. At balanced failure, i. e. , 
both the strains in masonry and steel reach their limiting 
values simultaneously and depth of neutral axis at balanced 
failure a^ can be calculated as 




where 

690 

]£ = — 

^ 1150 + fy 

Thus, moment capacity of an over-reinforced section is 
given by 

My “ (kp3^ ^^),b(d — ■■■ ) (4»8) 

Substituting from Eq, 4.7 to the above equation 

bd% (4.9) 

The condition that a EBB section is balanced or over- 
reinforced is given by 


or 


or 


a > 

^2^1im 


tsl 

% 

> 

kgkjd 

b 

fcjfw 

^t 

% 

> 


bd” 

f^ 

Ijkg 

and k^ 

depend on 


(4.10) 


stress block and strain limitauions. 


From Y«hitney*s theory applied to reinforced concrete, 

= 0,85 and k 2 = 0.85. For k^=0.85, k2=0.85 and fy=415 N/inm^, 
Eqs. 4.2, 4.5, 4.7, 4.9 and 4.10 can be expressed in a compact 
form 


My = Ag^f^dd - 0.59 


•^st 

bd 


w 


if 


£L 

bd 


0.259 f 


w 


-^st 


< 0,319 


“W 


if 


>0.319 (4.11) 


Dimensions and area of tensile steel of a RB section Vi^iich 
will be used for reliability analysis are shown in Fig*4»4« 

4*4 P 2 x>bability of Failure of ?33 Section 

Moment capacity of a RBB section Mj, is a function 
of strengths of masonry and steel, aiid geometric properties 
of the section and can be expressed as 

Mj, = g(f,.,,fy,b,d,Agt) (4.12) 

T?Aiere 

fy^ = strength of masonry at 23 days 

fy = strength of steel 

b = breadth of the section 

d = effective depth of the section 

Ag-j; = area of tensile steel. 

'Since f^^ and fy are random variables, 11 ^. becomes random 
variable and the section can fail either by yielding of 
tensile reinforcement or by crushing of masonry. 

let P be the 'Failure event' denoted by R-S < 0 or 
R/S < 1 where R and S represent generalized resistance and 
load respectively. Due to random variations, any of the 
following events can occur: 

U - the section is under-reinforced 
0 - the section is over-reinforced. 

Thus failure may occur under any of the following events: 
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- the section fails as uikier-reinforced 
i,e., (Fnu) 

Wq - the section fails as 0Te3>-reinfO2?ced 
i.e.; (?no) 

Since the events % and Fq are mutually exclusive (this 
assumption is justified because a section failing as under- 
reinforced cannot fail as over-reinforced at the same time), 
the probability of failure p^ can be expressed as 

p^ P(R-S<0) or P(B/S<1) =P(?u)+P(Fq) C4«t3) 

Defining p^ and p^Q as follows, from Baye*s Pheoron 

= p(?u) p(Fnu) =» p(f|d). P(U) (4.14) 

Pfo*= P(Fq) = P(Fno) = p(p|.0). P(0) (4.15) 

where P(f|.‘TJ) is the conditional probability of ? given U 
and P(Fio) is the conditional pirobability of F given 0. 

Thus, Eq. 4.13 can be written as 

^^f “ Pfu Pfo 

= P (F^u) + P(-Fh 0) 

= P(fJu).P(U) + P(F.l0).P(0) (4.16) 

Probability of failure p^ can be computed as the sum of 
(i) the probability of failure of the section as an under^ 
reinforced (Pfp^) stnd (ii) the probability of failure of 
the section as an over— reinforced (p^q). 
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4.5 Reliability Analysis Foiaulation for Probalistic 
Variations of Strengths of Materials 

4.5.1 General 

Bae ultimate moment - ' capacity ox’ a RB3 section 
may be expressed in a general form as 

and = K, hi% . -If ^ X 

v/here f^^^ and Ag^ are defined as in Eq. 4.12. 

K^,K 2 ,K^ and are positive constants which, depend on 
the assiimptions regarding the stress block and strain limi- 
tations. 

!rhe strengths of masonry and steel are treated as 
random variables whereas rest of the parameters are treated 
as deterministic constants in the following derivation, 

Erom Bqs. 4.17 and 4.18, the events are identified as 
follows : 

Event U : Under-reinforced Case 

Ihe condition, that a section is imder-reinforced 
is given by 

■^t ^y 

bF ‘ fJJ 


(4.19) 





and the corresponding moment capacity can be computed 
by Sq. 4.17. 

Event 0 : Over-reinforced case 


The condition that a section is over-reinforced is 


given by 


■^st fy 

a y K 

bd f, ~ 


( 4 . 20 ) 


-w 


and the corresponding nK>ment capacity can be computed 
by Eq. 4.18. 


In both the cases, it can be noticed that the 

'■w 


resistance is a function of random variables f^ and fy. 


and thus becomes a random variable. 

4.5.2 Computation of p^ for deterministic external moment 

Probability of failure p^, for deterministic external 


moment Hg is given by 


( 4 . 21 ) 


Pf = I’CMj.iMg) 

Eor convenience, Eqs, 4.17 and 4.18 is written in the 
following form 

= C-|fy(l ■■ °2 ^ "f^ ^ (4*22) 




where 

bd 

0 — K 


f °o 

fv 

if -f - 1 


(4.23) 


(4.24) 
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“ ^1 “^st ^ 

(4.25) 

■Agt 

(4.26) 

°2 " ^2 bT 

0 , = E, bd^ 

3 3 

(4.27) 


Following assumptions are made for the derivation of : 

(i) f and f are independent positive ramiom variables 

%/ 

ranging from zero to infinity, 

(ii) positive constants, 

(iii) value of c^ is such that M^, is always positive 
quantity. 

iLll the above assumptions are justified since strength 
of masonry and steel are positive random variables and all 
the constants Cq»c^,C 2 and c^ are functions of positive valued 
parameters. Value of Cq is determined from vdiich guarantees 
that llj, will always be positive. 

From Sqs. 4.22 and 4.23, p^ given by Sq. 4.21 is 
expressed as 



= Sfo 


(4.28) 
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It sbotild be noted that the first part of the above equation 
is Pj^ and the second term is p^q as can be seen from Sq»4«16* 
Ihe event {c^fy(l - ^ 2 ^ ^ 4-«28 can be written 

as 

f^ (c^fy - Mg) < c^C2fy (4.29) 

2 

The ri^t hand part c^C 2 fy of the above inequality is a 

positive quantity. Since (c^fy-Mg) can be negative or 

positive depending on fy smA Mg, the inequality given by 

Sq, 4.29 resolves into two conditional inequalities. 

c c f ^ 

The conditionCc^ f^-M ) > 0 implies f < ^ and 

(c^fy-Mg) < 0 implies f ^ > 0 which is trivially satisfied. 


‘Ihus, 




^y 

°1^y"^e 


if f ^ > — 

y c 


V > 0 


if f,. < 


1 

Me 


y 0 


1 


(Therefore, first part of Eq. 4«28 becomes 

c Cof^ M 


-ry^e 


. P(f > - 

J c 


) 


1 


Mg Mg 

+ P({ f^>0>n {fy<Cgf^ }l {fy< ~ }) . P(fy < — 

1 y e * 




(4.30) 
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^ 1 Q 

For » the event 

c^fy-Me Cq 

It < h!£L,n{f*> 

^ W- ^ Co 

vanishes • Hius the above event can be written as 
f. 


,2 


I C -1 C n^y 

{ -2 < f„ < laax ( ^ , --f— ) } 




Therefore, Sq. 4. 30 becoiaes 


_ p ( £ -Z < f < Qax ( , — 1_LJ — ) < f < °=>) ) 

Cq w - Co Cl y 


f.. 


+ P({-^ < < «}n{0 < f ^ < 


M, 


->) 


(4.51) 


The region of integration in f,,^ and fy plane 
is shov/n in Fig. 4.5. Sq, 4.31 can be expressed in 
integral form as 
CO max' 

u"^ ^ ^ ^fw^“w^ ‘^Sv] ‘^y 

iz 


. (£z , iildz- ) 


"e 

Cl 

Me 


i tXJ 




(4.52) 


where fp^i^Cf^) and ffy(fy) axe the probability density 
functions of independent random variables f^y and fy 
respectively. 
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Second part of Eg., 4,28 can be vnritten as 

}1 


je 

% 

'3 


5fo= ^ ( tfvi :r 


c, ' ' -w ” c. 


= P( 0 < f < min ) 

” Cq Cj 


(4.33) 


iOae region of integration is shown in Pig. 4.5. Eg, 4.33 

can be expressed in integral form as 
„ mln(h,|2_) 

PfO- /[/ ° ^ dfj ff;-(fy)dfy (4.34) 


where ff.v^Cf^) and ffy(fy) are the probability density 
functions of independent random Tariables f^^ and fj 
respectively. 

It shotild be noted that the expressions derived 
for P;^ given by Sq. 4.28 represents the probability 
distribution function % 3 ,(Mg) of random variable and 
thus can be expressed as 

lilr (*c> = ^ (Mj, < 

“ Pfu + Pfo (4.35) 

v^ere p^ and p^^Q has to be calculated from Eg, 4.32 and 
Eg, 4*34. 

4.5.3 Confutation of p^ for probabilistic external moment 

Moment capacity Mj, is a random variable because 
the strengths of masonry and steel are random variables. 
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loads coming on a structure is also a random variable. 

If load is taken as probabilistic, the ex-femal moment Mg 
becomes random variable. It will be assumed that load 
and resistance are independent, i.e, , Mg and are independent 
random variables. The probability of failure of a RBB 
section for probabilistic external moment Mg can be 
expressed as 

p^ = P (M^ < Mg) 

= ^/ %g(Mg) dM^: (4.36) 

where M^, and Mg are independent random variables with 
probability density functions % (M,,) and f^c (m^) 

respectively. Assuming external moment Mg as positive 
random variable ranging from 0 to 3q., 4.36 becomes 

oo 

= / %r(Mg) %e(Mg) diMe (4-37) 

where % 3 ,(Mg) has to be calculated by Eq 4.35. 

4.6 Computation of p^ for Deterministic External Moment 

Then Strengths of Materials Pollow Normal Distribution 

Expressions for probability of failure given by 
Sq.. 4.32 and Eq,4.34 were derived under the assumption that 
f^ and fj are independent positive random variables ranging 
from zeiro to infinity. Although strength of masonry f^ 
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and strength of steel fy cannot be negative, normal 
distribution was fotind to be suitable model for f^ and fy 
as discussed earlier. !Ehe probability of a random variable 
X,following normal distribution y/ith coefficient of 


variation 6, to be negative is given by 


o 1 x-Xt, 2 

P(X < 0) = / — — -exp[-*^ (-_) ] to 
3 ^( 2 %) *-2 8 -* 


= 0 C - 1/6) 


(4.38) 


For 6 as hirh as 0.25, the probability of the random 
variable being negative is 5.71x10“^ which is very small 
for practical piirposes. Thus, normal probability density 
functions can be used in Eq,. 4*32 and Sq, 4*34 for computa- 
tion of p^ v/hen the coefficient of variation is limited to 
25 percent without introducing much error. 


Substituting normal density functions for f™ and 


fyj Eq. 4.32 becomes 

“ maxA , ) 1 1 f -f 2 

Pfa = J" [ / °° “1V ^ ^ ""Pt- 5 Jlfw] 


Sy^ 


C-^ Co 


s^ 


( . 1 ( 

2 ^ 


1 CO i If —f ^ P 

■f / [ / exp{- - (-~^) ] 

o - . tSi 2 Svv 


^ f — f 9 

1 / -"v -mx^ 


.exp (- — — ' 

ySt Sy 2 3y 


) dfy (4.39) 
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Tidiere 


“ losan and standard deviation of 
respectively 

fyjjj,Sy = mean and standard deviation of fy 
respectively, 

Eq., 4*39 can be expressed in a simplified foim as 


1 


2 ^. 22 ) - iKz^)] . “ exp(- ^ (- -y 




df^ 


'1 M, 


+ / ^ [1- ^zj] (“ “ df^ (4.40) 


v/here 


z. 




f /o — f 
- ^o -^ma 


1 f p 

2 Sy 


( 


C^C 


2^y 


Cify- Mg 




^(.)= standardized noimal distribution function. 
Similarly, substituting normal density functions for f^ 
and fy, Eq. 4.34 becomes 


M. 


mm 


iy mg 

, — ) 


Pfo = /[ / ' °5 ' -=1 exp(- ^2^ )^) df^] 

o o T2^ s^„ ^ Syr 


1 


■'w 

exp ( - “ ) df^ ( 4 . 41 ) 


YSi" Sy * ^ ^ 

where f^, s^, fyj^ and Sy are defined in Eq. 4.39* 
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Hie above equation can be expressed in a simplified form as 

oo 

Pfo = (minCz^jZ^)) - 0 ] 


1 4 f -f^ O 

exp ( - I 


l2%s. 




(4.42) 


where = (Mg/c^-f^Vs^ 

and other variables are as defined in Eq, 4-40. 


Probability of failure pf thus can be obtained by 
sunEiing the probabilistics p^ and Pfo* Probability 
distribution function is given by Eq. 4.35 where p^ 

and p^g are to be calculated from Sqs. 4.40 and 4.42 
respectively T.hen both f^ and fy follow normal 
distribution. 


BBB section shown in Pig. 4.4 is selected to study 
the variation of probabilitj’’ of failure with coefficients 
of variation of strength of masonry steel (S^y). 

Following are the details of the RBB section: 

Q 

b = 350 mm , d = 175 mm , Ag^ = 150.8 mm (3 # 8 ^), 
Strengths of masonry and steel are taken as 

f^ normal f^ = 8.96 H/mm^ 

fy normal fy^ = 449.15 n/mm^ 

Using the mean values of f^ and fy, the mean moment 
capaci% is computed from Eq, 4.11, 

= 10.989 km . 


J 


160 


-03? 6^ - 0 . 05 , probability of failure versus for 
different M^/M^ is shovm in Pig. 4.6(a). Probability of 
failure increases with increase in 6^^ and decreases 

with increase in ^o.tio. When is increased to 

0.1142, increases for same f^ and M^/Hq and is shov.Ti 
in Pig. 4.6(b). The increase in is found to be dominant 
upto 5^=0. 14 , and after that increase in p^ is found to 
be marginal. 

Effect of load factor Pfu* ^fo Pf 

different combinations of 6^ and 6^^ is shown in Pi^,4.7 
and 4.8. It can be noticed that in some region Pj^ is 
dominated by under-reinforced probability p£>^ vdiereas in 
other region p^ is dominated by over-reinforced probability 
Pfo» 131 all the cases, a region is observed where p^ is 
dominated by both p^^ and Pfo* region varies for 

different combinations of &fy and 6^^ as shown in Pigs.' 4.7 
and 4.8. In this region, p^^^ and p^^ are of the same 
order whereas in other region is governed by either 
p^ or p^Q depending on their order of m^^gnitudes. It is 
found that p^^ and p^^ are sensitive to the relative 
magnitude of a^id ratio, other variables 

remaining same. 

The balanced area of steel for the section 

shown in Pig. 4.4 is calculated fi^m Eq.. 4.11 using the 
mean values of f^ and fj and is given by 




Probability 



Fig.4-7 Effect of load factor on probability of 
failure 


Probabilit 
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•^tb = 38S.16 mm and — — = 0,387. 

Vb 

Mean moment capacity is equal to 10,939 klin for 

2 

Agt “ 150,8 mm • Effect of increase in steel area on 
Pfo Pf external moment Mg = 10,989 klfta is shown 
in Eig. 4.9(a). If the area of tensile steel is incre^ed, 
the probability of under-reinforced failure p^ decreases 
and probability of over— reinforced failure p^g increases 
but the probability of failure of the section p^ decreases 
for the same external moment. Fig. 4.9(b) shows the 
effect of increase in steel area for Mg = 7,327 kfe 

= 1.3). Thus, if area of steel actually provided 
is more than that required to achieve a particular 
reliability of the section, the design will be on safer side 
but the probability of brittle failure will increase. 

It can be noticed from Figs. 4.9(a) and 4, 9(b) that if area 
of steel is increased to a very large value, the probability 
of failure becomes constant and is equal to the over- 
reinforced probability p^g in the limiting case. The region, 
where the probability of over-reir^orced failure vdll govern 
the overall probability of failure p^, depends on the 
relative magnitude of 6^, 6^^^, ratio and 


ratio 
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4.7 Gomputatioii of for Probabilistic External Mtraent 
(Iognoraial)w4i«a Strengths of Materials Pollow Iformal 
Distribution 


Statistical analysis of typical office building 
floor loads was carried out by Ranganathan and DayaratnazaCHO). 
The frequency distribution of floor load was found to follow 
lognormal distribution at 5 percent significance level. 

External moment Mg acting on floor slabs or beams thus 
becomes lognoimally distributed random variable . Assuming 
lognoimal distribution for Mgf Eq, 4.37 becomes 

1 


= f 




In "“e 


- exp [ - - (— ) ] Wq 

2 «5ln 


(4-43) 


where 


and = parameters of lognormal distribution 

= probability distribt.tion function of Mp 
= Pfu + Vfo 

where p^ and p^g are given by Eqs. 4-40 and 4.42 respectively. 
Parameters of the lognormal distribution can be calculated as 

®ln = ^ (4.44) 

^3n ~ ^em • (-0.5 (4.45) 

where and 6jj are the mean and coefficient of 
variation of external moment Mg. 
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Substituting 

1 

VT 




w 


O' In 


BQ.. 4.43 can be v/ritten in a simplified fonn as 


ao *1 


Pf = / 


•^if/ 

.e .F: 


<75 


Mr 


(Ml,.,, e 


^in ) dw 


(4.46) 


vdaere and are given by Eqs. 4.44 and 4.45 respec- 

tively. dThe integral given above is in standard form and 
is evaluated using G-auss- Heraite quadrature formula (ill). 

Probability of failure for probabistic variation 
of external moment is given in Table 4,2. External moment 
Mg is taken as lognormal ly distributed with, mean Mgjjj and 
coefficient of variation Mean value of is calculated 

by Sq. 4.11 using mean values of f^ and fy. Table 4.2 shov/s 
the variation of p^ for different and ratio, 

last column of Table 4.2 gives probality of f 3111106 for 
deterministic external moment, i.c., Failure 

probability increases as the variability of external 
moment increases and decreases as ^■^rm^^'em iiTcreases. 

It should be noticed that the ratio is the central 

safety factor and also called load factor applied to the 
mean values in a design problem. 



Table 4.2 ; Probability of Pailure for Probabilistic External Moment 



« 0,1142 and * 10.989 Klta 
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4.8 Monte Carlo Simulation 


Monte Carlo simulation technique (73) is used to 
generate random samples of moment capacity to study the 
probability'- behaviour and to estimate different parameters 
of iij,. Details of the simulation procedure is given in 
Chapter 2. Choice of sample size in simulation plays an 
important role. Generated samples of is used to estimate 
the mean and standard deviation. Iiarger the sample size 
used, estimates of mean, standard deviation etc. vd.ll be 
closer to their population values. Hie minimum sample 
size required to estimate a particular parameter depends 
on the desired accuracy of the estimate. The minimum sample 
size n, for the estimate of population mean with a' 
confidence level (l-a) percent, is given by (112) 



(4.47) 


where 


= sample standard deviation with n-1 as 
dencHEinator 

t = accepted error in the estimate of mean 


i ^/2 = - a/2) . 

For large samples, S 2 .(with denominator n-1 ) gives a good 
estimate of the population standard deviation. For an 
acceptable error of + 5 percent of the standard deviation 
and confidence level of 95 % , the minim-um sample size 

required is n == 1537 and for confidence level of 99 %, 
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n - 2663 * In the present work, 3OOOO samples are generated 
to study the parobability behaviour of lilj.. 

Samples of masoniy strength f^ and steel strength 
fy are generated independently from their parent distributions 
(normal in both cases). First fifty generated samples of f^ 
and fy following M (8.96, 1 , 26 ) N/nm^ and N(449.15 , 51.29)N/iBm^ 
are shown in Figs. 4 . 10 (a) and 4, 10(b) respectively# From 
each of the generated samples of f^ and fy, sample of is 
generated by Sq. 4,11 for the EBB section shown in Fig, 4,4. 
First fifty samples of using the generated sscaples of f^ 
and fy are shown in Fig,4, 10 (c), Since the section is highly 
under-reinforced deterministically (Ag^/Ag-|j]3=0,587), the 
samples of are dominated by the variation of steel strength 
fy and not by the masonry strength f^ as seen in Fig, 4 # 10 (a), 
(b) and (c). As the steel area is increased towards balanced 
steel area, dominance of both f^ and fy on is observed. 

The coefficient of variation of is found to be dependent 
on the relative magnitudes of coefficients of variation of 
f^ and fy. 

©wo typical histograms of the simulated samples of 
of the section given in Fig, 4.4 are shown in 
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26 3! 36 41 

p{« number 

es of fv, N (8-96, 1-26 
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Samplt nymbtr 

Cc) Samples of using genented sompks of f 
and fy for Ast =150*8 mrr? 


rig 4*10 r^andom sample 






Figs. 4 . 11 (a) and 4.11(b), Hie coefficient of variation of M^, 
is found to be around 11 percent in both, cases. For 
the saae C,0,¥, of steel strength, (^fj - 0.1142), the 
shape of the histogram changes vd.th the coefficient of 
variation of masonry strength Usual ohi-squars test 

is conducted to fit the data. Ihe samples shown in Fig. 

4 . 11 (a) is found to follow normal distribution at 0,5 percent 
significance level. Is 6^^ is increased to 0.20, the shape 
of the histogram changes drastically as seen in Fig,4« 11(b) 
and dses not follow normal distribution at 0,5 percent 
significance level. Fig. 4.12 shows another typical histo- 
gram of simulated samples of Mj. for 6 ^ = 0.05 and 6 ^^ y = 0 . 20 . 

It can be seen that the histogram is negatively skewed vath 
a coefficient of variation as low as 5.18 percent. Normal, 
lognormal, beta, tsrpe I extremal (anallest) and Type III 
oxtremal (smallest) have not satisfied chi-square test for 
the generated data even at 0,1 percent significance level. 
Samples of are generated using different combinations 
of and 6^ to study the shape of the frequency 
distribution. It is observed that for a f ix*' 1 &fyf the 
distribution becomes more and more negati/ely skewed 
with increase in For highl 3 r under-reinforced sections, 

the coefficient of variation of steel governs the coefficient 
of variation of moment capacity and normal distribution 
can be accepted as an approximate model to represent the 
probabilily distribution if 6f and are limited to 
15 percent. 
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2 

Y/hen the steel area is increased from 150*8 Emi 
( 3 Nos, 8 to 235*62 nm^ ( 3 Nos. 10 ^), is fotiM to 

increase from 5.18 percent to 6,25 percent for same 
and 6^ = 0,05. The effect of increase in area of steel 
is shown in Fig, 4*13. The frequency distribution is found 
to be more negatively skewed and none of the common standard 

f 

distributions satisfied chi-square test. Samples of are 
generated y/ith different steel area for different combinations 
of 6^ and It is observed that as the area of steel is 

increased towards balanced steel area, the frequency 
distribution becomes more and more skev/ed to the ri^t. 

The samples of as given in Fig, 4 . 13 acco separated 
into under-reinforced and over-reinforced cases using the 
condition given in Eq. 4.11. The samples for which the 
section is under-reinforced and over-reinforced are separated 
and histograms are shovm in Figs. 4.14 and 4,15 respectively. 
It can be seen that the frequency distribution of under- 
reinforced moment capacity has neglible skewness whereas 
the frequency distribution of over— reinforced moment capacity 
has a marked negative skewness, Kie coefficient of 
variation of is found to be equal to 5.18 percent 

whereas the ©oefficient of variation of found 

to be equal to 15 percent. Samples with different combina^- 
tions of 6^, 6^ and are generated. It is found that 
the coefficient of variation of over-reinforced moment 


Relative frequency tti Rclafivc frequency 



ig.4*14 Histogram of Mfu separated from 
shown In Fig. 4*1 3 
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capacity is governed mostly by the C.O.Y, of masoniy 

strength 6^^ idiereas the coefficient of variation of under- 
reinforced moment capacity is governed mostly by the 

C.O.Y. of steel strength Sfy* "fhe number of samples in 
over-reinforced case is found to be 726 out of 30000 sample 

as shown in Pig, 4. 15 which means that the probability of 

2 

the section given in Pig. 4*4 with Ag-|; = 235.62 mm being 
over reinforced is equal to 0.0242. The probability of 
a section becoming over-reinforced is a function of 
and area of steel It can be noticed from Pigs. 4.14 and 

4.15 that the mean of under-reinforced moment capacity 
is higher thaci that of over-reinforced moment capacily 
T-Ipom while the range of is higher than that of M 2 ^. 
Similar findings were observed in the simulated samples 
of prestressed concrete beams (85,86). This again depends 
on the relative magnitudes of 6^^^, 6^ and area of steel 
provided. 

Simulated samples of shown in Pig, 4.14 did not 
follow nonaal distribution at 0.5 percent significamce level. 
Similarly, other distributions like lognormal, beta, Type I 
extremal (anallest). Type III extrcmial (smallest) etc. 
did not satisfy chi- square test to the samples of 
Simulated samples of also did laot follow any of the 
above distributions. Samples of and are generated 
for different combinations of 6;^ and 


It is observed 
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that samples of follow normal distrihution at 2,5 percent 
significance level for 6£^=0.20 and 6^=0.1142 whereas 
for 5jf^=Q,20 and 6^==0,05 didnot fit normal distribution 
even at 0,5 percent significance l^-vel. Similarly, it is 
found that the samples of over-reinforced moment capacity 
,f or 6£^^=0,20 and 6^=0, 11 42, followed Type III extremal 
(smallest) at 1 percent significance level. Shapes of the 
frequency distributions of M 3 ,^ and I^I^Qwere found to change 
with the relative choice of S^y^jb^yand Ag^. Thus particular 
distribution like normal for under-reinforced moment 
capacity and Type III extremal (smallest) for over-reinforced 
moment capacity cannot always be accepted. 

Probability of failure p^- is calculated from the 
simulated samples of Mj. in the following way. Bie number of 
samples falling below is computed and p^ is calculated as 

Pf = P ( < Mg ) s ^ (4-48) 

vdiere 

Mg = external moment (deterministic value) 

= number of samples less than or equal to Mg 

n = total number of samples generated. 

Probability of failure computed from simulated 
samples by Iq, 4,48 for different values of Mg are shown in. 
Table 4.3 and compared with the results obtaiiu^d by analytical 
method throu^ Eqs, 4.28, 4.40 and 4.42. It can be seen 



Sable 4*3 5 Oomparison of Simulated Results 



Properties of the section are as shown in Fig 
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from fable 4-3 that the values of p^. by two different 
methods have very small difference. It is worthwhile to 
mention that minimum pixjbability of failure etjaputcd 
by Bq.. 4*48 is dependent on the number of sasples generated. 
In this case, it is equal to I/ 5 OOOO ,i.e. , 3.33x10 . 

4.9 Discussions and Conclusions 

Ultimate mcment capacity of a HB3 section is a 
function of geometric proportions , area of steel and strengths 
of masonry and steel. Due to random variations in strengths 
of masonry and steel, there exists a probability that the 
section will fail as an over-reinforced even though the 
section is designed as an urder reinforced based on 
deterministic analysis. Probability of failure p^ of the 
section is the sum of the probabilities Pf>^ and Sie 

computation of failixre probabilities p^ and p^^^ involves 
evaluation of multiple integrals and depends on the 
probability distributions of masonry strength and steel 
strength. The probabilistic variations of dimensional 
properties of the section are not incorporated in the 
present formulation. The expressions for p£>^ and p^^ 
are derived from the equations of omnent capacity given 
by Eqs. 4.17 and 4.18 which are also valid foe representing 
the mcment capacity of singly reinforced RO and prestressed 
concrete sections. Thus, the expressions of p£>^ and can 
be used for reinforced and prestressed concrete beam 
sections aia>. 
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The probabilities and Pj^^ are found to be 
sensitive to relative magnitudes of the coefficients of 
variation of strengths of masonry and steel and M^/M^ 
ratio* If the area of steel in a HBB section is increased, 
the probability of the section being over-reinforced 
increases* Consequently, under-reinforced failure probability 
(P:g^) decreases and over-reinforced failure probability (pfo) 
increases but the probability of failure of the section (pj) 
decreases for the same external moment* Therefore, if the 
area of steel actually provided is more than that required 
to achieve a particular reliability of the section, the 
design will be on safer side but the probability of brittle 
failure will increase. It is observed that p^ decreases to 
a limiting value eqial to p^Q as the steel area is increased 
to a very large value* This signifies that as the 

probability of failing as an under-reinforced is almost 
equal to zero for heavily oveivreinforced sections. 

The situations vdiere the probability of over-reinforced 
failure v/ill govern the overall failure probability p^, 
depends on the relative magnitudes of coef:^icients of 
variation of strengths of masonry and steel, -fi-gt/'^tb 
and 1^/Mg ratio. For probabilistic external moment, 
the probability of failure increases as the coefficient of 
variation of external moment increases. 

Samples giving ultimate moment capacity are 

generated by Monte Carlo simulation using the randomly 
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generated data of strengths of masoniy and steel. 

Histogram of the simulated samples has shown a marked 
negative skewness. As the percentage of steel area is 
increased to that of a balanced value, the distribution 
of Mj, becomes more and more skewed to the ri|^t. The 
simulated sample has not followed any standard probability 
distribution. The samples of M^. are aepeccatad into und®e- 
reinforced and over-reinforced cases. Normal and IVpe III 
extremal (smallest) distributions have been tried to fit 
the under-reinforced and over-reinforced samples respectively. 
Shapes of the frequency distributions of under-reinforce<i 
and over-reinforced samples change v/ith the relative choice 
of coefficients of variation of strengths of masonry and steel, 
and also with the Agt^-^tb Although normal and 

T3rpe III extremal (smallest) distributions are foxmd 
satisfactory to represent the under-reinforced and over- 
reinforced samples in some particular cases, but they 
cannot be acc jpted in gener-:il. For a highly under- 
reinforced section (deterministically) ,noi3nal distribution 
is found to be acceptable to represent the nK>ment capacity Mp, 
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RKDIABIHIY BiSED DESIGN OP REI2IPOHGED BRIGK BEAMS 

5.1 General 

Reliability based design procedure of reinforced 
brick beam section at ultimate flexural strengtbs'is presented 
in this chapter. Methods of computation of probability of 
failure for deterministic and probabilistic external mcment, 
considering the probabilistic variations of strengths of 
materials, are given in Chapter 4. Design of a RBB section 
for a specified probability of failure reduces to a problem 
of solving an integral equation to find the geometrical 
properties or the area of steel. She method is illustrated 
with examples. Semi-probabilistic limit state design method 
is also illiistrated, A set of partial safety factors for 
different variabilities of load and strength is also 
presented. 

5.2 Introduction 

Design methods can be broadly classified into two 
diffeirent categories : (i) deterministic design and (ii) 
probabilistic design. 
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In the detciiainistic design .-nethed, ell the design 
Tariables are trca,ted ns deterministic quantities. Hie 
safety of the structure is ensured by providing more resistsuace 
than the external forces. Working st .'oss method and ultinu».te 
load design method are based on deterministic .-spproach. In 
working stress method, the safety is ensured by limiting the 
actual stresses to a set of permissible or allowable stresses 
whereas in ultimate load design the safety is ensured throu^ 
proper sp<|cification of load factors. 

Loeds acting on a structure arc probabilistic rather 
then deterministic. Similcnly, strength of mat '.rials such as 
concrete, masonry, steel etc. are randcan variables. All 
these design variables are trea-ted as probabilistic in 
probabilistic design. IHie safety of the structure is e:q>ressed 
through a specified probability of survival (reliability) 
that the strucimire will not be unserviceable over a given 
period or throughout the expected life of the structure. 

Hiree degress of sophistications in probe.bilisric design are 
possible considering the statistical variabilities in resistance 
and lOcd(l13)» At the highest sophistication, termed as level III, 
the probability distribution of all variables and necessary 
convolutions must be known or computed and the members are 
designed for an accepted probability of failure prt^viously 
agreed upon by the profession. Hie reliability is measured 
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or specified by a safety index or reliability index P in a 
level II design criteria #iere tlie safety is checked at 
certain discrete points (i.e,, at selected values of mean 
live load to dead load ratio) of the limit state equation(72, 
75,11 3)« Probability distributions are not required in 
level II methods. Level I methods involve selection of only 
one set of loM and resistance factors vdiich can be applied 
to all designs regardless of mean live load to dead load 
ratio (113)* Cfurrent design practices are of Level I type 
although Levels I and II can be made equivalent by allowing 
load and •resistance factors to vary with different load 
ratios. 

5,3 Reliability Based Design of RBB Section 
5.3.1 Introduction 

Ihe probability of failure pj^ of a R!® section, as 
seen in Chapter 4, is the sum of under-reinforced failure 
probability p^^ and ove 2 >-reinforced failure probability Pfo* 
Both these probabilities are functions of mean values and 
coefficients of variations of masonry strenth and strength 

4 

of reinforcaaent and their probability distributions. Apart 
from these, p^^^^ and p^^ are functions of dimensional properties 
of the section, area of steel and the equation ’,#iich defims 
the section as under-reinforced or over-reinforced. Por 
deteimiinisiic external moment Mg, probability of failure can 
be expressed as 
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Sf = Pfu-^Pfo (5.1) 

utoere = g^(Mg, f^, 5^^, b,a, ig^.) 

Pfo “ S2^^e’ ^wD’ ®fw* b,d, -A^-t) 

Probabilities Pj^ and are also dependent on the type 
of probability distributions of strengths of the masonry 
and steel, Por probabilistic external moment, p^^i ®nd p^^j 
depend also on the type of probability distribution function 
of external moment (Mg) and its parameters. 

In a given design problem, the risk level p^ is only 
specified. Given the distributions along vd.th necessary 
parameters (like mean, standard deviation etc.) of Mg, f^ 
and fy, area of steel can be solved from Eq., 5,1 for an 
assumed set of values of b and d, Hius any set of b,d and 
Ag.jj, vdiich gives probability of failure less than the 
prescribed risk, is a feasible design. Breadth (b) ard 
effective depth (d) of the section can bo chosen frc© 
practical limitations and Ag^ can be solved for. As seen in 
Chapter 4, Eq, 5.1 is an integral equation, Design procedure 
for deterministic and probabilistic external moment is given 
in the following sections , 
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5*3*2 Reliability based design for probabilistic variations 
of strengths of materials and deterministic external 
moment 

It is assumed here that the dimension of the section 
is given and the area of steel is to be determined for a given 
probability of failure* If the external bending moment is 
assumed as deterministic and strengths of materials follow 
normal distribution, area of steel can be determined by solving 
the integral equation given by Eq* 4*28 , 4*40 and 4*42 for 
assumed values of b and d* The method is illustrated through 
the following example. 

Example 5*1 

A reinforced brick beam section ha,s to carry an 

external bending moment of 7 ktTm, Desi^ the section for 

-5 

probability of failure of 10 • Strength of masonry f^ and 

p 

strength of steel f„ are distributed as F (8,96,1,26) E/mm 

if o 

and E (449. 15> 51.29) E/mm respectively. 

G-ivenj f^ = 8,96 E/mm^ , E/mm^ 

fyjjj = 449.15 E/mm^, Sy =51.29 E/mm^ 

Mg =7 kEm (deterministic value} 

Pf = 10-5 

Breadth and effective depth of the section are assumed ast 
b = 350 mm , d = 175 mm. 
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Substituting above values in Eqs. 4.28, 4.40 and 4.42 , area 

of steel Ag^ is found by iterations. 33ie iteration is 
started with an initial value of Ag-^ and the value of is 
computed. Depending on Tfitoether the computed p^ is sm^ler 
or greater than the given probability of failure (i.e,, 
p^ =? 10 ^), Ag^ is decreased or increased correspondingly 
in the next iteration. Iteration is stopped vdien the 
difference between the computed p^ and given p^ is less 
than a specified tolerance. Using this procedure, Ag^. is 

p 

found to be 182,47 mm , Mean moment capacity for 
Ag^ = 182.47 mm^,. is equal to 13,076 kltn and load factor 
M^/Mg is equal to 1,87. 

In this example the dimensions of the section are 
assumed and Ag^ is computed. It is possible that for an 
assumed set of b and d, the iterations may not converge 
and the computed p^ is always higher than given p^ even 
for an increasing value of Ag^. This means that the section 
cannot be designed with those assumed dimensions of the 
section and accordingly either b or d has to be increased. 

IThe coefficients of variation in the example were given as 
6fw = 0,1406 and 0,1142, 

For an assumed set of b and d, area of steel required 
for a specified probability is a function of and Mg.' 

Table 5.1 shows the effect of increase of 5^ and Mg on the 
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P-5 

Table 5h1 : Area of Steel Required for ~ 10 , for 

Deterministic External Moment 


^e 

6^^=0,1406,6 

fy=0.05 

^f W~^ • ^ 4-*^^ » ^fy“^ • 

1142 

^st 

(xnm^) 

^rm 

iT” 

A* 

std 

(mm^) 

^st 

(mm^) 

Me . 

(mm^) 

7 IcRm 

119.49 

1.264 

128.03 

182.47 

1 .868 

144.68 

8 kRiQ' 

157.80 

1.264 

150.14 

211.12 

1.866 

169.66 

9 kRm 

156.83 

1.266 

173*75 

243.60 

1.871 

196.34 


steel by limit state of strength (see Appendix-B) 

Rote; b = 550 mm , d = 175 nmi 

f^ « 8.96 N/mm^ , = 449.15 R/mm^ 
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required Ag-(;» It can be seen from Table 5.1, that the steel 
area increases with increase in and 6^ for same 6^^ as it 
is expected. The area of reinforconent by limit state of 
strength is also given in Table 5.1. 

5,3,3 Reliability based design for probabilistic variations 
of strengths of materials and external moment 

External bending moment Mg is assumed as lognoimally 
distributed random variable, For f^ and f^ followir^ normal 
distribution, and Mg following lognormal distribution, area 
of steel can be found in a similar way as discussed in 
section 5.3*2 by solving the integral equation given by 
Els. 4.46 , 4.28, 4.40 and 4.42 for assimied values of b and d 
of the section. The method is illustrated through the following 
example. 

Example 5,i2 

^ -5 

Design a RBB section for p^ = 10 when external 
moment Mg is lognoimally distributed with mean 7 kNm and 
coefficient of variation 20 percent,- Other parameters are 
same as in Example 5.1. 

Given: 

Mgjjj = 7 kMa , % = 0.20 

^wm "" » s^= 1.26 R/ram^ 

fy^ = 449.15 R/mm^,Sy= 51.29 R/mm^ 

P£ =* 10 ^ . 
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Breadth and effective depth of the section are assmed as 
before 

b = 350 mm , d = 175 mm 

Mg is lognozmally distributed and the parameters of log- 
normal distribution are computed from Eq.s* 4-. 44 and 4*45 
and they are: 

01^ = 0.193 

%n ~ 6.864 Ma 

Ihe solution of the integral equation is obtained by 

iteration as discussed earlier. Por an assumed value of 

Ag^, P£ is computed by Eq.* 4*46 where Pjj^(x) is given by 

Eq.s* 4.28, 4.40 and 4.42. Comparing the computed p^ with 

given p^, Ag^ is changed accordingly in the next iteration 

till the convergence is achieved. Ag-t is found to be 
2 

275 mm and the corresponding mean moment capacity is 18,746 
kHn, Central safety factor or load factor ) is 

found to be 2,61S, Table 5.2 shows the effect of increasing 
and 6^ on the required steel area for same i't can 

be observed that for the same section, required ^st increases 
as the coefficient of variation of external moment increases 
which is expected, load factor is more in case of 

probabilistic Mg than that of deterministic Mg and increases 
with coefficient of variation of external moment as seen in 
Tables 5.1 and 5.2, This is expected because as the statistical 
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Table 5*2 : Irea of Steel Required for = •jo“^ , for 
Probabilistic External Moment 




^fw “ 

0.1406 




= 0.20 


0.25 


■^st 


■^st 

^'^rm 


(mm^) 

Mp>Tn 

(mm2) 

Mem 

0.0500 

254.0 

2,331 

297.0 

2.857 

0*1142 

275.0 

2.678 

347.0 

3.244 


Rote: External moment is lognormally distributed with 

mean = 7 klfei, 

ill other values are same as given in Table 5*1* 
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variation in exteomal load increases, higjier safety factor is 
needed for the sane reliability. 

External moment Mg acting on a section is equal to 
the sum of dead load moment Mjj and live load moment Mj^. Live 
load variation is significant compared to dead load variation. 
Eor small variation of dead load, dead load moment can be 
assumed as deterministic. Under this condition, external 
moment Mg can be written as 

“e = % + % (5.2) 

where is deterministic value (i.e., constant) and 

is probabilistic following probability densiiy function 
Assuming lognormal distribution for the 
probability density function of Mg can be written as 

1 .e 1 ^Iln. 

< Mg < ~ (5.3) 


where and are the parameters of lognormal 

distribution for live load moment Mj^* Probability of 
failure p^ can be expressed using Eq, 4.36 and is given by 


Pf = ^ %e(Me) <3Me 

IVinrv 


= fXK^- 


M. 


.dMo 




cJIln 


( 5 . 4 ) 
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Substituting 


1 l2l((M.-%)/ Mjij, ) 



Eq.. 5.4 can be written in a simplified foim as 

Pf %n n )* <irf (5.5) 

E<1. 5.5 gives the probability of failure ^en external moment 
consists of deterministic dead load moment and lognormally 
distributed live load moment. Por normal distributions of 
f^ and fyf ’to 5e calculated by using Eqs. 4.28, 

4.40 and 4.42. Parameters of lognormal distribution M jn-n 
and can be found from the mean and coefficient of 

variation of live load moment as usual. In an analysis 
problem, the parameters in the ri^t handside of Sq, 5.5 are 
known and the integral can be computed by using Gauss-Hermit e 
quadrature formula. On the otherhand, in a design problem, 
area of steel can be determined by solving the integral 
equation given by Eq, 5.5 for an assumed set of b and d. The 
solution procedure is same as given in example 5.2. 

5,4 Semi-probabilistic limit State Design of RBB Section 

5,4.1 Introduction 

A reinforced brick beam section can be designed for 
a specified probability of failure if ■ 
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uhe pr’Oba.bili’GAT dis'crlbii'';ions ot s'fcr'sn'i'ths of 

materials and load along with their parameters are known. 

This involves complicated procedure of solving integral 
equation which is not suitable for routine design in 
profession. Moreover, complete knowledge of the probability 
distributions of different variables are not known in real 
life* To avoid complicated design procedures, probability 
based deterministic code fomats were suggested by many 
researchers (70,71 ,72, 74j75)* limit state design methods 
can partly incorporate statistical variations in some of the 
variables in a deterministic w^’" and may be called semi- 
probabilistic approaches. Probabilistic variations Of 
strengthsof materials and loads are taken into account 
through characteristic strengths and loads, and the safety 
is ensured by proper prescription of partial safety factors 
applied to strengths and loads for different limit states 
considered. Codes based on limit state format for reinforced 
concrete are being used in the profession ( 39 , 40,41 )• Code 
of practice based on limit state format for reinforced brick- 
work has come in a draft fom (48) incorporating statistical 
variations in strengths of masonry and steel, and loads. 

Based on the formats suggested by others ( 70 , 71,72 , 74 , 75 , 113 , 
114 ), partial safety factors are arrived at using the results 
of statistical analysis. 
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5«^,2 Ciiaract eristic strengths of materials 


Characteristic strength, of a material is defined as 
the strength, below v^hich the test results are expected to fall 
with' an agreed probability. It is given 'oy 


% - 


where 


% - 


(5,6) 

characteristic strength of material 
Ejjj = mean strength 
s = standard deviation 

fc = a factor which defines the acceptable risk, 

BS 5628: Part t, 1978(45) defines characteristic strength of 
masonry as the 'strength below which the probability of test results 
•fculling is not more than 5 percent Assiming normal 

distribution ?/ith accepted risk of 5 percent, value of k 
works out to be 1,645(111), Thus characteristic strength 
of masonry f^^i^ is given as 


^wk ^wm s^ 

= (1 - 1.645 6j„) (5.7) 

where f^j and 6^^ are mean, standard deviation and 
coefficient of variation of strength of masonry respectively. 
The GBB-PIP Committee (82) recommends 5 percent accepted 
probability of test results falling below the characteristic 
value for concrete and steel. Thus characteristic yield or 
proof strength of steel assuming normaJ. distribution is 
given by 
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V “ -3® “ Sy 

= f^d - 1.645 6^) (5.8) 

7/here , Sy said 5^y are mean, standard deviation and 
coefficient of variation of strength of steel respectively, 

5 • 4 . 3 Che r act eristic 1 o ads 

Characteristic load Sjj. is defined as the load which 
has a 95 percent probability of not being exceeded during the 
life of the structure (41 >43). Since data is not available 
to express loads in statistical terms, the loads specified 
by I.S. Code (115) can be assumed as the characteristic load (41). 
Characteristic load may be expressed as 

Sk = ( 1 +k 6) (5.9) 

where = value of the most unfavourable load with a 

50 percent probability of its being exceeded 
once in the expected life of the structure 
6 = coefficient of variation of load 

1c = coefficient depending on agreed probabxli'uy 

of loading greater than Sj^. 

Dead load on a structure may be assumed as norau-tlly distributed 
with small coefficient of variation for practical purpose, 

Dor an agreed probability of 5 percent that the dead load 
will be greater than characteristic dead load Di^, 3q,5.9 becomes 
% = D^ (1+1.645 6jj) (5.10) 
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where = mean value of dead load 

~ coefficient of variation of dead loa,d. 

Ihe frequency distribution of floor live load W3.s found to 

follow lognormal distribution from a load survey conducted 

at I* I, 'I, Kanpur(1l0)* Assuming live load follows lognormal 

distribution, characteristic live load can be expressed as 

ln(Li./lii]a ) 

p (1 < Ljj.) = 0 -) = 0.95 (5.11) 

^In 

where I n ri and are the parameters of lognormal distribution. 

Using the relations between parameters of lognormal distribution 
vlth mean and coefficient of variation, Iq, 5.11 can be 


written as 


exp [ 1.645 f(ln(l+6^)) ] 

V 1 +6][j 


( 5 . 12 ) 


where = mean value of live load 

= coefficient of variation of live load. 

Eq. 5.12 gives the relation between the characteristic live 
load vvdth the mean live loal. 


5,4,4 Partial safety factors 

Probability of failure of a RB3 section can be stated as 

Pf = ^(Mj, < Mg) = P(F^ - Mg < 0) (5.13) 

where and Mg are moment capacity and external moment 
respectively. If the probability distribution of the safety 
margin M^-Mg is known, p^ can be easily calculated. Ihe 
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distribution of is important in the calculation of 

probability of failure end not the distribution of other 
Yariables which affect the safety margin 1/I^-Mg(70). Structural 
design is generally done for low probability of failure 

f- 

(i.e,, Pf < 10 " say). At very lo\7 risk level (for P£<10~^), 
the calculated p^ is sensitive to the distribution of the 
safety margin. R.eg8.rdless of the distribution of 
and for independent Ivlj, and probabi3.ity of failure can 
be vrritten as (70,72) 


Pf = - 


,rr 2" 

4 


■) 


(5.14) 


where distribution function of the normalized 


Z 


rcUidom variable X defined by 
11 ^) 

^ 


(5.15) 


Y 7 here “ ^-an of ILj, and respectively 

% , - standard deviation of and. Mg respectively. 

Sq. 5.14 can be written as 

^Pf^- 4^ (5.16) 

viiere (. ) is the inverse function of Pj(. ) • 

Assuming the distribution of the safety margin (M^^-Mg) follows 
normal distribution, the distribution function P2£(3:) becomes 
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standardised normal distribution function and Eq,,5.'16 becomes 

Mim = “em + P 4 > (5.17) 

i 0 

— .i 

where p = - 0” (p^) 

= for < 0.5. 

The safety index P relates mean values of moment 
capacity and external moment in terms of standand 
deviation of the safety margin, for normal distribution , 

p is rela,ted to the probability of failure as 


Pf 

p 

10“^ 

3.09 

lo""^ 

3.72 

10“5 

4.26 

10“^ 

4.75 


To ensure that the probability of failure v/ill be less 
than p^, Sq. 5.17 takes the foim of following design 
inequality • 

> Mem ^ (5.18) 

If the distribution of is not normal, fq, 5.16 can 




still be expressed as 3q, 5.18 vi/ith P defined as P — (p^ 
Defining a quantity as 


a 


1 


St.t + SJI 


(5.19) 
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Eg.. 5»18 takes the follovdng fom 


or hij > + P«i %, ) 


(5.20) 


'The central safety factor required for the design is given 

v = (5.21) 

^em (l~pa^ ) 

'The design inequality given by Bq. 5.20 can be e:cpressed 
in the convenient form. 


^rm E M0 jj3^ (5.22) 

where y = (1 ~ pa.j 

S = (1 + pa., 6 ^^) 

'The coefficients r and 1 are the strength reduction factor 
and load fantor applied to the respective mean values. The 
two sides of the design inequality given by Bq. 5.22 are 
coupled by a^. Talue of a.j given in Bq, 5.19 depends on the 
rela,tive values of %/r_^ and sjg;^. Bor any combinai,tion of 
values of ana g^. , Vcilue of a. lies betv/een 0,707 to 1,0. 
Bor a pra-cti cal range of s^i^ and Sjj , value of a. v/as found to 
be 0,75 by minimi sing the error in the linearization given 
by Bq. 5.19 (72 , 75 ) . Value of a that minimize s 'ihs error 
depends on the range selected for % /s7,£ ratio. Assuming 


1 
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= 0,75, th.e resistajLce and loa-d side of Iq. 5*22 becomes 
independent of eacb other but interconnected by P, Thus, y 
and ]? in Sq. 5.22 becomes 

r = (i - 0.75 P 6jj^) ( 5 , 23 ) 

P = (1 + 0.75 P %„) 

External moment Mg can be expressed as the the ston of dead 
load moment Mq and liye load moment Hj^ 

Mg = + (5.24) 

‘Hie mean vs-lue and standard deviation of Mg is given by 



M- 


Ilm 


+ M 


Im 




(5.25) 


where , M j^ = mean values of and respectively 

% ’ = standard deviations of Mjj and respectively. 

Substituting a.j = 0.75 , and using the relations given in 
Eq. 5.25, Eq, 5.20 becomes 

I.L_ > lLj)„ + + 0.75 P (3 m^+ + ) (5.26) 


•■"rm 

Esfining another qacntity ag 

Eq. 5.26 becomes 


(5.27) 


2 S’!) 


Y Mrm 


(5.28) 
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where 


r 

= (1 

- 0.75 


(5.29) 


= (1 

+ 0.75 

^“2 %) 

(5.30) 

~L 

= (1 

+ 0.75 

Pag 

(5.31) 


Bq, 5.28 gives general foim of design inequalrby Y/here Ij) 
and are the partial safety factors applied to the mean 
values of dead load and live load moments respectively and 
r is the strength reduction factor to be applied to the 
mean value of moment capacity, Partia,! safety factors 


and are coupled by which is given in Bq, 5.27. fhe 
coefficient given. by Bq, 5.27 can be written in a conveiuent 
foim as 

a 

As seen in the above equation, oi2 depends on the magnitudes 




(5.32) 


of htQ ’ ratio, 

moment to mean live load moment 


Effect of mean dead load 
ra-tio, for different values 


of and partial safety factors Bp and Bp are shown 

in Bigs. 5.1 and 5.2. Range of ratio of mean live load , 
moment to mean des-d load moment is selected as 0,5 to 2,5. 
Bigs. 5^1 and 5.2 are plotted using Bqs. 5.30 to 5.32 for 
coefficient of variation of dead load moment = 0.10. ir 
can be seen from Bigs. 5.1 and 5.2 that the safety/ factor Bp, 




Ffg.5t Variction of partial safety factors 
with ratio 




206 


required to aciiieire a particular reliability (i.e., for constant 
^), increases with. also with coefficient 

of variation of live load moment Since partial safety 

factor I'd is coupled to 1-^ by a^, a slight increase is 
observed with increase in and 3:atio. For a 

particular value of P, the change in Fjj is insensitive to 

and ratio as can be seen from Figs. 5.1 and 5.2. 

Variation of strength reduction factor y with ooeffioient of 
variation of moment capacity is shown in Fig. 5.3 
for different values of reliability index p. Strength 
reduction factor y decreases linearly as increases for 
a particular value of P« As the risk level decreases (i.e., 

P increases), strength reduction factor y decreases for 
any particular value of 6j|^ as can be seen from Fig. 5.5. 


* 


From the strength reduction factor r , the partial 
safety factor to be applied to strength of masonry is derived 
as follows. -The design moment of an under-reinforced section 


can be written as 






(5.33) 


where f^^ = characteristic streng"ch of steel 

^wk“ eristic strength of masonry 

Yj = partial safety fa,ctor applied to steel 
7 ^ = partial safety factor applied to masonry. 
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Substituting mean values of f^ an.d fy in Eq. 4.11, the 


mean moment capacity is given by 

^t 


= ^t V ^ ~ 


bd 


ym 


‘■van 


(5.54) 


Since r is the design moment, partial safety factor 


TC^ can be found by equating andrM^ 
^d “ 


(5.35) 


Substituting and given by Eqs. 5.33 and 5.34 
respectively in the above equation and using the following 
relations between characteristic and mean strengths 

(5.36) 


^yk ~ 1.645 ^fy) fyjj^ - 0-j 

^vdc ~ 1«^45 ^vzm 

Partial safety factor is given by 


(5.37) 


w 


0.59 C 


(1 - 


rr-, 


vdiere 


3 


0 


1 


(1-0.59 C, -^)) 
^ °1 


(5.38) 


'yk 


3 bd f-yylC 

Pajrtial safety factor for steel ( is established as 1.15. 

Pig. 5.4 shows the plot of ^ computed from Eq. 5.38 for 

4st ^yk 

different values of ^ 3 - and strength reduction . 

bd f^^ 

factor y • Coefficients of variation of strengths of 
masonry and steel are chosen as 15 and 5 percent respectively. 
For other, combinations of coefficients of variation of 
strengths of masonry and steel, is found to change 
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marginally and is not reported here. 


She value of reliability index p is 4-* 26 for an 

accepted probability of failure of 10 . She coefficient 

of '/ariation of dead load moment Si'.t is chosen as 0.10, 

Partial safety factors Pp and depend on p, the coefficient 

of variation of live load moment and the ratio of the mean 

values of the live to dead load moments, Por P = 4*26, 

^M])“ and = O.3O, partial safety factors Pjj and Pj^ 

are selected from Pig, 5,2 as 1.3 and 1,8 respectively, 

2he coefficient of variation of moment capaci'fcy depends on 

coefficients of variation of strengths of masonry and steel 

and percentage of reinforcement, is expected . to be 

in the range of 0,05 to 0,10, Por P = 4*26, the strength 

reduction factor is found to vary from 0,68 to 0,84 as 

seen from Pig, 5»5 £or ~ 0,05 to 0.10, Therefore, is 

^ -^st 

selected as 0*7. In RBB, is usually less than 0.5 

percent. Por fyjj- = 415 U/mm^ and f^ in the range of 8 to 
15 IT/mm^ j ,£2^— varies in the raiige of 0.1 38 to 0,26, 

In this range, the partial safety factor for masonry strength 
y is found to be 1.75 to 2,5 for'y= 0,7 from Pig, 5#4» 

^ Ag-t- ^vk 

If -.1^ , = 0,2 is taken, is equal to 2,0 as can be 

bd fwk » w - 

seen from Pig, 5,4. BS 5628: Part 2 original draft code (48) 
has recommended partial safety factor for masonry equal to 
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2*5 to 2,8. The revised draft code BS 5628: Part 2 has 
now suggested the partial safety for masonry ^ as 
2.0 to 2,3 ( 134 ). A partial safety factor of 2 to 2.5 
for strength of masonry is recommended, when partial 
safety factors 1.3 and 1.8 are used for dead and live 
load moments (to be applied to the mean values) 
respectively. 
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DTJRIBILITI OF REIliFORCES BRICSC MiiSOHRI 
6.1 Introduction 

Durability of a structural system is defined as the 
ability to endure against the actions of destructive environment 
The destructive process could be of chemical, physical or 
mechanical origin, Oiemical factors are the most influenti€il 
in causing corrosion of building materials. Humidity with 
oxygen is an important factor for corrosion under normal 
building conditions. The G-ovemment of Bihar at the request 
of National Building Organisation carried out an extensive 
survey (117) covering aroui^ 250 houses and other buildings 
in the state, A considerable deterioration was observed in 
the structures, the age of which varied from 20 to 40 years*^ 
Similar deterioration wsis observed in the case of about 5000 
houses of age varying from 50 to 54 years in New Delhi (117)*' 
Some of the major defects and type of deterioration found in 
reinforced brickwork are given below (117,118); 

(i) Cracks were observed inmost of the buildings of 
age varying from 20 to 25 years. The damage was 
found to be less for the roofs which were protected- 
from direct exposure of sunlight. 
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cracks had appeared in most of the bxiildings 
after 20 to 25 years^ there were some instances where 
the plaster had been replaced just 12 years after 
construction. 

A network of cracks was observed in the bottom face 
of the roofs along the main reinforcements and along 
the distribution bars. 

The reinforcements were severely attacked and a 
considerable volume of partially oxidized material of 
dark colour was observed. In a good number of cases, 
it was noticed that the reinforcements of l/2** to 3/Q* * 
diameter were reduced to half their diameter. 

Dayaratnam (119} presented some results on durability case 
studies of reinforced brick masomy structures*' It was 
concluded that the corrosion of reinforcement appears to be 
the most important factor that governs the serviceability and 
strength of reinforced brickwork. Durability aspect of 
reinforced brickwork was discussed by Poster (3). Corrosion 
of reinforcement was found to be related to the d^th of cover 
and even 40 years old reinforcement bars were little affected 
provided they were not touching the bricks having a cover of 
about 50 mm of dense brick and mortar. Spalling of mortar in 
the joints and splitting of the bricks in horizontal planes 
at the level of reinforcanents were reported by Dayaratnam(1l9). 


(ii) 


(iii) 


(ir) 
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Durability estimation is rather difficult because 
many uncertain environmental factors enter into the probl«ii,i 
It is time consuming and a satisfactory conclusion m^ not 
be easily related to all possible circumstances. Rajagopalan 
et al, (120) reported results of an investigation into the 
causes of corrosion of reinforcements in reinforced concrete 
and reinforced brick constructions. Attenpts have been made 
by different researchers to study the durability aspects of 
reinforced brick masonry through possible short term measuraaents 
and observations, but it does not give any quantitative assessment 
of durability itself, (The electrochemical process of corrosion 
of steel ^ihedded in concrete was discussed by Bazant (121), 

The time to cracking or spalling of concrete due to corrosion 
was formulated by Bazant (122) and is given by 

"^cr “ ^p ^cor (6.1) 

where t^^ = time at which corrosion would produce cracking 

or spalling of the concrete cover 
t = time of depassivation 

ir 

^cor ~ duration of steady state corrosion. 

Duration of steady state corrosion t^Q^ was related to the 
rate of rusting and the increase in diameter of the reinforcement 
bar due to corrosion. Different possible causes of deteriora- 
tion of reinforced masonry with particular reference to 
corrosion of reinforcements are reviewed and discussed 



in the following sections. An attempt is made to estimate 
the time to craclcing or spalling of masoncy due to coonrosion 
of reinforcement following similar f oimtilation given by 
Bazant (122). finally, reliability analysis considering the 
decrease in moment capacity of reinforced brick beams Y/ith time 
for probalistic external moment is presented* Moment capacity 
of RBB at a particular time is assumed as deterministic 
function of initial random strength. 

6,2 Causes of Deterioration 

Deterioration of reinforced brick masonry structure 
can be broadly classified into three aspects (i) chemical 
attack on the mo 2 ?tar and brick, (ii) corrosion of embedded 
reinforcement and (iii) mechanical wear and tear due to 
repeated loads, 

6,2,1 Chemical attack 

Mortar present in brick masonry may be affected due 
to the attack by different aggressive media like acids, salts 
with ex'ohangable ions, Ihe degree of attack depends on 
the strength and concentration of the aggressive medium. 

The cement present in mortar is attacked more easily, The 
main components of ordinary portland cement are (123): 

(i) tricalcium silicate 3Ga0,Si02 

(ii) dicalcium silicate 2 Ga 0 ,Si 02 



215 


(iii) tricalci-um aluminate 3CaO.Al2*^5 

(iv) tetracalcium aluminofe2?rite 4Ca0.^2*^3*^®2^3 

(v) gypsum CaS0^.2H20 

Hie de'fceriora'tioiL of cement mortar due to chemical attack 
can be either dissolving attack or swelling attack. Strong 
mineral acids such as sulphuric acid, hydrochloric acid 
and nitric acid dissolve all components of set coaent whereas 
weak acids and many organic acids only form water soluble 
compounds present in the cement (124). The degree of acid 
attack is specified throu^ range of pH values in G-erman 
standard DIN 403O (124). pH value less than 4,5 has a very 
strong attack, between 4.5 to 5.5 a strong attack and between 
5.5 to 6,5 a weak attack* Galcium, aluminium and iron salts, 
and silica gel form when strong acids react with set cement (124). 
Some typical reactions of sulphuric acid with different cement 
compounds are: 

3CaD,2Si02.3H20+3H2S0 ^ 30aS0^+ 2(Si02.3H20) (6.2) 

3Ca0.4l20^.6H20+6H2S0^- 3(GaS0^.2H20) + il2(SO^)._5 

+ 6H2O (6,3) 

3Ca0..I'e2O5*6H2O+6H2S0^ - 3(GaS0^.2H20) + He2(S0^)^ 

+ 6H2O (6.4) 

The typical reaction of hydrochloric acid on calcium silicate 

hydrate is; 

3Ga0.2Si02.3H20+6HCl 3GaGl2+2Si02-t6H20 


(6.5) 
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Carbonic acid is capable of dissolving lime from cement 
mortar or concrete and the reaction occurs partly as free 
carbonic acid and partly from the maintenance of carbondioxide- 
limestone equilibrium (124). Calcium hydroxide reacts with 
carbondioxide and forms calcium carbonate (CaCO^) which 
further reacts with GO 2 ijfi presence of water and foms 
csilcium bicarbonate; 

Ca(0H)2 + CO 2 ^ CaCO^ + E^O (6,6) 

GaCO^ + H 2 O + GO 2 ^ Ga(HG0^)2 (6.7) 

[Ohis process is known as carbonation, GaCO^ is sparingly 
soluble whereas Ga(HC 05)2 is hi^ly soluble in water. Thus 
it is capable of dissolving lime from cement mortar. Due to 
carbonation, pH value of coaent mortar or concrete reduces,' 

Sparingly soluble and voluminous reaction products may 
form inside the set cement during corrosive reactions and can 
exert pressure upon the surroundings. Calcium ‘aluminate 
present in portland cement when attacked by sulphate ions can 
form sulpho aluminate hydrate (5CaO,Al2®5*3CaSO^, 32H2®^» silso 
called ettringite (124,125)* 

3GaO.Al205 + 3(GaS0^.2H20) + 26 H 2 O - 3Ca0.4l205,3CaS0^.32H20 

( 6 , 8 ) 

Magnesium sulphate may also react with calcium aluminate 
hydrate to form ettringite (126); 
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2(3Ca0.il20^.12H20) + 3(MgS0^.7H20) 

- 3Ca0..il205.3GaS0^.31H20 + aAlCOH)^ + 3?fig(0H)2 + 8 H 2 O 

(6.9) 

Ettringite ( 3 Ca 0 „il 20 ^* 3 CaS 0 ^. 30 - 32 H 20 ) crystalizes in the 
fom of needle shaped crystal and has a "very large Tolume 
compared to the components it replaces because of the 
introduction of huge quantity of water of crystalli 2 ation(l 255 l 26 ). 
As the space it has to occupy is limited, the voluminous 
compound thus formed subjects to its surroundings an enormous 
pressure and therefore produces cracking and spalling. 

Evidences of sulphate attack on brick wall is reported by 
Benningfield (127). If sulphates are present in the aqueous 
phase for prolonged periods or arise from ground water in 
situations below damp proof course, sulphate attack is possible. 
Eormal building bricks of low permeability and suction rate 
offer good resistance against aggressive medium and are 
absolutely resistant to chomical action from the atmosphere (1 24). 
Bricks must be free of swelling components which may cause 
flaking off or destructive leaching action. 

6,2.2 Corrosion of reinforcement 

Reinforcement is normally placed in mortar joints of 
reinforced brickwork. The environment surrounding the 
reinforcement changes with time. Before the construction, the 
steel bars are exposed to atmospheric rusting which is due to 
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the simultaneous presence of water and oxygen (air)-. At- the 
time of construction, the bars are siorrounded by freshly 
mixed cement mortar which is normally so alkaline that it 
prevents further corrosion of embedded steel. The protective 
role of the cement mortar depends on its ability to maintain 
an alkaline environment round the reinforcing steel and its 
impermeability to oxygen(air), moisture and other corrosive 
chemicals. 

Corrosion of reinforcement is an electrochemical 
process. Olie likelihood of corrosion for a metal Immersed 
in an electrolyte solution depends on the pH value of the 
solution and the electrical redox potential of the metal (128). 
Fig, 6..1 shows the simplified Pourbaix diagram of iron (129). 
This shows the equilibrium regions where the metal is in a 
state of immunity, passivity or where corrosion will occur. 

The Pourbedx diagram for carbon steel is essentially the 
same as that of iron shown in Fig. 6,1, In the pH inteirval 
of 9.5 to 12.5, as can be seen from Fig,^ 6.1, steel remair^ 
in a passive state. A layer of ferrous hydroxide forms on 
the metal surface which protects the steel from, further 
corrosion. Hydration of Portland cement produces calcium 
a nd other hydroxides which causes the final alkalinity in 
excess of pH 12(124,128,130) and offers a natural protection 
against corrosion in cement mortar or concrete. Typical 
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reactions of hydration of portland cement are (128) 

2(3GaO.Si02)+6H20 -* 5Ca0.Si02.3H20+ 3Ca(0H)2 (6.10) 

2(2Ca0,Si02)+4H20 ^ 3Ga0.2Si02.3H20+Ca(0H)2 (6.11) 

Generally alkalinity of pore fluid of fresh caaent mortar or 
concrete lies within pH 12 to 13 (I30). Carbondioxide . 
present in air dissolves in the moisture present in mortar 
to form weak carbonic acid #iich over a long period of time 
slowly reacts with the alkalis in the pore water of mortar. 
Thus, the pH may come down below 9.5 which no longer ensures 
passivity against corrosion. At constant temperature and 
average humidity the progress of carbonation is app2?oximately 
proportional to the square root of time (124). 

Simultaneous presence of water and oxygen is the main 
cause of corrosion of reinforcement, Normally anode areas 
develop on steel surface and iron ions pass into the solution 
at the anode surface after the depassivation of steel (t21, 1 3I ) 

Pe -* Pe^^ +2 6” (at anode) (6.12) 

Adjacent to anode acts as cathode and dissolved oxygen in 
pore water reacts with incoming electrons to form hydroxyl 
ions in presence of water 

O2 + + Ae" -*• 4 oh” (at cathode) (6,13) 

Electric current is created in steel due to liberation of 
electrons. The negative hydroxyl ions liberated at cathode 
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aJ-rive at anodic area and electrically neutralizes ttie 
dissolved Ee forming a solution of ferrous hydroxide*^ 

S' e*"*" + 2 OH -♦ 5'eC0H)2 at anode (6.I4) 

The total cell reaction is 

21e + aHgO + O2 - 2fe(0H)2 (6.15) 

Herrous hydroxide is unstable in oxygenated solution(l51 ) 
and further reacts with available oxygen and water to form 
ferric hydroxide (Pe(OH)^) which precipitates from the 
solution, 

4J'e(0H)2 + 2H2O + 02^ 45*6 (OH)^ (6,16) 

This final product constitutes hydrated red rust (5'eO(OH.)+H20) 
and has a voliune of approximately four times thart: of steel(l2l). 
There are other forms of rust, for example Hl'e0(0H),HF‘e02j 
HeSO^ and 5'e^0^(black rust). The volume of black rust is 
approximately twice the volume of steel. The typical reaction 
is 

3Pe + 8 0H“ - He^O^ + 8 e*” + 4H2O (6.17) 

Generally, the tensile reinforcanents in reinforced brick 
beams and slabs are subjected to tensile stresses under 
service condition. Thus the danger of corrosion increases, 
because of several flexural cracks vhich allows moisture and 
oxygen to reach the reinforcements more easily. 
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6 * 2*5 Mechanical wear and tear due to repeated loads 

Effect of repeated loads on BB beams were studied 
by Srivastava (1 , 132)* Ihe increase in deflection due to 
repeated loads was rapid upto about 0*2 mega cycle and 
stabilized after about 0*5 mega cycle* Cumulative residual 
deflection of the order of I/13OO to I/900 of the span was 
observed under repeated loads* Experiments on RB slabs 
under repeated loads (2) showed that the ratio of recovery 
and the ratio of additional deflection to cumulative 
deflection damp out with more number of cycles which show 
ductile nature of RB slabs. The upper limit of repeated 
load was kept at about 80 percent of the ultimate load 
capacity of the beams in both experimental investigation(l,2jl32)* 
Due to limited experiments the decrease in ultimate strength 
of RB beams and slabs could not be related with magnitude and 
range of repeated loads and number of cycles. However , 
cracks in brick and mortar joints were observed after 
application of repeated loads* Propagation of cracks through 
mortar joints were also observed in some of the beams*- 

6*3 Cracking of Masonry due to Reinforcement Corrosion 

A typical RBB section is shown in Pig* 6, 2 (a)* Det D be. 
the original diameter of the reinforcing bar* Since the volume 
of rust is more than the volume of steel which is 
converted to rust, the average diameter increases and as a 
result exerts pressure to the surrounding brickwork , liet 
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mass of steel converted to rust per tinit length of 
the reinforcement and be the mass of rust thus produced.^ 

The increase in volume per unit length A Y can be written as 


A V 

where 


_ ^rust 
Prust 

^ rust 
^st 


(Y 1) M 

Pst ^ ^st 


density of rust 

density of steel =7.85 gm/cm- 
%ust ^ st 

Prust 


( 6 , 18 ) 


M, 


st 


( 6 . 19 ) 


~ volume expansion ratio which depends on the 
type of compound produced. 

Let aL be the average increase in diameter and can be 
found from Eq, 6,18 


AY 


= -^‘[(D + AD)^ - D^] = (Y^-1) . ^ (6.20) 


Assuming Ap << j)^ 311,3^ neglecting' higher orders of Ah, . 

Eq. 6.20 can be written as 

- 2 M3-{; 

As = (V„ - 1) . , .{6.21) 

^ Its ‘"st 

Eq, 6.21 gives the relation between AD and the mass of 
steel Mg3j which is converted to rust. Let p^ be the 
radial pressure developed on the surrounding brickwork due to 
increase in diameter a B as shown in Eig. 6, 2(b), Two 
possible modes of cracks can develop due to the corrosion 
of reinforcement as shown in Eigs. 6.2(c) and 6, 2(d). 

If spacing of the bars (s) is large compared to the diameter 
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of bax, failure mode consisting of 45° inclined cracks can 

be assumed as shown in Fig. 6.2(c). let f^ be the average 

tensile strength of brickwork. Considering the equilibrium 

of vertical forces, can be related to f^ as 

D ^ 

“ = 2 fwt ~ ^ 

vhere d^ = clear cover to the reinforcement as shown in 
Fig. 6.2(c) 

Inclined cracks of 45° are found mostly in reinforced . 
concrete constructions. If the spacing of bars(s) is not 
large compared to D and parallel bsirs rust simultaneously, 
parallel cracks as shown in Fig. 6, 2(d) are expected which 
in turn may cause cover peeling. Splitting of the bricks 
in horizontal plane at the level of reinforcement at some 
places were observed (119). Equating the vertical forces, p^, 
can be related to f^ as 

Pj, B = (s - B) (6.23) 

Although average tensile strength of brick work in parallel 

cracks will not be in general equal to the average tensile 

strength in inclined cracks, no distinction is made in the 

present formulation. Ihe condition that a parallel crack 

will form is given by 
s - 2D 

^c > 


2 


(6.24) 
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Brickwork; is not a homogeneous and isotropic rnediim# 

It is assumed here that the surrounding brickwork is homogeneous 
and isotropic medium for the following derivation. Assuming 
a thick cylinder of thickness d- as shown in Fig, 6, 2(b), AD 
can be related to the internal pressure p^, as 


ad - p^ (6.25) 

where is the flexibility coefficient which depends on 
Poisson's ratio. Young's modulus and other parameters. 
Under the assumption of plane strain, a™ is given by 

{ 1 - B ^ ^ 

= [1 + ]•— (6.26) 


2do(a„+D ) 


where 


V 


E 

d. 


= Poisson's ratio 
= Young's modulus of brickwork 
= Thickness of assumed cylinder, i,e,^ clear 
cover to the reinforcement. 

Assumption of plane strain is justified when the length of 
the member is large compared to the other dimensions and 
corrosion occurs in a considerable length. If there are 
cracks perpendicular to the reinforcements (as may be the 
case in KB beams at bending), plane stress assumption is 
justified. Under the assumption of plane stress, is 

given by (122) . 


a^ - [1+ V + 


D 


2d„(do+ B) 


-]• 


D 

B 


(6,27) 



227 


If a cylinder of infinite space (d^ °°) is asstimed, 

(1+v)D 

~ "g (6,28) 

Eq., 6,28 is valid under both plane strain and plane stress 
assumptions, Thus, coefficient a^ in Eq_, 6,25 is bounded as 

< a^’ 

where a^ is given by Eq. 6,28 and aj,' is given by Eqs. 6,26,. 
or 6,27 for plane strain or plane stress cases respectively. 
For practical purpose, the flexibility coefficient may be 
taken as 

« (a^ + a'^ ) / 2 (6.29) 

When parallel bars rust, the effect of adjacent bars can 
be taken into account and an approximate formula is given 
by Bazant (122), For s > 10 D, the effect due to adjacent 
bars on Uj, is around 2 percent and can be neglected. 
Substituting value of from Eq, 6,25, Eq, 6,22 and Eq, 6,23 
become 

D 

ad = 2 f^ ^^C 2 ^ • d” (inclined crack) 

= f^ (s-D) . (parallel crack) 

“ (6.30) 

Eq. 6.30 gives the increase in diameter AD required to cause 
cracking. The mass of steel Mgt which has to be converted 
to rust to cause cracking or spalling of brickwork is given 
by (from Eq. 6.21 ) 
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“at’' j 

(Tj,- 1) ' 



D 

+ ^ (for inclined crack) 


2(y3,-1) 


f^(s-D). (for parallel crack) 


(6*31 ) 

^sunning tliat the steady state corrosion begins after the 

depassivation time tp, the time to cracking due to 

corrosion is given by Eq.. 6.1, where t^^^ = 

and Dqqj, is rate at v»hich the steel is converted to 

rust. 


The above formulation is based on the assumption 
that there is no void or clearence between the bar and 
surrounding brickwork, i.e., a slight increase in diameter 
will start developing pressure to the surrounding medium. 

In a practical situation, a thin flexible film surrounding 
the reinforcement or any other t 3 rpe of void may exist and in 
that case t^j, can be estimated by modifying Bq., 6,1 as 


'cr 


where 



^oor2“ 


^p ^cor1 ^cor2 


(6.32) 


depassivation time 

time taken to compress the flexible film before 
exerting pressure on the surrounding brickwork 
time taken to cause cracking after the flexible 
film is compressed 

Mg^/jgor wkere is to be computed from Eq..6.31. 
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Estimation of Yr^ 

Yolxme expansion ratio given by Eq.. 6.19* depends 
on the type of rust compound produced, Y^ calculated by 
Eg., ,6.19 for some typical compottnds are shown in lable 6.1. 
Densities of different compounds shown in Table 6»1 are 
taken from Handbook (133)« It should be noted that one 
atom of Ee produces one molecule of Fe( 0 H) 2 ,Ee( 0 H)^ 
and EeO(OH) whereas three atoms of Ee produce one 
molecule of Ee^O^, The volume expansion ratio is 
approximately 2 for black nnist whereas Y^ varies from 3*85 to 
4*42 for red mst depending on the extent of hydration as 
can be seen from Table 6*1, Thus production of red rust is 
more dangerous than black mst for causing crack to the 
brickwork. However, the proportions of different compounds 
formed in actual situation is not known. 

Depassivation Time tp 

After the alkaline environment surrounding the 
reinforcement is lost, the steady state corrosion begins.- 
If the protective layer to the reinforcement is broken 
repeatedly by some other means (e.g., repeated drying and wet 
condition) or repeated loads, steady state corrosion begins. 
Depassivation time tp depends on the following factors: 
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Table 6,1 : Yolnme Expansion Ratio for 

Different Rust Compounds 


Rust 

Compound 

Density 

brust 

(gm/a2i3) 

Mol.weight 

of 

compound 

^st 

^‘^st 

P st Y 

— - 'ij, 

Prust 

Ee(0H)2 

3.4 

89.85 

1.609 

2.309 3.71 

Ee(OH)^ 

3.4“ 3.9-»^ 

106.85 

1.913 

2.013-2.309 3.85-4.42 

EeO(OH) 

4.28 

88.85 

1.591 

1 .834 2.918 


5.18 

231.55 

1.'382 

1.520 2.100 


■» Depends on the extent of hydration 

Note: At,, wt. of Ee = 55.85, Pgt ” gm/cm^ 
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(i) Rate of carbonation of cement mortar Yiiich. causes 
acidic environment, 

(ii) Atmosphere having high concentration of sulphur 
dioxide which also destroys the alkaline environments 

(iii) Sulphate content of bricks. Sulphates leached out 
from brick can destroy the alkaline envii^onment. 

(iv) If the reinforcement touches the brick, there is no 
alkaline environment which v/ill protect the reinforcement 
from corrosion, 

(v) Percent of water absorption of bricks. 

Quantitative assessment of depassivation time is very difficult 

because it depends on numerous factors as given above. 

Depassivation time can be increased as follows; 

(i) Using bricks of low permeability and low suction rate, 

(ii) Using dense mortar which will give low permeability, 

(iii) Careful supervision that reinforcements donot touch 
bricks, 

(iv) Increasing clear cover to the reinforcements. 

(v) Plastering of all exposed surfaces. Water proof 
plastering or protective coating in cases where the 
brickwork is exposed to high humid or moist condition. 

Estimation of t^Q^,-} 


let c’ be the thickness of flexible film around the 
reinforcement. Assuming all round rusting occurs, the time 
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taken to compress the flexible film before exerting pressure 
to th.e surrounding brickwork can be expressed as 



(6*33) 


where 



average penetration rate per unit time 
volume expansion ratio. 


Estimation of thickness of flexible film has considerable 
influence on the initial time taken before developing 

pressure on the surrounding brickwork. Estimation of 
would involve a very sophisticated instrumentation. The 
value of c ' at the most can be assumed as the permissible 
crack width of water retaining structures. It is only an 
rough estimate* c' 0,01 mm is assumed in the numerical 
illustration. 


Estimation of tQOp2 


After the flexible film is compressed, further 
corrosion of reinforcement will start developing pressure 
on the surrounding brickwork iidiich in turn will cause cracking. 


Assuming all round rusting as before, 


t 


cor2 


Mst _ ^st 
3r ^ 


(6.34) 


where = average penetration rate per unit time and is 
to be calculated from Eg,. 6.31* Eq.s.6.33 ^nd 6,34 are based 
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o n the assumptions that c • and are both small compared to 
the diameter (D) of the reinforcement, Application; of the 
above equations is illustrated by the following example*^ 

Young's modulus of masonry (E) is found to be 
between 400 to 1000 times the masonry crushing strength ( 134 )* 
Hendry ( 36 ) gives an approximate relationship 

E = 700 f^ 

where f^ is the crushing strength of masonry. 

Long term Young's modulus can be taken as one half to one third 
of the instantaneous Young's modulus (34,36). Poisson’s ratio 
of brickwork can vary from 0,11 to 0,20(34), The proportion 
of brick and mortar enclosed by the assumed cylinder (shown 
in Eig, 6.2(a)) is impoirfcant to select the appropriate value 
of E, If the area of mortar in the assumed cylinder is 
dominant as compared to the area of brick, it may be reasonable 
to take Young’s modulus and Poisson's ratio of mortar in the 
numerical calculation. However, these values will of the same 
order of magnitude to the values of masonry. Let, 

2 

D = 12 mm, d^ = 25 mm., s = 95 mm and f^ = 10 H/mm 

If 1:3 mortar is used, average tensile strength of mortar 

2 

f^=1.1 N/mm • Assuming red rust is produced, 

TO fw P 

Let, E = = 3500 H/mm'^ and v = 0,15 • 

£ 1 ^ = 35,5 mm > dg = 25 mm , inclined cracks are 
expected from the condition given by Eq., 6 , 24* Assuming 
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plane stress condition, Eq* 6.29 gives = 0.004076. 
therefore, Eq. 6.32 gives the amount of steel to be consumed 
to develop inclined cracks as 

Pst^ 12 

^st = (25+ -^) (0.004076) 


= 0,04633 

Eield exposure studies on concrete cubes having reinforcement 
at varying cover thickness have shown all round rusting at 
13 nmi cover thickness (135)» The average corrosion rate was 
found to be 0.002 to 0*006 mm/year at different places for 
concrete free of salt (135). Eet the corrosion rate (jj.) 
be 0,001 mm/year for the 25 mm cover and calculated 

from Eq. 6.34j 


^cor2 


0.04633 Pgt® 

0.001(12) Pg^15 


= 3.86 yrs. 


Assuming c' = 0.01 ram, is calculated by Eq, 6.33» 

0.01 


0.001 ( 4 - 1 ) 


= 3.33 yrs* 


Carbonation could reach 50 ram in less than 20, years with 
the types of mortar likely to be used in reinforced brickwork (136), 
Although there is uncertainity about the rate of penetration, 
let 15 years be the time to carbonize 25 mm cover thickness*^ 

Thus, 


15 yrs 
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azid time to cracking or spalling of brickwork is given by 
Bq.. 6,32, 

^cr “ 3*33 + 3*86 = 22.19 yrs» 

used in th.is example are correct to the order of 
magnitude but differs from situation to situation. Even if 
ttie existance of thin flexible film is ignored, a conservative 
estimate of t^^. would be 

^cr “ ^5 + 3*86 = 18,36 yrs. 

It is clear that after the depassivation occurs, the duration 
of steady state corrosion to cause cracking of brickwork is 
only 3*86 years which signifies that the depassivation 
time (tp) is the controlling factor. Thus, care has to be 
taken to increase the depassivation time to avoid cracking 
or spalling of brickwork. However, after the mortar is carbonized, 
simultaneous presence of v/ater and oxygen is necessary for 
rusting of reinforcement. Continuous wetting and drying 
can give more rusting. Some of the measures which could be 
taken to increase the depassivation time were discussed earlier, 

6*4 Reliability Analysis of RBB when Strength Deteriorates 
with Time 

Moment capacity of a RBB section is a randcm variable 
since it is a function of other random variables f ^ and fy 
as seen in Chapter 4.' moment capacity of the beam can 
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decrease with thae if the beam is exposed to corrosive 
environment. The decrease in moment capacity may be either 
due to corrosion of reinforcement or due to decrease in 
strength of masonry or both depending on the type of corrosive 
environment, 

Iiet M(t) ) . the moment capacity after time t ,be a 
deterministic function of initial (random) moment capacity Mj, 
and expressed as 

M(t) = g(t) (6.55) 

where = initial moment capacity at t = 0 

(random variable) 

g(t) = a decreasing positive function of time t, 

Bifferent forms of deterioration function g(t) are possible 

depending on the type of deterioration,. When the loads are 

applied either in equal interval or at prescribed instants, 

the life of a deteriorating member can be measured in teims 

of ntimber of load applications. let the moment capacity of 

the beam decreases with time as given in Eq., 6,35 and 

independent of the niimber of load applications. let 

be the probability distribution function of the ith applied 

e^rtemal moment Me(i) at prescribed instants t^v Bet the loads. 

at prescribed instants be independent of each other. Therefore, 
n 

[TT % (i) %^(yo) dyo is the probability that the 
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m^ber with initial moment capacity between y^ and y’o'^^yQ will 
survive first n applied loads(58), when y(ti) is the moment 
capacity of the same member at time tj_> i»e«, after i— 1 load 
application. %j,(yQ) is the probability density function of 
Mj,.' lollo-wing Sq, 6*35, y(tj^) is given by 

y(ti) = yo sC'ti) (6.36) 

Summing over all possible values of y^, the probability 
of the member, with resisting strength specified hy the 
probability density function fM_(y„)»to survive a series of n load 

JL w 

applications, i.e,, the reliability function is given by 

oo n 

%(“) = /[ H (y(tj^))] % (jo) ayo 

° i=1 

The probability of failure is given by 

oo ■ n 

F„(n) = 1 - %(n) = / [1- n % %p(yo)^yo 

0 i“1 ® . 

(6*38) 

' Eq. 6.38 can be written as 

OO ^ 

= /[I - 1 (1) %j.(yo^ *^0 

O i=1 (6.39) 

where 

Now, n (1 (i) (y('ti)) can be approximated by 

i=1 ® 

[1- ^ % (i)^^^’^i^^ 1 ^ written as 

issi , 8. ■ ■ ■ 
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%(>!) - / SMgd) (yrti)) fMj.(yo) ayo <6-40) 

The approximation is based on the assxmption that the ma 30 r 
contribution of the integral in Eq.* 6»39 comes from those 
values of y^for ^ch « 'I * 

It follows that 

%(ti) = %3.(yo) ay„ (6.41) 

Hence, Eq., 6.40 becomes 
n, 

Pjj(n) * ^ %Q(i)^y^'*'i)) ayCt^) 

= (6-+2) 

where p^ ("ti) is the probability of failure of the member 
resulting from ith application of external moment and is 
given by 

Pf(ti) = s\u) (y(ti)) %(t.)(y(^^i» ay(i=i> (6-+3) 

If all the external moments MgCi) japplied in a sequence , 
belong to one and the same distribution %g(x) * ^9.* 6.43 
becomes 

oo 

Pf(ti) = / [1 -%g(y(ti))] ■%(*!) (y(q» 

° (6.44) 
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Using 6«41 > ills above e<iua'tion can be wriiienr as 

oo 

= {[ ’- ( 6 -+ 5 ) 

P:f('fc^) is tile probability of failure of the member resulting 
from i oh application of external moment and can be expressed 
in an alternative form as 

^ (M(t^) < Mg) 

= P (\ < M^gCti)) 

OO 

where 

= probability density function of external 
moment Mg 

- probability distribution function of 
moment capacity 

The probability distribution function given in 

BQ.. 4.35* When strengths of masonry and steel follow 
normal distribution ■, Vfv, Pfo 

to be computed from Eqs. 4.40 and 4.42 respectively. 


% ^ 
%^(x) 
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The reliability function becomes 

%(“) = ■' - 1=^ (6.47) 

The application of the above equations is illustrated 
through the following example. 

Example 6«1 

Strength of masonry f^ and strength of steel f_ 
are distributed as 1T(8,.96, 1,26) l/mm^ and N(449. 15, 22.457 )2J/mm^ 
respectively. Ihe external bending moment on the section in 
every year is indopondent and lognomally distributed with mean 
7 klSfei and coefficient of variation 20 percent • Moment capacity 
of the EBB section (details of the section are given below) 
decreases with time to 75 % of its initial strength at the, 
end of 50 years period, Pind the probability that the section 
will fail in 40 year period, when: 

(a) moment capacity decreases asymptotically mth time 

(b ) moment capacity decreases parabolcally v/ith time 
having a zero gradient to start with, 

Details of the section are as follows: 

b = 550 mm , d = 175 mm 

= 255*62 mm^(5 Boa, 10 bars) 
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The external moment is given by: 

^em “ ^ and = 0.20 

Using the mean values of the strengths of the 
masonry and steel, the mean moment capaoity is computed 
from Eq. 4.11 and it is 

Mct = 16,415 kHin. 


The parameters of lognormal distribution are obtained from 


Eqs.4.44 and 4.45 as; 


M' 


'In 

In 


= 0.198 

= 6.864 


The ratio of the mean values of the moment capacity and 
external moment at the initial stage is 


^ 16.415 
M, 


2.5447. 


em 


(a) Asymptotic deterioration 

Tho deterioration function is assumed as 
g(t) = 

At t - 50 yrs and for g(t) = 0.75, the value of a works out 
to be a =0.00575. 

-0.00575 t 

Therefore, g(t; = e 

The graph of the above deterioration function is shown by 
curve A in Eig.6.5. The probability that the section will 
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fail in 40 year period is computed from Eq.,6.42 and it is 
given by 

. 40 

P(1T<40) «i; p^(t^) = 7.336x10*“^ 

i=1 

^f^'*^i^ computed from Eq, 6.46. Hie probability 

distribution function used in Eg., 6,44 is given in 

Chapter 4» The probability of failure at the end of first 
year is found to be 1,195x10 It can be noticed ISiat the 
probability of failure in forty year period has increased by 
6 14 times that of the first yeair’s probability of failure* 

Hie probability of failure of the section in forty year 
period in case the strength of the section remains constant 
with time (i,e,j no deterioration) will be about forty times 
that of the probability of failure in the first year. Hie 
twenty five percent deterioration has increased the probability 
of failure by maiiy times* 


Curves A, B and C in Eig, 6.3 represent 25% , 

15 % and 5 % decrease in strength in 50 years re^ectively. 

Hie probability of failure at every year and the cumulative 
probability of failure with number of years cori?esponding 
to these degrees of deteriorations (shown in Eig, 6,3) are 
shown in Eig. 6.4. When the degree of deterioration is 
increased from 5% to 15%, "tbe cumulative probability of • 
failure at 51st year increased from 9*3182:10 to 4*008x10 ^i«e*» 

3 times increase in degree of deterioration has increased 
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the probability of failure by 4.3 times., cumulative 

probability at 51st year is 2,156x10"^(i.e. I8O4 times 
increase with, respect to the 1st year*s probability of 
failure) when the degree of deterioration is 25 %as can be 
seen from Pig, 6,4. As the degree of deterioration is 
increased by 5 times (i,e,, 5 %to 25 % ),the cumulative 
parab ability of failure has increased by 23 times* 

(b) Parabolic deterioration 

Asstmie the deterioration function to be 

S C t ) “ a^ "t" a^t "i" a^t 

d 

At t = 0 , g(0) = 1 and ■— g(t) ] =0 

'at t = 0 

give a^ = 1 and a2 = 0, Hence, the deterioration function 
becomes 

g(t) l+a^t^ 

At t = 50 yrs and for g(t) =0,75 , the value of a^ works 
out to be a^ = - 0,0001, 

[Cher ef ore, g(t) = 1 - 0,0001 t^ 

This type of deterioration can be expressed in a general 
form as 

g(t) = 1 - at^ 

The graph of the above deterioration function is shown by 
curve A (25 % decrease) in Pig, 6,5. ®ie probability that 
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Numblf of years 
Figs. 64 ' of failure with number of years 
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■tlLe section will fail in 40 years period is computed fisjm 

Eq., 6*42 and it is given by 

40 

P(ir<40) -a 2 Pf(t.) = 2.233x10“^ 
i=1 

IThe probability of failpre at tbe end of first yean is found 
to be 1. 058x1 0“5. The probability of failure in 40 years 
period has increased by 211 times that of the first yearns 
probability of failure. 

Curves A, B and G in Pig. 6,5 represent 25%, 15% and 
5 % decrease in strength in 50 years respectively. Ihe 
probability of failure at every year and the cumulative 
probability of failure with number of years coinre spending 
to these degrees of deteriorations (shown in Pig, 6,5) are 
shown in Pig, 6,6, TlVhen the degree of deterioration is 
increased from 5% to 15%, the cumulative probability of 
failure at 51st year increased from 8,022x10“^ to 2,671x10'“^, 
i.e., 3 times increase in the degree of deterioration has 
increased the probability of failure by 3*5 times, The 

—2 

cumulative probability of failure at 51st year is 1,287x10 
(i,e,| 1216 times increase with respect to the 1st yearns 
probability of • failurG)when the degree of deterioration is 
25 % . As the degree of deterioration is increased by 5 times 
(i.e.', 5 % to 25 % ), the probability of failure has increased 
by 16 times. It is to be noted that in case of assrmptotic 
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deterioration, the probability of failure increases by 23 
times whereas for parabolic deterioration it is 16 times* 

Thus, the type of deterioration function and the degree of 
deterioration have a direct impact on the cumulative proba- 
bility of failure. 

6.S Discussions and Conclusions 

The alkaline envi 2 X>nnient caused by the hydration 
of cement mortar gives a natural protection to the reinforce- 
ment. Steel remains protected and immune to corrosion as 
long as the pH value of the surrounding medium is around 
10 to 13*' Due to action of salts present in the brick 
or due to carbonation, the pH value of mortar decreases and 
consequently the danger of corrosion becomes inevitable. The 
volume of rust is two to four times that of the corresponding 
amount of steel reacted depending on the type of rust compounds 
formed. The increased volume of the material due to rusting 
exerts pressure on the suirroiuiding brickwork causing cracking 
and spalling of brickwork. Spalling of mortar and splitting 
of bricks were observed in some cases (119). Estimation of 
time of cracking or spalling of brickv/ork due to corrosion 
of reinforcement is formulated* The time of cracking depends 
on depassivation time, the diameter and spacing of reinforce- 
ment, Poisson^s ratio and Young’s modulus of brickwork, average 
tensile strength of brickwork, type of rust compound produced 
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ax^ the rate of corrosion. Simultaneous presence of water 
and oxygen is necessary for corrosion of reinforconent after 
the mortar is carbonized,: Once the steady state corrosion 
starts after the depassivation period, the time required for 
cracking or spalling of brickwork is estimated to be about 
three to four years for a corrosion rate of 0,001 mm/yr, lEhe 
depassivation time which is the period before steady state 
corrosion stairbs, is the dominating factor that governs the 
total time elapsed before cracking or spalling of brickwork 
takes place. Thus, care has to be taken to control the 
depassivation time by the measures discussed earlier to 
increase the durability of reinforced brickwork. 

Due to corrosion, the effective area of reinforcement 
decreases which results in decrease of moment capacity. 
Deterioration in bond between reinforcement and brickwork due 
to rusting is not yet understood. The moment capacity of a 
RBB section can also decrease with time due to decrease in 
brickwork strength by chemical reaction,; Moment capacity of 
the section at any time is assumed to be a decreasing positive 
function of the initial (random) moment capacity but 
independent of the load applications. The reliability analysis 
considering the decrease of moment capacity with time is 
formulated when there are series of independent load 
applications* The probability of failure can be estimated by 
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sinmning the probabilities of failure at each load application 
considering the deteriorated strength of the section, !Pwo 
types of deterioration models namely asymptotic and parabolic 
deteriorations are considered in the illustrative example. 
However, the actual type of deterioration will depend on 
the ts^pe of exposure condition* Cumulative probability 
at 5"Ist year is used for comparison purpose, Por 25%deg3?ee 
of deterioration, the cumulative probability of failure has 
increased by I8O4 times and 1216 times with respect to first 
yearns probability of failure for asymptotic and parabolic 
deteriorations respectively, for a 5 times increase in 
degree of deterioration (from 5 % to 25 % ), cumulative 
probabilities of failure are found to increase by 23 times 
and 16 times for asymptotic and parabolic deteriorations 
respectively, Ihus, the type of deterioration function and 
the degree of deterioration have direct impact on the 
cumulative probability of failure, 

Ihe probability of failure at each year with 
different ratio is plotted in fig. 6.7 for both 

asymptotic and parabolic deteriorations. It can be seen 
from fig.i 6.7 that the probability at eveiy year depends on 
the type of deterioration model and decreases with 
ratio,. Ihe emulative probability of failure with number 
of years for different Mrn/^em ^ fig,. 6.8 




^ 5 deer fast !n 50 yr« 
A sym ptotic 
Parabolic 
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for both types of dsterloratiors. As ratio increases 

the probahillty of survival also increases as e:.peoted. 

Mm/’V is ills initial load factor applied to the mean 

values. ?or = 2,5^ the oum\a.atlve probabilities of 

failure at 51st year are 7*288x10"^ 4*4.25x10”^ for 

asymptotic and parabolic deteriorations respect i-vely. 

Similarly , for = 3,0, the cumulative probabilities 

of failure are 2*350x10 ^ and 1*415x10*’^ for asymptotic and 
parabolic deteriorations respectively. Hence, an increase in 

1*2 times has decreased the probabilities 
of failure by 30#6 and 31 times for asymptotic and parabolic 
deteriorations respectively*^ A small increase tt> 
ratio has decreased the cumulative probability of failure 
considerably. If ratio is increased by 1,4 times 

(i.e,, from 2*5 to 3.5), the probabilities of failure reduces 
from 7.288x10“^ to 1 *01 7x1 0*^ and from 4.425x10*"^ to 6.404x10*^ 
for asymptotic and parabolic deteriorations respectively as 
can be seen from Pig* 6,8, Thus, selection of initial load 
factor (i*e,, ratio) has got a direct impact on the 

probability of failure*. The probability of failure is also 
sensitive to the type and degree of deterioration. In a design 
problem, the probability of failure for a given design life 
has to bo checlced, Por a preassigned reliability j the initial 
load factor has to be increased depending on the t3rpe and 
degree of deterioration.! 



CHiPiEBR 7 

GONGIiCrSIOHS AHD SECCMMiSlDATIOI© 

7.1 General 

Statistical analysis of different properties of 
brick, mortar and masonry is carried out and presented 
in Chapter 2. Chi-square test is conducted on different 
random variables to fit suitable probability distributions* 
Reliability of chi— sqpare test is discussed presented in 
Chapter 5. Reliability analysis and reliability based 
design procedure incorporating statistical variations in 
various parameters are presented in Chapters 4 5 

respectively. Chapter 6 presents a study on the durability 
aspect of reinforced brickwork, in particular, corrosion 
of reinforcement. The conclusions and recommendations 
based on the work carried out in this thesis are summarized 
in the present chapter. Suggestions for further research 
are indicated at the end. 

7*2 Conclusions and Recommendations 

Nine different brick samples each containing about 
hundred bricks were studied for statistical variations in 
their properties. The first eight samples were taken from 
eight different manufacturers at random while the ninth 
contains several sets supplied by different building 
supervisors for quality analysis. The ninth sample can be 
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considered completely random* Seven different properties 
namely length, breadth, height, dry density, wet density, 
percentage of water absorption and flatwise compressive 
strength of each brick were measured and statistical 
analysis was carried out. Chi-square test was applied to 
test whether the data fits a particular probability 
distribution. 

The coefficient of variation in length and breadth 
of brick was found to lie within 2.8 percent indicating a 
high consistency in the dimensional accuracy of the 
made bricks* The coefficient of variation in length of 
brick from a single manufacturer was within 2 percent and 
that associated with many manufacturers was within 2*1 
percent* The variation in height of brick appears to be 
more than that of length and breadth. Even in this case, 
the coefficient of variation was limited to 5.4 percent* 
Dimensions of bricks manufactured from a single manufacturer 
showed a very small variability with a tendency towards 
deterministic value. The dimensions can be assumed to be 
normally distributed for practical purposes. 

Density of dry bricks of different individual 
manufacturers followed normal distribution* The mean 
values of dry density of bricks of different sanqsles 
differed reasonably but iaie coefficient of variation of 
each aemple was found to be within 5 percent. This again 
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reflects the dimensional accuracy of the hand moulded hricks 
and also the compaction and homogeneily in preparing the 
soil for brick making, The mean value of the percentage 
water absorption of brick samples varied from 13.41 to 
18.47 with the coefficient of variation ranging from 4.T to 
17.7 percent. Normal and lognormal distributions did not 
fit the water absorption data in five samples but in other 
three samples normal distribution can be accepted at 5 
percent significance level. In these samples, nc-pnal 
distribution appears to be a better model than lognormal 
distribution. Normal distribution can be accepted to 
represent the variability of percentage water absorption 
for practical purposes. 

Compressive strength of brick followed normal 
distribution at 2.5 percent significance level. Within the 
first ei^t samples, first three samples were collected 
from different manufacturers of Kanpur zone and the rest 
five were collected from Faroke in Kerala State. Mean 
compressive strength of brick samples from Kanpur zone 
varied frcan 14*57 N/mm to 25.69 N/mm with a coefficient 
of variation ranging from 22.5 to 24.6 percent whereas 
mean compressive strength of Paroke zone varied from 7.6 N/mm 
to 14.16 N/mm^ with a coefficient of variation ranging frcm 
15,1 to 18.6 percent. This shows that bricks of Kanpur zone 
are much stronger than those made in Paroke. Normal 
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distribution did not fit the random lot (the sample which 
contains several sets supplied by different building 
supervisors) at 1 percent significance level vdiereas 

lognomal distribution fits at 18 percent significance 
level* 

Cement— sand moirtar cubes of three different mixes 
(1 .5» 1 *4 and 1 :5) were cast in the laboratory and tested 
after 28 days, to study the variability of mortar strength. 

The coefficient of variation of the strength of mortar 
cubes in three different mixes varied from 10 to 18 percent. 
The coefficient of variation of strength of field cubes 
may be expected to be about 20 percent more than that 
observed in the laboratory, and could be anywhere between 
1 2 to 22 percent. This is comparable with the quality of 
Ml 5 concrete produced from nominal mix. Compressive 
strength of mortar (laboratory specimens) of each mixes was 
found to follow normal distribution at 10 percent significance 
level, Field data should be collected for a proper under- 
standing of the variability of mortar strength in actual 

practice. 

Thickness of mortar joint was measured to study the 
variability in joint thickness. Three different buildings 
of exposed brick walls were selected for study. Both bed 
joint and vertical joint thicknesses were measured at 
different locations of the wall selected at random,- The 
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coefficient of vari nv, -, 

d. bed joint thickness vanied 

between 12 to 14. -nAT-«av,4. -u 

4 percent whereas that of vertical joint 

tiHotoiess varied from 1 3 to 20 percent. 35ie loseet and 

highest measured thidmesees were observed to lie within 

+ 5 mm Of the mean which ^ 

wmcn signifies a very good quality 

control at the time of brick laying. Ihickness of mortar 

joint followed normal distribution at 1 peroent slsnifioanoe 
level . 


The strength of brick masonry depends on the 
quality of brick, mortar, joint thickness and joint layout. 
Due to non availability of field data, the strength of 
masonry was generated by Monte Carlo simulation on a 
digital computer from the randomly generated data of brick 
and mortar strength using their individual statistical 
properties. Strength of brick and moirbar Y/ere taken as 
normally distributed. Tyto deterministic formulae relating 
strength of brick, mortar and masonry were used for 
simulati^^n, and given by 

^w ^ ^^b ^m^ • ^ ) 

f^ = A(r..„758 + 0.155 f^ + 0.0082 f^fj^) (7.2) 
where f^j, fjjj and f^ are strengths of brick, mortar and 
masonry respectively, Eq, 7.2 was derived from the fommila 
given in SOEI recommendation (15) whereas Eq, 7.1 is based 
on limited experimental investigation (2). 
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Siniult, i;0(3. basad on "Pn 7 i * • -’ 

n Eq, 7,1 jjQ-j follov/ noimal. 

distribution even at 0 1 a. 

• P icenir signi^icancs level, Tiie 

central region oS the simulate! data matched vdth nomal 
distribution but mismatching was obseiwed at either tall. 
Simulated data based on 3q. 7.2 iollowed normal distribution 
at 0.5 percent significance level. Strength oi masonry may 
still be assumed as normally distributed random variable. 

Ohi-square test and Kolmogorov - Smirnov test (K-S 
test) 8X0 the two popular tests used to examine whether a 
data of a psxticular random variable follows an assumed 
probability distrioution. Usual ohi-square test with equal 
class interval for testing the suitability of a hypothesized 
prob<a,bili ...y dist.i’ibubion depends on three arbitrary decisions 
such ao ohe choict' oi the starting point, number of classes and 
width of classes. Hence this test can lead to doubtful 


accept unco o.c rejection of the null hypothesis, -fha three 
arbitrary choices in usual chi-square test can be made to 
a single varii'.blc if the test is conducted with dlasses of 
equal probability, Even if the classes of equal probability is 
used, %c computed value of chi-square statistic is observed 
to chtai;.^e rciidoialy with the choice of number of classes (k) 
and const. quoni;ly the maximum significance level (p-lcvel) also 
chgaiteSf Hence for a dcpcndablo decision, a series off 

tests with different number of classes is recommended. 


chi-squoxe 
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Based on the different p^levels obtained by the series of 
tests, characteristic p-level may be computed by Ea. 5.15 

and a decision can be taken for acceptance of the null 
hypothesis. 

Selection of a particular probability distribution 
requires engineering or scientific judgement considering the 
actual physical phenomena. It is possible that two or more 
probability distributions may appear equally competent 
models for a random variable considering the actual phenomena.i 
To test vdiich one of the probability distributions is a better 
model to represent the random variable, the characteristic 
p— levels for the competent distributions can be compared. 

The distribution which gives a higher characteristic p— level 
may be accepted as a better model as compared to the other 
distributions.' 

Kolmogorov-Smimov test (K-S test) does not require 
grouping of the data into classes as in chi-square test and 
hence arbitrary choice of number of classes etc.' are not 
involved.' Vl/hen the parameters of the null distribution are 
completely specified, K-S statistic is independent of the 
null distribution.' But the K-S statistic becomes dependent 
on the choice of the ndll distribution when the parameters 
are not specified and have to be estimated from the data.- 
In this case, standard tables available for use with K-S 
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test cannot be used, Sables showing the critical values 
Of K-S statistic for different level of significance are 
presented for uniform, nomal and exponential distributions. 

Ultimate moment capacity of a reinforced brick 
beam (BBB) section is a function of geometric proportions 
of the section, area of reinforcement, and. strengths of 
the masonry and reinforconent. She section can fail as an 
under-reinforced or as an over— reinforced. Due to random 
variations in the strengths of masonry and steely there 
exists a probability that a RBB section will fail as an 
over-reinforced even though the section is designed as an 
under— reinforced based on deterministic analysis* Probability 
of failure (p^) of a BBB section is the stm of (i) the 
probability of failure of the section as an under— reinforced 
(p^) Qhd (ii) the probability of failure of the section as 
an over-reinforced (Pfo)* 

The computation of failure probabilities and 
involves evaluation of multiple integrals and depends 
on the probability distribution of strengths of masonry 
and steel. The probabilities p^^ and p^^ are found to be 
sensitive to the relative magnitudes of the coefficients 
of variation of strengths of masonry and steel.^ The 
probability of failure of the section increases with the 
increase in coefficients of variation of strengths of masonry 
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and steel for a given evt^r-non 

n external moment, if the area of 

3^061 in a, RBB seo+iAvi a 

section IS increased, the probability of 

failure as an 'undeivreinforced (n ) 

vPfn’ decreases and the 

It IS observed that p^, decreases +n = to- •4-- 

rf vi«Q-reases to a limiting value eijual 

to as the area of q+^oi • 

^ steel IS increased to a very large 

value# which signifieq +ha+ 4-v,^ v i. 

sxgniiies that the probability of failing as 

an under-reinforced is alr.ost e®al to aero for heavily 

over-relnfofced sections. Por probabilistic ertemal moment, 

the probability of failure Increases as the coefficient of 

variation of external' moment increases. 


A large number of values of moment capacity are 
generated by Monte Carlo simulation on a digital cbn5)uter 
from randomly generated data of strengths of masonry and 
steel conforming to their individual statistical properties. 
Histograms of the simulated samples showed a marked negative 
akewnesst’ As the reinforcement ratio (Ag^A>d) is increased 
towards the balanced section value, the distribution of Mj, 
becomes more and more skewed to the right,' For a highly 
undor*-reinforcGd section (determini 3 tip. fl.ny ), moment capacity 
Mj, is found to be normally distributed. 


!Bhie design of a RBB section involves calculation of 
steel area and geometrical proportion such that the section 
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can resist the external design bending moment,: In a 
probabilistic design, area of steel and geometry of the 
section are so proportioned that the probability of. 
moment capacity being less than external bending moment 
is within a preassigned probability of failtire. A method 
which involves solution of an integral equation to find 
the area of steel or geometrical proportions of the section 
for a prescribed propability of failure ov reliability is 
presented in Chapter 5» It is found that the coefficients of 
variation in strengths of masonry and steel have direct 
impact on the design, 

limit state design method, based on semi-probabilistic 
approach, incorporates statistical variations in some of 
the design variables throu^ partial safety factors applied 
to loads and strengths of materials* Partial safety factors, 
required for a target reliability, are found to be dependent 
on the relative magnitudes of the coefficients of variation 
associated with loads and strengths of materials, Por a 
practical range of coefficients of variation (as discussed 
in Chapter 5) and for an accepted probability of failu3:^e 
of 10*^, the partial safety factors applied to the mean 
values of dead and live loads are found to be 1,3 and 1,8 
respectively* Corresponding partial safety factor for 
strength of masonry is found to be 1.75 to 2.5, Partial 
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safety factor of 1.15 13 used for strength of steel.. 

Original draft code BS 5628s Part 2.(48) rec<miiends a 
value Of 2.5 to 2.8 for the partial safety factor applied 
to strength of masonry at limit state of strength in bending. 
The revised draft code (BS 5628 s Part 2) has now suggested 
a value of 2 to 2.5. A partial safety factor of 2 to 2.5 
for strength of masonry is recommended for limit state 
design of RB beam section in flexure, when partial safety 


factors 1.3 and 1,8 are used for dead and live loads (to 
be applied to the mean values) respectively,' 


Corrosion of reinforcement in reinforced brickwork 
construction appears to be a serious problem from 
serviceability point of view.' Rusting of reinforcement 
increases the volume and causes thrust to the surrounding 
brickwork resulting in cracking and spalling of brickwork, 33ie 
volume of rust is two to' four times of the corresponding 
volume of steel reacted depending on the t3rpe of rust 
compounds formed* Estimation of time of cracking or spalling 
of brickwork due to corrosion of reinforcement is formulated.' 
Tho cracking of brickwork depends on depassivation time, the 
diameter and spacing of reinforcement, cover to the reinforce- 
ment, Poisson's ratio and Young's modulus of brickwork, 
average tonsils strength of brickwork, type of rust compounds 
produced and rate of corrosion. The depassivation time isdiich 
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i,S p03?XO(i l^©f^OlC0 ei+»-k4. 

steady state corrosion starts, is the 

dominating factor that governed +«+ -, 

s the total time elapsed before 

cracking or spalling of brickwcyir -nv, 

^ uricKrwork, fh.us, care has to be 

taken to control the depassivation time,: 


corrosion, the effective area of reinforcement 
decreases which results into decrease in the moment capacity, 
fhe moment capacity of a RBB section can also decrease with 
time due to decrease in brickwork strength by the action of 
aggressive environment. Reliability analysis considering 
the decrease in moment capacity with time is formulated and 
presented in Chapter 6, (Two types of deterioration models 
namely asymptotic and parabolic deteriorations are considered,' 
However, the actual type of deterioration will depend on the 
typo of exposure,' As the moment capacity of a EB section 
decreases with time, the probability of failure at every 
yeax' due to external moment increases and as a consequence 
the cumulative probability of failure increases. Bor 25 percent 
decrease in moment capacity in 50 yrs, the cumulative 
probability of failure at 51st year is found to increase by 
1804 times that of the first year's probability of failure 
for aisymptotic deterioration. For parabolic deterioration 
the cumulative probability of failure increases by 1216 
times that of the first year's probability of failure,: It 
is found that the type of deterioration and the degree of 
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deterioration have a direct Impact on the emulative 
probability oi failure. Cumulative probability of failure 

at a particular year depends on the Initial load factor 
chosen. An increase of load factor from 2.5 to 3.0 has 
reduced the emulative probability of failure at 51 st year 
from 7.288x10 ^ to 2.358x10“'^ for asymptotic deterioration. 
Similar findings are observed for parabolic deterioration* 
Selection of initial load factor has got a direct impact on 
the probability of failure,; In a design problem, the 
initial load factor ...las to be increased depending on the 
tjrpe and degree of deterioration to achieve a preaissigned 
reliability.. 

7 #3* Suggestions for further Research 

Dimensional variability of EBB section was not 
taken into consideration for reliability analysis. 
Investigation of dimensional variability in actual practice 
is suggested. Effect of variability of joint thiciciess 
on the strength of SBB can also be investigated. Ultimate 
strain and Young’s modulus of masonry can be taken as random 
variable in the reliability analysis and data may be collected 
to establish mutual dependence of the above variables.; 

Extensive data need to be collected on the strengths 
of mortar and masonry from field, and the decrease in strength 
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with time in different 

rent ejtposure conditions, ia^mptotic and 

parabolic types of deterioratinr, m.;, 

ration models are considered in 

the present work, it io -v-, , 

It 13 possible that the actual deterioration 

data may fit into one of these models . Collection of this 
type Of data retires considerable time. National Organisations 
can take up a long tem programme for collection of data on 
the deterioration under different erposure conditions. 


Research may be directed on the luantltative 
estimation of strength deterioration with time under 
different e.tpoaure conditions. Effect of corrosion of 
reinforcement on the bond between reinforcement and masonry 
may bo invostigaied to study the exact mechanism of 
deterioration of moment capacity of RB beams. Strength 
deterioration due bo repeated loading is not considered 
in the present work. Reliability analysis and design 
of RBli under repeated load conditions when strength 
deteriorates with each load application along with deterio- 
ration with time due to aggressive environment could be 
an important direction of further research. 
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^ 13 a function Of several rardoi. variable Z„x,„ 

Siv*G2i by . ^ 

^ = g (x^, X2, 

By erpandirg the function around the mean values of the 

ra^om variable and taking erpeotations and variance on 

both sides neglecting the higher order terms, the 

approxmate mean and standard deviatinn k 

deviation can be computed 

as ( 96 ) 
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where the partial derivatives are to be calculated at the 
mean values of X,.. x^^ and s^^ are the me^ »rd standard 
deviation of the variable is the correlation 

coefficient of the random variable and Xj. If XJs are 
independent random variable, variance of Y becomes 

5 > ^ On 
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Introducing partial safety factors and for 
strenj^ths of masonry and steel respectively, the design 
moment can be expressed as 

“a = , a(i-o .59 ^ ) 
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Prom Eq., B-1 the area of steel Ag^ for an under-reinforced 
section can be computed from the following quadratic equation 
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design bending moment. 


fflao lowest root given by Eq, B~3 gives the required area 
of steel. Using = 2.0, Yy = 1.15, area of steel 
(refoCTOd as is computed and is given in fable 5,1. 
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